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MEETING OF 1951 JANUARY 12 
Professor W. M. Smart, President, in the Chair 
The election by the Council of the following Fellows was duly confirmed : — 
*Bertrand John Harris, Royal Observatory, Greenwich, London, S.E.10 
(proposed by A. Hunter); 
J. A. Jacobs, M.A., Ph.D., A.F.R.Ae.S., 10 Woodland Avenue, Windsor, 
Berkshire (proposed by W. H. McCrea); and 
Ewen Adair Whitaker, Royal Observatory, Greenwich, London, S.E.10 
(proposed by A. Hunter). 


Ninety-three presents were announced as having been received since the last 
meeting. 


The President announced that the Council had awarded the Gold Medal of 
the Society to Professor Antonie Pannekoek, late Director of the Astronomical 
Institute of the University of Amsterdam, for his work in astrophysics and on the 
structure of the galactic system. 


* Transferred from Junior Membership. 





AN ATTEMPT TO DETECT A GENERAL MAGNETIC FIELD OF 
THE SUN BY A SPECTROGRAPHIC METHOD, USING A 
LUMMER PLATE 


H. von Kliiber 


(Received 1950 November 9) 


Summary 


The general magnetic field of the Sun has been examined, using the 
McClean spectrograph of the Cambridge Observatories, crossed with a 
quartz Lummer plate of excellent quality which gives a resolving power 
of about 600,000. This equipment was combined with a polarimetric 
device, similar to that previously suggested by Wood, consisting of a thick 
quartz plate and polaroid. Using two carefully selected Fraunhofer lines, 
one can in this way measure with great accuracy possible small wave-length 
displacements caused by weak Zeeman effects, the measurements being 
strictly differential and not affected by possible local Doppler effects. This 
arrangement was used in the Summer of 1949 to obtain spectrograms from 
two points on the Sun’s central meridian at +-45° latitude. There were about 
60 plates with altogether about 1000 interference points to be measured. 

The m.s.e. for a single determination of the magnetic field strength is 
about +6 gauss, and for each set of two corresponding plates about +3 gauss. 
For each of the two points on the Sun’s meridian 15 pairs of plates have been 
measured. ‘The internal m.s.e. of all these measurements is -+-0-8 gauss 
for each hemisphere. 

The final discussion shows that at the time of the observations (1949 
August-September) any possible polar field strength was too small to be 
measured within a m.s.e. of +0°9 gauss. The scatter of the observations, 
shown in Fig. 3 (a) and Fig. 3 (6), is within the random error of measurement. 
A field of 1-2 gauss should have made itself evident beyond all doubt on 
these diagrams. 





The existence of a general magnetic field of the Sun, in the form of a dipole, 
was made highly probable in 1913 by G. E. Hale and his collaborators (x, 2, 3). 
Further extensive investigations of this problem were not published for a long 
time. Smaller observational series at the Mount Wilson Observatory during 
recent years (4, 5, 6, 7, 8, 9) make it, however, very doubtful whether a dipole 
field of the order of about 50 gauss at the poles, such as was first found by Hale, 
really exists at present. A still unpublished photographic-interferometric series 
by the present author, in collaboration with H. Miiller (obtained in Autumn 1944), 
with the tower telescope at Potsdam during the last sunspot-minimum, seemed 
to show that the existence of a dipole field of even a few gauss was improbable 
at that time (10). G. Thiessen (11, 12) in Spring 1945 derived from visual- 
interferometric measurements the old value of about 50 gauss for a general dipole 
field of the Sun, but later (13, 14) he was unable to confirm this by his more 
recent observations, which appear to be consistent at most with a field of a few 
gauss only. 

Knowledge of the structure and intensity of the general magnetic field of 
the Sun has become of increasing importance for wide fields of astrophysics, 
including radio-astronomy, and also in connection with new fundamental 
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conceptions concerning the connection between rotation and magnetic field 
strength, recently advanced by P. M. S. Blackett (15). All this seems to justify 
attempting once again a photographic determination of the Sun’s field, as 
outlined in the present paper. 

The basis of the present investigation is the reasonable assumption that the 
general magnetic field of the Sun is a dipole field, similar to that suggested by 
G. E. Hale, and to that actually existing in the Earth. Furthermore, it is assumed 
that the field also exists in those layers of the Sun in which the Fraunhofer lines 
originate. Only then is the application of the Zeeman effect possible, and our 
measurements therefore refer to these layers only. Special series of measurements 
would be required if the field were not of a dipole type, as some theoretical 
investigations suggest (16, 17). 

Another assumption is that the suspected dipole field approximately coincides 
with the axis of rotation of the Sun. In this case the line of sight of an observer 
on the Earth runs parallel to the lines of magnetic force at two points at latitudes 36 
on the central meridian of the Sun. Therefore, if a magnetic field exists, we expect 
a purely longitudinal Zeeman effect at these two points. On the other hand, the 
magnetic field strength H, of a homogeneous magnetized sphere is a function of 
latitude obeying the law H,=H,)V1+3 sin?¢, where H, is the field strength 
on the equator (18). F. H. Seares (2) has already shown that in consequence 
the most suitable points for observations on the Sun’s disk to establish the 
existence of a longitudinal Zeeman effect are situated at +45° latitude on the 
central meridian. ‘There we may expect a field of about ? of the polar field 
strength. All observations of the present investigation have therefore been 
restricted entirely to these two points on the solar disk. 

The difficulties and uncertainties in the measurement of the Sun’s general 
magnetic field are essentially due to the extraordinary smallness of the expected 
Zeeman effect. ‘The classical Lorentz-formula for the amount of the normal 
Zeeman-splitting for a field H in gauss gives 


AA =C.. HA*p with the constant C = 4-7 x 1075 


for the wave-length in centimetres. The values for the Landé-factor p are generally 
between © and 3. It is possible to use for field strength measurements some 
lines in the solar spectrum which belong to the simple splitting-type and have the 
large Lande-factor 2-5 and 3 (19, 20). Even under the most favourable 
circumstances a field strength of 25 gauss corresponds to a Zeeman effect of the 
order of only 0-001 A. Since, however, the expected field-strength may easily 
be much smaller still (<5 gauss), it is quite out of the question to measure the 
real line-splitting even for the best usable Fraunhofer lines. One of the most 
suitable lines is FeA5250-22, which even in the normal solar spectrum has 
already a half-width of the order of o-1a. Using the well-known device of 
suppressing by polarimetric means alternately one or the other of the two 
circularly polarized Zeeman components, produced by the longitudinal Zeeman 
effect, only quite minute displacements of the centre of gravity of the line are to 
be expected. Under favourable circumstances it appears possible (11, 13) that 
sufficiently sensitive methods will disclose changes in the position of the centre 
of gravity of less than ;4, of the half-width. In principle one can hope that a 
sufficient accumulation of measurements will disclose field-strengths of a few 


gauss down to about 1 gauss. To achieve such accuracy it appears essential to 
1* 
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have a spectroscopic resolving power of the order of at least 0-01 a. for visual 
wave-lengths. Interferometric methods are therefore indispensable. 

Instrumental equipment.—The instrumental equipment for such observations 
consists of a spectroscopic and a polarimetric arrangement. 

Spectrograph: The McClean solar spectrograph of the Cambridge 
Observatories was used for the present investigation (21). In this horizontal 
arrangement the Sun’s light enters a doublet telescope lens via a coelostat and an 
auxiliary mirror. This lens has an aperture of 12 inches and a focal length of 
60 feet, and gives an image of the Sun about 16cm. in diameter (I mm. =11 seconds 
of arc) on the spectrograph slit. An otherwise excellent set of glass prisms showed 
during a preliminary examination very strong and irregular polarization properties. 
They were therefore replaced by a plane grating of 600 lines per mm., and an 
area of g by 12cm. (here not fully utilized), which was used in auto-collimation. 
The spectrograph lens has a focal length of 450cm. The light of the most 
disturbing reflection images was eliminated by giving the lens a small tilt and 
by the use of suitable small diaphragms. The whole spectrograph was kept at 
constant temperature by an electric thermostat. 

The careful adjustment and testing of the grating spectrograph was carried 
out in the usual way with an Fe-arc and a cooled Hg-lamp. ‘The observations 
were made in the brighter second order, where the dispersion is 2 a./mm. ‘The 
resolving power of this arrangement was about 100,000. (We note that the 
grating was not fully illuminated.) It was, for instance, possible to separate 
easily and with certainty the close pair of Fe-lines A 4482-25 and A 4482-17, as 
produced in the ordinary Fe-arc. From the position and intensity of the ghosts 
of the grating one can estimate that they will not interfere seriously with the 
present measurements. 

After the experience previously gained by the author (10, 22, 23) with a 
similar arrangement, it appeared of interest to employ a Lummer plate as 
interferometer. Through the kindness of Professor S. Tolansky of the Royal 
Holloway College, London, a large quartz Lummer plate of excellent quality 
was obtained on loan. It was made by A. Hilger Ltd., London, and belongs 
to the Royal Society. The plate was mounted in an adjustable metal case and 
protected against temperature variations by a thick asbestos cover. It has a 
length of 200mm., a width of about 30 mm. and a thickness of 3-42mm. ‘This 
gives a range of 0-33 A. at 5250, which corresponds to an air gap of about 4:2 mm. 
in a normal Fabry—Perot interferometer. For the same wave-length the resolving 
power is about 600,000, which is one of the highest ever used for the solar 
spectrum. Inside the plate about 28 reflections with continuously decreasing 
intensities are produced, equivalent to about 20 reflections of equal strength. 

For the investigation of a continuous spectrum such an interferometer can 
only be used when crossed by an auxiliary spectrograph, in this case the McClean 
grating spectrograph just described. It is known that an interferometer can be 
combined with an auxiliary spectrograph in three different ways: externally 
mounted in front of the slit; behind the focal plane of the spectrograph; or 
internally mounted within the spectrograph itself (24). The last-mentioned 
arrangement cannot be used here because of the type of spectrograph construction 
and for other mechanical reasons. We may note in passing that with an ordinary 
spectrograph (not a Littrow mounting), a quartz Lummer plate should not be 
used between the dispersing element and the camera lens, because its anisotropy 
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then produces serious distortion of the interference fringes (25). To place the 
interferometer behind the image plane of the auxiliary spectrograph, which as a 
rule is one of the best places for it, had been attempted already during the 
Potsdam observations of 1943 by the author (10). This arrangement originally 
appeared particularly advisable as it was possible to place the whole interferometer 
as well as the polarization optics, both quite temperature sensitive, inside the 
perfectly temperature-controlled auxiliary spectrograph. However, the use of 
Wood’s procedure (see p. 7) then entailed bridging the rather large spectral 
range of about 120A., and also compelled the work to be done strictly 
differentially. ‘This arrangement was therefore abandoned in favour of the more 
orthodox method, using the interferometer in front of the slit. In the present 
paper this last arrangement was retained in order to achieve completely 
equivalent light-paths in the interferometer for the two wave-lengths to be 
compared (see p. 8). In the meantime, however, P. J. Treanor, S.J. (26) has 
made successful experiments by placing the interferometer behind the focal plane 
of the spectrograph. The author hopes to use an arrangement similar to 
Treanor’s for a future new series of observations of a close pair of lines, but for 
the present investigation the Lummer plate was placed in front of the slit of the 
spectrograph, in a section of the light beam made parallel by two auxiliary lenses. 
Both parts of the interference system were imaged on the slit with a high-quality 
lens of 25cm. focal length. A suitable field lens was used in front of the slit in 
order to concentrate also those light beams which have an inclination to the 
optical axis of the objective of the spectrograph (a procedure which is especially 
important in interferometric work) (22). 

P. J. Treanor has shown that the range of the interferometer should be at 
least three times the instrumental line-width of the auxiliary spectrograph if 
troublesome overlapping of false light from adjacent orders is to be avoided. 
This condition is fulfilled here. The range should also be related favourably to 
the intensity of the Fraunhofer lines investigated and was actually somewhat too 
large for the lines used. This is, however, not a serious disadvantage in what 
follows. 

The adjustment of the Lummer plate and the continuous control during the 
period of observations was performed by exposures on the hyperfine structure 
of the green Hg-line 45461, whose wave-length is near that of the observed 
Fraunhofer lines. An enlargement of such an exposure is shown in Plate 18). 
The light source was a Hg-lamp working at 4 amps. and cooled in running water. 
Although the sharpness of lines produced by this source does not approach 
that of a hollow cathode cooled with liquid air, it suffices, nevertheless, to provide 
all necessary information about adjustment and resolving power of the optics. 
Comparison with the known hyperfine structure of the Hg-spectrum (36) shows, 
even on the reproduction of Plate 1(6), that in the innermost interference fringes 
the two Hg-components A and (C +4) (36) with a distance of 0-04 a. can be 
resolved very easily and completely and that even up to about 16 orders. This 
proves that a resolving power of more than half a million is really reached without 
difficulty. ‘The component C of the one order, and b of the other, practically 
coincide. ‘They can be distinguished, however, thanks to the relatively high 
dispersion of the auxiliary spectrograph, by the fine spike protruding from the 
left end of the line (C+4) on Plate 1(6). From this we can read off for the 
range at this wave-length the value of about 0°36 a. 
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If sunlight is fed into the instrument and if the spectrograph slit is narrow 
(<0-15 mm.) one sees in the plane of the plate the continuous spectrum, crossed 
by the well-known system of parallel bright and dark parabolic arcs (22, 27-35). 
Contrary to the case of the Fabry-Perot interferometer, the parabolas with the 
Lummer plate are concave towards the red end of the spectrum. A suitable 
width of the parabolas can be achieved by adjustment of the slit width. 

A further easily recognizable visual and photographic test of the adjustment 
and resolving power of the instrument is offered by the line pair, 5883-91 Wv of 
intensity 5 and (5883-84 Fe of intensity 4 (37, 38). Indeed, it is possible to use 
the Doppler effect due to the solar rotation and to the relative motion of the 
Earth and the Sun to produce a doublet of variable wave-length separation, which 
is very valuable for spectral tests. 

The rotation of the Sun causes a Doppler shift of about 2 km./sec. = + 0-037 A.; 
the orbital motion of the Earth one up to 0-51 km./sec. = + 0-008 a. (according 
to the season), and the Earth’s rotation one up to 0-25km./sec. = + 0-004 A. 
(according to the time of day). From the East limb to the centre of the Sun’s 
disk one can separate this pair of lines easily and with certainty, using one of the 
large grating spectrographs of the usual tower telescopes, giving a resolution of 
about 150,000, From the centre to the West limb the separation becomes 
increasingly difficult and provides a sensitive test of the resolving power for 
Fraunhofer lines of any spectroscopic arrangement. At the West limb the 
distance apart of the two lines is only about 0-04 a. and their profiles overlap 
considerably. Plate 1(c) shows these lines, which have also been used by Treanor 
in testing his interferometric arrangement. 

Polarimetry.—The crystalline Lummer—Gehrcke plate is cut with its optical 
axis on the surface and perpendicular to the long edge of the plate. The double 
refraction of the quartz then produces two interference systems due to the 
ordinary and the extraordinary ray, which are linearly polarized perpendicular 
to each other. Since the refractive index differs for the two rays, these systems 
are separated by about 4 of the range. ‘The Fresnel reflection coefficient of the 
extraordinary ray is appreciably larger than that for the ordinary. As a result, 
during practical observations under comparable conditions, the number of 
effective reflections (i.e. of those corresponding to equal intensities) can be about 
twice as large for the extraordinary ray as for the ordinary one (39). ‘This 
number determines in the well-known manner the resolving power. Therefore 
the sharpness and the resolving power of the two interference systems differ 
appreciably. This can be noticed at once by visual observations of emission 
spectra as well as of channelled continuous spectra. Under these circumstances 
the original intention to use simultaneously the two interference systems for 





LEGEND TO PLATE 1. 


(a) Normal solar spectrum: one part artificially right-handed circularly polarized, the other part 
left-handed. Shows position of absorption bands caused by the thick quartz plate and polaroid. 
(Positive) 
(6) Hg 5461 hyperfine structure: McClean spectrograph crossed by quartz Lummer plate. 
Range 0°36 A., the distance from A to (C+5) is 0°04 A. (Negative) 
(c) Solar spectrum: McClean spectrograph crossed by Lummer-plate. Region of D-lines. 
(Positive) 
(d) As (c): Region of investigated Zeeman lines, range 0°33 A. (Positive) 
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relative measurements of the Zeeman effect of the same line on the same plate 
was abandoned (22). Instead the measurements were carried out with a device 
described elsewhere (40) based on Wood’s method of a thick quartz plate (41, 42). 

If we restrict ourselves to the observation of lines which are simply split 
according to the classical Lorentz triplet, we have to expect a practically pure 
longitudinal inverse Zeeman effect (43) at the two previously mentioned points 
on the Sun’s meridian. If there exists a magnetic field, such lines would consist 
of two circularly polarized components in asymmetrical position about the original 
line. We now consider two such lines A, and A, not too far apart in the spectrum 
(Fig. 1). It is then possible to cut a plane-parallel quartz plate parallel to the 
axis down to a thickness calculated to make it act for A, as a }A plate and at the 
same time for A, as a }A plate of high order. If we call m, and n, the refractive 
indices of quartz for the ordinary and for the extraordinary ray, respectively, 
and if the subscripts 1 and 2 refer to the wave-lengths A, and A,, then the thick- 
ness D of the quartz plate must fulfil the two conditions: 


D(n,- Ni) =A,(R+ }), 
D( Nog — Mem) =Aq(k +1 + 3). 


de A 


| | 
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Fic. 1.—Investigation of Zeeman effect with thick quartz plate and polaroid (schematic). 


Here k =0, 1, 2, 3,... and for A,>A, alsol =o, 1, 2, 3,.... If we permit a value 
of D of the order of a few millimetres, it is not difficult to find a sufficiently close 
approximation for D which can be adopted for the practical manufacture of the 
plate. ‘The required values of n,, and n, and other relevant optical data concerning 
this arrangement are compiled in (40). 

Amongst the numerous lines in the solar spectrum giving simple splitting 
and having large Landé factors there are unfortunately very few really close 
pairs suitable for this procedure (23, 40), and actually it proved impossible to 
find a sufficiently close pair with completely identical Landé factors. 

After a series of experiments with different combinations of lines the pair 
A5250:22 Fe, intensity 2, and A5263-32 Fe, intensity 4 of the Zeeman types 
0(3)/1 and 0(5)/2, respectively, was selected. These lines give interference fringes 
of especially fine quality on the fine grain plate used. An excellent crystalline 
quartz plate of the thickness D = 10-222 had been made according to the author’s 





8 H. von Kliber, An attempt to detect a general magnetic Vol. 111 


specification already in 1944 for this line pair, and it was used for the present 
investigation. The small corrections during the final figuring of the plate, 
necessary to get the exact phase-wave-length relation, were checked by a large 
grating spectrograph. This plate was placed, just as was the interferometer, 
within the parallel light beam before the slit of the auxiliary spectrograph. It 
converted the circular polarization of the Zeeman components of the two lines 
into linear polarization. But now the two outer components are linearly polarized 
in the same plane, and two inner ones in the perpendicular plane. This is 
illustrated in Fig. 1, which further shows a polaroid (Pol) inserted between the 
quartz plate and the spectrograph slit. According to the orientation of this 
polaroid only the two outer or only the two inner components are transmitted. 
We denote by d the normal difference in wave-length of the two lines in the 
absence of any Zeeman effect. If, however, a longitudinal Zeeman effect A, 
is present, we observe a wave-length difference (d+A, +A, or in the other 
position of the polaroid the difference (d—A,—A,,). Finally the difference of 
these two differences gives (if e.g. A,,=A,,) four times the value of the pure 
Zeeman effect, and this quite unfalsified by any Doppler effects. 

The complete spectroscopic-polarimetric arrangement is schematically 
represented in Fig. 2. The Sun is imaged by the two coelostat mirrors C and A 
and the lens L, on the screen S. The latter serves at the same time for focusing 
and for finding the right orientation of the Sun’s image. ‘Through a small 


L, #4 Lu USL = 


SH Q@ P G 


Ls pL 


Fic. 2.—Sketch of optical arrangement. 


opening on this screen, followed by a heat filter H, only the light from one of the 
above-mentioned points on the Sun’s central meridian enters the lens L,, at 
whose focus the screen S is placed. From there the light passes the thick quartz 
plate Q, a rotatable $A plate, and a polaroid P, and reaches the Gehrcke 
entrance-prism of the Lummer plate Lu. The two interference systems of Lu, 
one on each side of the optical axis, are imaged by the lens L, on the slit of the 
grating spectrograph. A field lens L, concentrates those rays which are inclined 
to the optical axis on the camera objective L,, bringing them also on the grating G. 
The grating spectrum, crossed by the two interference systems, is finally imaged 
on the photographic plate. This arrangement provides sufficient protection 
of all optical parts against direct heat radiation of the Sun, an important point, 
as the phase differences within the quartz plate Q are highly sensitive to 
temperature variations. The relevant difference (n,—1,,) of the refractive indices 
of the quartz at room temperature decreases by one unit in the sixth digit for 
every 1 deg. C. temperature increase. ‘This causes a displacement of the 
“‘ phase-positions "’ in the spectrum by about 0-7 a. towards the violet (40). 
It was desirable to keep the temperature of the quartz plate constant within 
1 deg. C. during the observations and this was achieved without difficulty. The 
Lummer plate is still more sensitive to temperature variations, but jacketed in 
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its asbestos box, of several centimetres thickness, the variations could be reduced 
to a small fraction of a degree. 

The rotatable $A plate behind the quartz plate serves when required to turn 
the direction of the linearly polarized Zeeman components by 90°. This makes 
it possible to keep the polaroid P which follows the plate in the same unchanged 
position. In spite of this, according to the position chosen for the 4A plate, 
either the two outer Zeeman components of the line pair only are transmitted 
or the two inner only. The result is that from here onwards only light of 
unvarying polarization passes through the whole apparatus (including Lummer 
plate and grating spectrograph), so that errors due to a changing state of 
polarization are completely avoided. The polaroid P is placed in such a way 
that only the extraordinary component of the Lummer plate is transmitted. 

Observations.—Before and after each exposure the proper adjustment of 
the whole apparatus was checked by visual observation of the sharp pattern of 
parabolas at different points of the interference system, using large magnification. 
This furnishes an experienced observer with an extremely delicate criterion of 
the adjustment. Furthermore, a small test object was placed within the light 
path near S, introducing artificially circularly polarized light of known sense 
into the incoming sunlight. It consisted of a polaroid followed by a }A plate 
for the proper wave-length, held in the right orientation. The neat device of 
A. E. H. Tutton (44) for the determination and checking of the right sense of 
the circular polarization should be mentioned in this connection. If one places 
this test object near to S in the light path, then the whole solar spectrum in the 
neighbourhood of the two lines will be circularly polarized. The presence of 
the quartz plate Q and of the polaroid P then produces in the spectrum at PL a 
series of dark absorption bands. Before every exposure it was confirmed visually 
that, for a given adjustment of the optical paths, one of these dark bands really 
coincided with one of the two lines under investigation. If the temperature 
necessitated it, a small correction was easily achieved by tilting the quartz plate 
through a small angle. This increased the optical thickness of the quartz plate 
and thus shifted the bands a little. If then the }A plate was turned by 45°, the 
bands were displaced in such a way that the other of the two Fraunhofer lines 
now coincided with a dark band. This important test was repeated at the end 
of each exposure in order to detect at once shifts due to temperature variations. 
In fact, even for exposure times of 30 min. or more, the position of the bands 
remained sufficiently constant. Exposures taken with the same emulsion as 
used for all exposures showed that the absorption was complete for all practical 
purposes. Plate I(a@) shows such a photograph, obtained for illustration 
purposes, after removing the Lummer plate. 

The observations started by moving the image of one of the two points on 
the Sun’s central meridian to the diaphragm of the screen S. It is well known 
that with such a coelostat arrangement it is impossible to rotate the Sun’s image 
in its plane by any required amount and to make its central meridian parallel to 
the spectrograph slit. A rotation of the image by means of additional reflections 
should be avoided 1n view of the polarimetric optics. ‘To a restricted degree, 
however, rotation is possible with such a horizontal mounting by choosing the 
position of the coelostat mirror C relative to the auxiliary mirror A (45, 46, 47). 
Let us denote by p (without regard to sign) the angle between the vertical 
drawn in the plane of the spectrograph slit and the hour circle which halves the 
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Sun’s disk. Further, let p be the angle between the Sun’s axis and this hour 
circle. We now want to give p such a value that the Sun’s axis becomes parallel 
to the slit. By spherical trigonometry it is easily shown that for an auxiliary 
angle 6 
° 
» wet . 
Cts = 90° 5-¢ veer 
cos —————— 
2 

Here 5 denotes the declination of the Sun and ¢ the geographical iatitude of 
the observer. We now call a the azimuth of the mirror A as seen from the 
mirror C and we obtain 
_ cos ¢ tan 6—sin ¢ cos 6 
4 sin 0 
The proper choice of a, within the geometrical limits of the whole mirror system, 
enables us thus to produce a suitable rotation of the Sun’s image. In so far as 
this was impossible, a small inclination of the Sun’s axis to the direction of the 
slit was unavoidable. ‘This is not a serious disadvantage, however, as only a 
strip of the Sun’s disk of about 31mm. (35 seconds of arc) in height is 
observable through the entrance of the Lummer plate. 

Next, the adjustment in the image plane of the spectrograph was performed 
by testing with a strong magnifying glass, as accurately as possible, the 
interferences and the proper positions of the absorption bands of the test object, 
caused by the quartz plate as mentioned above. Immediately afterwards the 
first exposure of one of the two points on the Sun’s central meridian at + 45° 
was made. ‘Then the photographic plate was changed, the } A plate turned by 45° 
and, without further touching the instrument, a second quite identical exposure 
was made. After the exposure of every plate the positions of the above-mentioned 
bands were tested again. In the subsequent measurement and discussion, 
plates were always paired together in this manner. After these exposures 
the visual test was repeated, if necessary taking a Hg-plate, keeping the 
optical adjustment completely unchanged. Experience gained during earlier 
observational series had led to the decision that in the present series the slit 
of the spectrograph was not to be opened so far as to bring about a contact of the 
various interference parabolas, creating a quasi-normal continuous spectrum 
on the plate (22, 32, 33, 34, 35). It is, indeed, evident that the contrast on the 
plates is remarkably increased if the parabolas are separated by bright intervals. 
This contrast is essential for attaining the highest accuracy in setting and 
measurement. Furthermore, a narrower slit also reduces the general scattered 
light within the spectrograph. The width of the slit finally adopted was about 
o-Imm., which corresponds to a width of the parabolas of about 0-2 a. on the 
plate. 

Backed Ilford thin film half-tone plates were used. For this contrasty 
but rather insensitive emulsion, the exposure times for the Sun were 
15 to 20 min., according to seeing. 

Care was taken to distribute the exposures as uniformly as possible for the two 
points of the Sun’s central meridian. Particularly favourable weather conditions 
during the Summer 1949 helped in obtaining during the months August and 
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September a very uniform and rather extensive series of some 60 plates of this 
kind, which forms the basis of this discussion. 

Measurements and reduction.—For the measurement of the plates a Zeiss 
spectrocomparator of the Cambridge Observatories was used. In this instrument 
the plate and a fixed measuring scale are mounted on a common movable table. 
The setting and reading is done with fixed microscopes in units of 0-000I mm., 
although the last digit is not really significant. ‘The measurements were 
arranged to be strictly differential and the same parts of the scale and drum 
were always used to avoid possible scale errors. The temperature variations 
during the measurements of one set hardly ever exceeded +1 deg. C. Six settings 
were made on each interference knot in symmetrical positions. When using a 
Lummer plate two symmetrically placed knots of the interference system play 
the same role as, the diameters of the interference circles in the case of a 
Fabry-Perot interferometer; together they furnish one particular wave-length 
value. In general, for every line and on every plate four or sometimes five pairs 
of knots were measured. ‘The best ones, which are those near to the plane of 
symmetry of the interference system, were selected. In order to derive one 
Zeeman-displacement according to Wood’s scheme, the difference for the lines 
A5250°22 and A5263-32 have to be formed on each plate of each pair, and then 
one of these two differences has to be subtracted from the other. This means 
that each value of a Zeeman displacement originates from the measurement of 
eight interference knots. If, by some small plate defect, etc., the measurement of 
even one of these knots was found to be uncertain, this whole Zeeman value 
was omitted. This was necessary in a few cases only and is the reason why on 
Fig. 3 sometimes only three values appear for a pair of plates instead of the usual 
four. ‘The measurements gave millimetre-values, which were converted into 
angstrom units in the usual way by using the value of the range of the 
Lummer plate. The determination of the latter at the position of the knots 
is conveniently done by counting the interference arcs between some known 
wave-lengths of the grating spectrum (23, 28). ‘The values of angstrom units 
finally give (apart from unavoidable random errors) a direct measure of the 
derived Zeeman-shift and thus of the magnetic field strength. We find from the 
diagram in Fig. 1, applying the Landé factor 3 and # for our two lines, that the 
angstrom unit values for a given field strength are nine times as large as those 
which would follow for the field strength from the classical Lorentz-formula 
with g=1. A value of 0-001 a. here corresponds therefore to a magnetic 
field strength of about 8-5 gauss. 

It may be advisable to give here some data relative to the internal accuracy, 
as they may also provide useful guidance in future similar measurements. The 
standard error (m.s.e.) of one setting of one interference knot was found to 
be between +0-002 and +0-005mm., according to the quality of the knots. 
This is in agreement with former measurements (10, 22). A field strength value 
derived from eight means, of six settings each, shows a m.s.e. +6 gauss. Each 
pair of plates furnishes, in general, four such values for the field. This brings 
the m.s.e. for each pair of plates to about + 3 gauss, and the m.s.e. for each of the 
two points on the Sun’s central meridian, each of which was derived from 
15 pairs of plates, to +0°8 gauss. The final result for these two points has this 
:m.8.e., as far as internal consistency alone is concerned. 
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The distances of any pair of interference knots (corresponding to the 
diameter of an interference ring in a Fabry-Perot system) are about 16-7, 15-3, 
13°8 and 1I2:Imm., respectively, counting from the outside inwards, the 
interference system being imaged with the lens L, of 25cm. focal length. 
The distances between two neighbouring knots, which decrease outwards, 
had a mean value of the order of 0-75mm. Since the range in the region of 
the observed lines is 0-33 A., this means that a decrease in wave-length by 0-2 A. 
increases the distance of two corresponding interference knots by I mm. 

The two lines examined are only 7 mm. apart on the plates. ‘The differences 
in the distances from the symmetry plane of the interference system are, for 
corresponding knots of both lines, on the average only 0-05mm. It can easily 
be verified that the possible systematic influence of differential line curvature 
remains quite negligible (40, 48). ‘Two interference knots of the one line and 
the two corresponding knots of the other, refer therefore practically exactly 
to the same point on the Sun’s disk. The light of all interference knots 
together comes from an area on the Sun of about 3xImm., or 35x12 
seconds of arc. We are therefore justified in considering all these knots as 
representative of one and the same small part of the Sun’s disk. The influence 
of Doppler effects of all kinds, in so far as they affect simultaneously both 
these two neighbouring Fe-lines, is completely excluded by the differential 
nature of the method. This feature of the work should also considerably 
reduce, or possibly entirely eliminate, other sources of error besides Doppler 
effect. The largest likely source is a change of scale between two plates of a pair, 
but the differential method should make this negligible. 

Results.—We have seen that the remaining residuals, after forming the 
differences between the measurements of the knots, represent—in so far as they 
are not simple random errors of the observations—the actual Zeeman effect 
and therefore the field strength in gauss. Fig. 3 shows the latter for the two 
selected points on the Sun’s surface. The single values for all pairs of plates 
are plotted at equal abscissae-intervals, i.e. without regard to the actual time 
intervals between the pairs. In the upper part of the figure each point is the 
result of the eight combined measurements of the interference knots of one pair 
of plates. Each pair gave, as a rule, four such independent values. The 
lower part of the picture shows the corresponding means for the plates. 
The simplicity and uniformity of the material did not seem to justify a more 
complicated reduction by weighted means, or least squares. 

Formally, at least, it should be possible to recognize a field strength of 
I gauss in the mean values for each of the two hemispheres. The plotted 
means for the one hemisphere should deviate systematically from the zero line 
as compared with those of the other hemisphere. If, however, the mean of 
all 15 plates for the northern hemisphere is taken, we obtain a south polarity 
of 0-08 +08 gauss (m.s.e.). The southern hemisphere also gives a south 
polarity, namely o-18+0-8 gauss. ‘The sign and size of these small quantities 
are such that they are evidently without physical significance. From the figures 
we can therefore safely conclude that this method, applied to the measurement 
of about 1000 interference knots, has disproved the existence of a magnetic field 
of more than about one gauss at the two observed points. 

The above-mentioned scatter of the single values is considered to be 
exclusively due to the unavoidable errors of measurement. It hardly appears 
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(b) 
Fic. 3.—Results from measured plates. 
(a) for the point +45° on the Sun’s meridian. 
(b) the same for —45°. 
Upper part : values for each pair of fringes. 
Lower part : means. 
Abscissae : dates of observation and plate numbers. 
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justifiable to conclude from the scatter that there exist small irregular fields in 
the region of the Sun measured, although, of course, such a possibility (within 
the indicated field-strength limits) cannot be ruled out completely. Stronger 
fields, e.g. those originating in single granules which have sometimes been 
suggested, would be completely smoothed out in view of the large observational 
area as compared with the small size of a granule, and also by the long exposure 
time. 

It should also be remembered that the observing period of the present 
series, contrary to the majority of the previous observations (1, 4, 5, 6, 10, 11), 
did not coincide with a sunspot minimum. ‘The last minimum was in 1944°2 
and the last maximum in 1947°5._ During the present series the Ziirich relative 
numbers exceeded 150 on some days and altogether the solar activity was not 
inconsiderable on most of the observing days. We have not tried to correlate 
the scatter in our field strength values with the momentary spot conditions of 
the Sun, but such a possibility has been kept in mind, in case it should be needed 
in future work. For this purpose a map has been sketched giving the sunspot 
positions for every observing day; the magnetic field strengths of some of the 
larger spots were determined in the usual manner (20), using the McClean 
grating spectrograph. ‘The resolving power and the dispersion of this 
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Fic. 4.— Synoptic chart for the time of observations (1949). 

@ active centres on the Sun’s disk. 

E occasional flares. 

+ points on the Sun’s meridian at which the present observations were made. 
instrument are, however, not favourable for this purpose. On the other hand, 
the regularly published values obtained by the Mount Wilson observers (49) 
give for every observing period only one value per spot group. It thus appeared 
advisable to tabulate the actual photographically measured field strengths for 
the more important spot groups visible on our observing days, as they became 
available from unpublished observations at the tower telescope in Potsdam 
(Einstein ‘Tower). Fig. 4 shows in a synoptic chart (§0), the most important 
centres of solar activity for the period of the Cambridge observations. We 
remember of course that such synoptic charts do not give us sufficiently accurate 
information as to the real effective condition of the Sun on a particular day. 
To achieve this they have to be supplemented by tables characterizing the 
development trends, which are regularly published (§1), where we also find 
additional details concerning the kind of activity, etc. In Fig. 4 the signs (+ ) mark 
the two observed points on the Sun’s central meridian. We note, after additional 
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examination of the daily sketches, that there was no active centre nearer than 
30° to one of the (+)-points. The following table compiles values for those 
groups for which effective field strength measurements are available at Potsdam, 
together with the classification of the group in the Ziirich nomenclature, the 
Ziirich relative number of the particular day, the Mount Wilson number and 
the published maximum field strengths of the Mount Wilson observers (49) 
(which, as mentioned above, do not necessarily refer to the same day). A number 
of unimportant spots with partly uncertain measurements of the field are omitted. 
It may, however, perhaps be of interest to study in future all available 
measurements of the general magnetic field of the Sun with a view to a possible 
correlation with the simultaneous solar activity. 


TABLE I 
Measured magnetic field strengths of sunspots visible during period of present observations 
5.8.49 109 N 10°-W 4° Dy Nio, Sis 9936 15 
8.8.49 5° N28 -E 43 J, Se: 9940 34 
12.8.49 56 19 -W11 J, Sec 9940 34 
13.8.49 82 19 -We2s_ J, Ses 9940 34 
-E 18 'C,, Sie, Se 9946 21 
W360 Jk Sas 9940 34 
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14.8.49 168 


15.8.49 


5 6 Cy Nu 995° 
5 32 Nu 
-E 37 pC S 9948 
5 42 
2 48 OE, 9949 
25 E;; I 9949 
20 26 i 9953 
E 23 ' S, 9955 38 
, , ) 9958 27 
9964 20 
9976 28 
9976 28 
9980 28 


2.9.49 
4.9.49 


9977 21 
997 35 
= 9982 16 
19 -E 45 a Sie, Sis 9983 18 
23 -E 53 Nis 9984 2 
14 -E 61 Nie 9987 32 


N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
S 
S 

S 

S 

Ss 

S 
N 
5S 

S 

N 
N 
Ss 


If, in conclusion, we recall the formula (18) already given for the field 
distribution H, on the surface of a homogeneously magnetized sphere, then 
the present series shows that at the time of the observations in 1949 August 
and September, within the formal m.s.e. of + 0-9 gauss, the polar field strength 
was too small to be measurable. Independent of the formal value of the m.s.e., 
Fig. 3 shows at once that a field of 1-2 gauss would have made itself evident 
beyond any doubt. 
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It is well known that the observations of the general magnetic field of the 
Sun, as far as they are available now, are ambiguous. The most recent series 
obtained at the Mount Wilson Observatory in 1940, etc. (§§) differ amongst 
themselves appreciably. G. Thiessen’s series of 1945 (11) gives a polar field 
strength of 53+12 gauss (m.s.e.), his later measurements (13, 14), however, 
all indicate a very small or unmeasurable field. In 1949 his observations with 
an interferometer and a photomultiplier give the value H, =1-5 with a m.s.e. 
of +0-7 gauss, which in view of the observational restrictions need not be taken 
as contradicting our present results. 

Recent researches into cosmic radiation (§2, §3) show that the existence of a 
general magnetic field exceeding a very few gauss is unlikely. Also a recent 
paper of P. A. Sweet (§4) makes the existence of a measurable magnetic field of 
the Sun seem very improbable on other grounds. 

On the other hand, the characteristic appearance of the solar corona and 
the periodic variability of its structure insistently suggest to us that a periodically 
variable general dipole field may exist, even if only a very weak one. ‘lo decide 
the question whether such a general and variable solar field has reality we 
need a still broader observational basis. Even a very weak dipole field of the 
order of less than one gauss, which would be in agreement with present-day 
observations, appears able to explain the various forms of the solar corona and 
other phenomena possibly connected with a general solar magnetic field. 


The writer wishes to express his gratitude to Dr A. Beer for many helpful 
discussions, and for his valuable assistance in preparing this paper. 


The Observatories, 
Cambridge: 
1950 November 6. 
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THE SUMMER DAYTIME METEOR STREAMS OF 1949 AND 1950 
I. MEASUREMENT OF THE RADIANT POSITIONS AND ACTIVITY 
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Summary 

Radio echo observations of the summer daytime meteor streams, first 
investigated in 1947 and 1948, have been continued with an apparatus 
specifically designed for continuous operation and routine radiant 
determination. Details of the radiant position and development of each 
daytime stream in 1950 have been recorded, together with information 
on the great sporadic activity which occurs during the same period. It is 
now possible to identify three permanent streams which have appeared 
regularly each year since their discovery. In addition to the permanent 
streams there are numerous streams which, although rich in meteors, are 
of a transient or long-period character. 





1. Introduction.—A high rate of transient radio echoes associated with meteors 
was found during the summer daytime observations made by Hey and Stewart 
in 1945 (3), and by Prentice, Lovell and Banwell in 1946 (4). ‘The significance 
of these results was not then apparent owing to the lack of knowledge of the 


relation between the number of radio echoes and the number of meteors. 
Subsequent work by Clegg, Hughes and Lovell in 1947 (1) showed, however, 
that these results represented an abnormally high level of meteor activity, 
associated with several well-defined radiants situated in the daylight sky. These 
results were confirmed and extended in 1948 (2). This work was carried out 
using a single directional aerial which was unsuitable for a continuous 24-hour 
survey. The present paper describes the continuation of the work on the 
daytime streams in 1949 and 1950 using an apparatus specifically designed for 
continuous operation and routine delineation of radiants. The results form part 
of amore general survey (§), but are presented separately here in juxtaposition with 
the velocity measurements (Part II) and the orbital calculations (Part II). 

2. Apparatus and technique.—In the original apparatus a single aerial was used, 
movable in azimuth, on which pulses of radio energy were received after reflection 
from ionized meteor trails (4, 6). Reflections occurred when a meteor stream 
was travelling at right angles to the aerial beam, and the radiant coordinates 
could be determined by observation of echoes along two azimuths using the 
method described by Clegg (7). This apparatus was not suitable for a general 
survey of meteor radiants, as the results obtained from one aerial are difficult 
to interpret, except in the case of single, well-defined streams. ‘lo minimize 
these difficulties, and to permit radiant coordinates to be determined in one 
operation, a new apparatus has been developed. ‘This apparatus, which has been 
described in detail elsewhere (8), consists essentially of two fixed, narrow-beamed 
aerials directed along azimuths 242° and 292° E. of N., operating on a wave- 
length of 4 metres. An active meteor shower produces echoes in each aerial in 
turn as its radiant moves across the sky. ‘These echoes are recorded by continuous 
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photography, and from their range-time variations it is possible to determine 
the right ascension and declination of a major shower to within +1}°. Before 
accepting a radiant position, a stream is required to satisfy certain conditions, 
which have been found necessary to separate the known night-time showers 
from random fluctuations in the sporadic background. On this system it is 
convenient to place a radiant in a class A, B or C weighted in order of importance 
as 2, 1, 0. The parameters of the apparatus have been adjusted so that the 
hourly rate of echoes corresponds closely to the visual hourly rate of meteors. 
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Fic. 1.—The day- and night-time meteor activity during the summer of 1950. 


Cuda Upper part: Maximum hourly rate of echoes during daytime 06 o0oM—18) oo™ U.T. 
Lower part: Maximum hourly rate of echoes during night-time 18% ooM-o6)} oo™ U.T. 


Abscissae: Longitude of Sun. 
3- Observation of the summer meteor activity of 1950.—Complete 24-hour 
records have been obtained for each day in 1950 May to August, which yield 
details of the development and extent of the daytime activity. Fig. 1 shows the 
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general activity during the summer, where the maximum number of meteors 
recorded in any day- or night-time hour is plotted. Certain increases in the 
echo rate correspond to the appearance of a meteor stream, and the radiant 
positions determined for individual streams are given in Table I. On days 
not included in the table no radiants were evident. In the case of a weak 
fluctuating stream, the weighted mean position of the radiant is given to indicate 
the centre of activity of the stream. For the intense Arietid stream an additional 
position is given for each day, representing the smoothed motion of the radiant. 
These radiant positions are shown on a Mercator projection of the sky in Fig. 2. 
The hourly echo rate given for each stream is a mean taken over the period of 
the range-time envelopes, excluding sporadic meteors. ‘This method produces 
values somewhat lower than those quoted in previous work, but gives a more 
reliable measure of the relative strengths of the streams. 


Degrees 
declination 


+20 


+10 —t- 
o-CETIDS oy TAQUARIDS 
0 od —$hep 
. | aguaRtus 
—10 5 vf 
CANIS 
—20 = 7 
110 100 90 80 70 60 50 40 30 20 10 OO 350 340 330 320 
Degrees Summer Radiants 1950 Mercator projection 
R.A, Epoch 1950 


Fic. 2.—The radiant positions of the 1950 daytime meteor streams shown on a Mercator 
projection. Individual radiant positions are plotted for the Arietids, n-Aquarids and o-Cetids. 
The diameter of a radiant represents the probable error in its position. For the ¢-Perseids and 
B-Taurids the weighted mean position is given and the area shows the scatter of the radiants in 
each group. 

4. The development of the summer activity 1950.—After the vernal equinox 
the sky became active in Aquarius and Pisces, and a single radiant was detected 
on May 6 near » Aquarii. This was immediately masked by the high background 
rate, and the activity spread to regions above and below the ecliptic near Pisces. 

o-Cetids. On May 14 a group of meteors was observed with a radiant near 
o Ceti. The hourly rate of 18 was comparable to the sporadic background, and 
it was not possible to detect a succession of radiants on following days. ‘Three 
similar radiants were found from May 16 to May 23, but the stream was weak 
and intermittent in character. After May 23 the radiant area again became large 
and followed the Sun along the ecliptic. 

Arietids. On June 2aremarkable progression of radiant positions commenced 
in Aries. For 12 days the radiant was observed to move with a definite but 
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abnormal motion away from the ecliptic, and was not deviated to any great 
extent by stream irregularities. The actual range-time plots covering this 
period, and from which the radiant positions were found, are given in Fig. 3. 
This diagram also illustrates the major character of the Arietids, which subse- 
quently proved to be the dominant stream of the sequence, reaching an hourly 
rate at maximum of 57. 

On May 31 and June 1, kefore the development of the main Arietids, two 
radiants were detected near 41 Arietis, but it was not certain that they indicated 
an independent stream. Pending further data these radiants have been placed 
with the Arietids in ‘Table I. 

(-Perseids. ‘The ¢-Perseids appeared during the Arietid period, but reached 
a maximum rate on June 4, four days before the Arietids. The radiant positions 
observed on each day for this stream are scattered over a region of sky 20° x 10°, 


TABLE I 


n-Aquarids 

Mean Position 
of Radiant 
R.A. Dec. 


May 6 12 338°5 3°0 << B 338'5 +-3°0 


Mean Hourly Radiant Position Radiant Class of 


Date Beko Base «2A. Dee. Diameter Raedient 


Daytime o-Cetids 


Mean Hourly Radiant Position Radiant Class of Weighted Mean 


> tase : te MRE: mens Position of Radiant 
Echo Rate R.A. Dec. Diameter Radiant aA Gen 


Date 


8 3°5 
8 3° +0°'o 
2 5‘0 

+1°O 


May 14 I 
May 16 I 
May 21 2 
May 2 I 


27°8 

31°0 

Daytime Arietids 

Smoothed Position 


of Radiant 
R.A. Dec. 


— Mean Hourly Radiant Position Radiant Class of 
Echo Rate R.A. Dec. Diameter Radiant 

May 31 18 43°0 +26°0 

June 20 43°0 27°0 

June 36 39°0 +16°0 

June 39 39°5 + 185 

June 41 42°5 +17°5 

June 31 41°5§ +20°0 

June 4! 42°0 +21°0 

June 30 44°5 + 23°5 

June 57 45°5 21°0 

June 31 46-0 25°0 

June 37 48°5 23°5 

June 23 46:0 +26°5 

June 23 46°5 +25°0 

June 32 48-0 27°0 

June 38 

June 22 No position determined : 

June 19 records masked by 

June 16 intense radio-frequency 

June 14 emission from the Sun. 


39°7, +171 
40°5 18°0 
41°3. +189 
41°9 +198 
42'8 +20°8 
43°5 +21°7 
44°3 +22°6 
45°I +23°5 
45°38 +24°4 
466 +25°3 
47°3 +262 
48'°0 +27°0 
48-7. +28-0]* 
49°5 +28°9 
50°3 +29°9 
50°9 _+30°8 
516 +31'8 
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* Extrapolated positions. 
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Fic. 3.—Diagram showing portions of the range-time plots covering the period of the Arietid and 
{-Perseid streams. 
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TABLE I (continued) 
Daytime ¢-Perseids 


Mean Hourly Radiant Position Radiant Class of Weighted Mean 
Echo Rate R.A. Dec. Diameter Radiant Position of Radiant 


A 
B 
A 
B 
B 
A 
A 
B 
B 
B 
B 
A 
B 
B 


Date 


° ° 
20 58:0 +28°0 
20 §7°5 +27°0 
39 53°5 + t9’5 
29 §8°5 +14°0 
18 60°5 
20 54°5 
June 22 62-0 
June 2 63°5 
June 9 66-0 
June 10 17 63° 
June 12 62°5 
June 14 69°5 
June 15 * 
June 16 71's 


June 
June 
June 
June 
June 
June 
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Daytime B-Taurids 


Mean Hourly Radiant Position Radiant Class of Weighted Mean 


Date Echo Rate R.A. Dec. Diameter Radiant Position of Radiant 


5 


June 26 14 80°% B 
June 27 23 87:5 
June 28 15 , 86° 
June 29 14 88: 
June 30 12 86: 
July 1 19 85° 
July 2 29 86:5 
July 3 21 86- 
July 4 20 88+5 
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although individual radiants were seldom more than 3° in diameter. ‘There was 
a general tendency for the radiant to move eastwards but this motion was too 
irregular for a smoothed position to be given as in the case of the Arietids. 

B-Taurids. After June 18, a diffuse radiant again appeared near the Sun, 
until on June 26 a centre appeared to the East of the main activity. Further 
observations indicated a radiant which moved at random from day to day in an 
area of sky 8° x 14° in Taurus. This stream reached an hourly rate of 29 on 
July 2 and appeared to be similar in character to the (-Perseids. 

From July to August the daytime activity gradually subsided and no further 
radiant positions were delineated. 

5. Observation of the summer meteor activity of 1949.—During the summer of 
1949 only one fixed aerial of the new apparatus for determining radiants had been 
completed, so that the streams could not be studied as effectively as in 1950. 
Subsidiary information was obtained during part of the summer, using the 
original movable aerial, but the data proved to be inadequate for determining 
radiant positions. The broad features of intense meteor activity, dominated by 
the Arietids, were, however, still apparent. 

6. Discussion of the observations from 1947 to 1950.—As a result of four years’ 
observation it is now becoming possible to distinguish the permanent streams 
which have recurred regularly at the same solar longitude each year, with 
approximately the same radiant positions. A comparison is made in Table II 
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of the radiants of each year, and four permanent streams are identified in the 
final column. ‘The 7-Aquarids, which have been previously established by 
visual methods, begin the sequence, and indicate the yearly variations to be 
expected in a stream. ‘The f-Taurid stream has variations greater than the 
n-Aquarids. ‘These differences are attributed in part to irregularities in the 
stream itself, and partly to errors in the early observations when declination 
measurements were vitiated by radiant movement between aerial settings. ‘The 
other permanent streams, the Arietids and ¢-Perseids, occur simultaneously, 
and were not at first successfully interpreted with the single aerial. With the 
new apparatus, however, these streams are clearly separated as shown in Fig. 3. 
Of these permanent streams it now appears certain that the Arietids produce the 
richest meteor shower, surpassing the Perseids and Geminids in duration and 
intensity. Fig. 4 illustrates the grandeur of this shower, where the hourly 
rates of Arietids and Geminids, as determined by the new apparatus, are plotted 
for comparison. 
May June 1950 June 
wot 2 Fe 4@ § SS 8 9 10 I! 12 13 14 15 16 17 18 





DAYTIME 
ARIETIDS 











7 8 9 10 Wb 12 13 14 15 
Dec, 1949 Dec. 
Fic. 4.—Comparison of the mean hourly rates of the night-time Geminid shower in 
1949 December, and the daytime Arietid shower in 1950 June. 
Ordinates: Mean hourly rates of echoes. 

In addition to the permanent streams there are certain streams which have 
not recurred regularly each year. Some of these—the 54-Perseids at © =93°, 
for example—have produced high hourly rates in the past, and may have a 
pronounced periodicity. Others, such as the system between 46°<(©<59° 
are weak and ill defined and may represent spurious groupings in the sporadic 
background. * 

* Note added in proof: The development of the activity in 1951 has been similar with the 
Arietids and {-Perseids appearing prominently at the end of May and early June. 
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The general eastward trend of both permanent and transient streams, and 
the apparent focusing of radiants near the solar point, is explicable in terms of 
elliptical orbits, and will be dealt with more fully in a later publication (5). 

The remaining component of the daytime display consists of the numerous 
meteors which cannot be assigned to definite radiants by the apparatus so far 
used. In each year this sporadic rate has correlated approximately with the rise 
and fall of the ecliptic, and it remains to be seen whether this indicates that the 
orbits of the meteors are of low inclination or whether there is a definite asymmetry 
in the distribution of meteoric matter round the Sun. 

7. Acknowledgments.—Both authors wish to acknowledge the advice and 
encouragement given by Dr A. C. B. Lovell throughout the work. ‘This work 
forms part of the general research programme at the Jodrell Bank Experimental 
Station of the University of Manchester, to which the Department of Scientific 
and Industrial Research have contributed substantial financial aid. Both 
authors are also indebted to the Department individually for the award of 
maintenance grants. 
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THE SUMMER DAYTIME METEOR STREAMS OF 1949 AND 1950 
II. MEASUREMENT OF THE VELOCITIES 


FJ. G. Davies and F. S. Greenhow 
(Communicated by A. C. B. Lovell) 


(Received 1950 December 22) 


Summary 


This paper describes the measurements of the velocities of the summer 
daytime meteor streams using the radio echo diffraction technique. The 
observations were made simultaneously with the radiant determinations 
described in Part I of this series of papers. Mean geocentric velocities 
of the o-Cetid, Arietid, ¢-Perseid and f-Taurid meteor streams have been 
determined, and it is shown that there is a considerable spread in the 
velocities of the individual meteors belonging to each shower. 





1. Introduction.—The radio echo measurement of meteor velocities by 
observation of the diffraction pattern as the ionized column forms in the atmo- 
sphere has been described previously (1, 2). During the summer months of 
1949 and 1950 this technique has been applied in an attempt to determine the 
velocities of the daytime meteor streams (3, 4,5). ‘The results of this work, 
which was carried out simultaneously with the measurements of the radiant 
coordinates described in Part I, are described in the present paper. 

2. Technique.—A steerable aerial array (6) with a beam width of +7’, working 
on a wave-length of 4-2 metres, was used throughout the experiments in conjunc- 
tion with electronic recording equipment similar to that previously described (2). 
The use of this directional array facilitated the measurement of the velocities 
of meteors from specific radiants. For this purpose two techniques were 
employed. 

(a) In the case of multiple or unknown radiants it was necessary to maintain 
the aerial at a fixed azimuth and low angle of elevation in order to obtain range- 
time plots similar to those described in Part I. By using the theoretical 
range-time distribution as a criterion it was possible to select the measurements 
appropriate to the shower under investigation. This procedure led to the 
routine rejection of non-shower meteors lying outside the theoretical range- 
time envelopes. In this method there is no criterion whereby non-shower 
meteors inside the envelopes can be identified. 

(6) In the case of a very active radiant of known position the azimuth of the 
aerial was adjusted over a number of hours so as to maintain the maximum echo 
rate. 

From 1950 May to July, systematic velocity measurements were made in 
this manner. During 1949 it was not possible to carry through a regular 
programme, and the results are restricted to the period May 30 to June 18. 

3. Results—The complete list of velocities obtained in 1950, together with 
other relevant information, is given in Table I; the 1949 velocities are given 
in Table II.* The results, related to the principal radiants described in Part I 
are discussed in detail below. 


* For Tables I and II, see pp. 31, 36. 
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(a) o-Cetids. According to the data given in Part I, the Cetid stream was 
the first clearly defined radiant to be observed in the summer of 1950. ‘This 
shower was active from May 13 to May 23, and the velocities measured during 
this period are given in Table I. 

The results described in Part I show that a single isolated radiant was present 
with a rate comparable to that of the background activity. Measurements were 
made at three different aerial azimuths, chosen in such a way that the radiant 
was able to pass completely through the beam in each position. ‘The range- 
time envelopes for the Cetid radiant were computed for each aerial direction, 
so that the appropriate selection could be made for each position. 

The velocity distribution, showing a single group between 28-9 and 
44:3 km. )sec., is given in Fig. 1. The mean velocity is 36-8km./sec. with a 
standard deviation of 4-2 km./sec. 


104 


8« 














[ino 


30 40 
Fic. 1.— Distribution of o-Cetid velocities. 


Ordinates: Number of echoes. 
Abscissae: Velocity (km./sec.). 














(b) Arietids and ¢-Perseids. ‘The Arietid and ¢-Perseid meteor streams were 
both active from approximately May 31 to June 18. Since the two radiants 
were separated by only 20° in right ascension, it was necessary to make 
the measurements with the aerial at a fixed azimuth in order to resolve the 
individual radiants. This technique was applied between May 31 and June IT. 
Fig. 2 is a range-time plot of the 78 echoes for which velocities were determined 
during this period. Theoretical range-time envelopes for the Arietid and 
¢-Perseid radiants are also shown in the figure. ‘The large concentration of echoes 
inside the Arietid range-time envelope is very marked; there is also a smaller 
disiribution corresponding to the (-Perseid radiant. ‘The velocity deter- 
minations of the echoes lying in or near each theoretical range-time envelope 
are plotted separately in Figs. 3a and 6 and superimposed in Fig. 3c. The 
division into two distinct groups, corresponding to the Arietid and ¢-Perseid 
radiants is evident. The 44 Arietid velocities have a mean of 37-8km./sec., 
with a standard deviation of 4-0km./sec., and the mean of the 21 ¢-Perseid 
velocities is 28-8km./sec., the standard deviation being 3:2km./sec. Thirteen 
velocities lying outside the theoretical envelopes were rejected as non-shower 
meteors. In addition one velocity with a value of 53-7 km./sec. lying inside the 
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f-Perseid envelope was rejected, since it was outside the probable limits of the 
velocity distribution. 

After June 11 the aerial was moved continuously in azimuth, in order to 
increase the number of results for the Arietid shower. A further 78 velocities, 
with a mean of 37-6km./sec. were obtained over this period. ‘The distribution, 
shown in Fig. 4, has a single peak at approximately 37-5 km./sec. 
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Fic. 2.—Range-time distribution of echoes giving velocity determinations (May 31—JFune 12). 
Ordinates: Range in km. 
Abscissae: Hours in U.T. 
x : Arietid echoes. 
: ¢-Perseid echoes. 
: Theoretical range-time distributions. 
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Fic. 3a.—Distribution of Arietid velocities. Fic. 36.—Distribution of €-Perseid velocities. 
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The mean of all the Arietid velocities obtained between May 31 and June 18 
is 37°7 km./sec., with a standard deviation of 4-3km./sec. There is no significant 
difference between the mean velocities of the stream determined over the periods 
May 31 to June 11 and June 12 to June 18. 
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Fic. 3c.—Distribution of Arietid and ¢-Perseid velocities. 
Ordinates : Number of echoes. 
Abscissae : Velocity (km./sec.). 
Dates of observations : May 31—Fune 12. 






































i 


Fic. 4.—Distribution of Arietid velocities (Fune 12-—Fune 16). 
Ordinates : Number of echoes. 
Abscissae : Velocity (km./sec.). 
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In 1949 there was also a large increase in the meteor echo rate during this 
epoch. Although no radiant positions were determined in that year, on applying 
the theoretical range-time envelopes corresponding to the Arietid radiant 
delineated in 1950, a group of 22 velocities with a mean value of 38-5 km. sec. 
was obtained (‘Table II). 

(c) B-Taurids. The 8-Taurid radiant was active over the period June 27 
to July 4, and velocity measurements were made between these dates. Owing 
to the comparatively low hourly rate, the measurements were made at a fixed 
aerial azimuth in order to apply the range-time criterion for selection of the 
B-Taurid meteors. ‘Twelve velocity determinations from meteors satisfying 
the range-time conditions for the 8-Taurid radiant, were obtained (Fig. 5). 
The group of 10 velocities between 23-7 and 36-7 km./sec. has a mean value of 
31'4km./sec. and a standard deviation of 4:1km./sec. The two other velocities, 
19°8 and 44:3km./sec., have been rejected as probable sporadic meteors. 


— 














, LJ 
20 30 40 
Fic. 5.—Distribution of B-Taurid velocities. 
Ordinates : Number of echoes. 
Abscissae : Velocity (km. /sec.). 





4. Discussion of results.—The results of the measurements of the velocities 
of the daytime meteor streams, made during 1949 and 1950, are shown in Table III. 

The velocity distribution for each shower shows an appreciable spread about 
the mean. The standard deviation for all cases is approximately 4km. sec. 
This is partly due to the presence of sporadic meteor velocities in the groups, 
but even in the case of the Arietid radiant, where the shower rate was greatly in 
excess of the background activity, the spread in velocities remains large. ‘The 
experimental error in each individual velocity determination is given in Table I. 
The r.m.s. value of the error in the Arietid group is 2-4km./sec., and the spread 
in the velocity distribution due to errors introduced experimentally is therefore 
of this order. As the standard deviation of the group is 4-3 km./sec., it is clear 
that there is a true spread in the velocities of individual meteors of approximately 
3°6km./sec. ‘The significance of this spread is discussed in Part IIJ of this paper. 

5. Acknowledgments.—The authors are indebted to Miss M. Almond for her 
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Director of the Jodrell Bank Experimental Station, where this work was carried 
out, for his help and encouragement, and his advice in the presentation of this 
paper; and to A. Aspinall and G. S. Hawkins for their ready cooperation. 
One of the authors (J. S. G.) wishes to thank the Department of Scientific and 
Industrial Research for financial assistance. ‘The work forms part of the general 
meteor programme at Jodrell Bank which receives substantial financial aid from 
the Department of Scientific and Industrial Research. 
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TABLE I 
Date 1: Range Aerial Velocity 
lime ' ' 
1950 (km.)  azimuth(°) (km./sec.) 
May 1 0749 2 ogo 35°14 
May 6 0609 ogo 30°54 
0728 ogo 66-7 
May 0846 040 36° 
0917 050 
0957 3: 070 
May Inknown : ogo 


Shower 


Inknown ogo 

May 11 0945 ‘ o60 
May 1005 ‘ ogo 
1107 § ogo 

1129 : ogo 

1139 55 ogo 

1141 ogo 

1152 ogo 

Og14 § 060 

Ogi5 o60 

0926 o6o0 

0927 c 060 

0929 5 c6o 

0931 o60 

og50 o60 

1053 ogo 

1101 ogo 

1103 : ogo 

1106 ogo 

1133 ogo 

1153 ogo } 

0857 060 *3+3° o-Cetid 

ogo8 o60 o-Cetid 

0923 o60 o-Cetid 

1018 ogo o-Cetid 

1030 ogo o-Cetid 

1033 ogo o-Cetid 

1057 ogo o-Cetid 

1058 ogo o-Cetid 

1121 2 120 o-Cetid 

o-Cetid 
o-Cetid 


o-Cetid 
o-Cetid 
o-Cetid 
o-Cetid 
o-Cetid 
o-Cetid 
o-Cetid 
o-Cetid 
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o-Cetid 
o-Cetid 
o-Cetid 
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0932 o6o0 
10271 ogo 
1128 ogo 
0935 ogo 
0940 ogo 
1004 ogo 
1032 ogo 
1039 ogo 
1047 ogo 
1102 ogo 
1025 ogo 
1032 3 ogo 
1112 ogo 
1122 ogo 
1127 ogo 
0950 ogo 
1018 ogo 
1156 090 
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0926 280 ogo 34°31 
1044 550 ogo 360°4+2 
1045 370 ogo 349+ 
1112 430 ogo 30°4-+4 
1028 130 ogo 30°3-+ 
1057 440 ogo 28-3 
1124 280 ogo 36°14 
1125 300 ogo 30°1-+4 
1138 290 ogo 25°32 
1143 240 ogo 28°8-+ 
1145 390 ogo 40°1+ 
1157 240 ogo 19°2+ 
1008 445 090 47°04 
1044 ogo 34°84 
1117 230 ogo 26°5-+ 
1118 325 ogo 38°34 
1129 ogo 2°4-++ 
1132 330 ogo 38°04 
0920 ogo 48°3-++ 
1024 ogo 
1045 ogo 
1047 ogo 
1132 ogo 
1134 ogo 
1153 090 
1226 090 
1232 ogo 
1012 4 ogo 
1040 ogo 
T1I5 ogo 
1136 35 090 
1210 ogo 
0938 ogo 
1038 ogo 
1039 35 ogo 
1042 ogo 
1132 090 
0946 ogo 
0956 090 
1057 ogo 
1112 ogo 
1145 ogo 
1020 ogo 
III ogo 
0939 ogo 
1032 ogo 
1053 ogo 
1055 ogo 
1059 ogo 
1133 ogo 
1142 ogo 
1148 ogo ¢-Perseid 
1155 ogo C-Perseid 
1156 ogo C-Perseid 
1033 ogo Arietid 

1047 090 Arietid 

1109 ogo Arietid 

1109 ogo Arietid 

1112 ogo Arietid 


o-Cetid 
o-Cetid 
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1022 410 ogo 
1046 350 090 
1116 385 120 
1120 580 120 
1125 410 120 
1135 410 120 
1145 350 120 

June 20 0956 330 ogo 
0959 370 ogo 
1002 430 ogo 
1010 440 ogo 
1023 330 ogo 
1041 360 ogo 
1105 370 ogo 
1206 405 ogo 
1219 340 ogo 
1220 300 ogo 
1227 480 ogo 

June 28 1042 350 ogo 
1045 400 ogo 
1051 330 ogo 
1056 310 ogo 
1110 430 ogo 
1131 390 ogo 
1156 240 ogo 
1204 370 ogo 
1220 430 ogo 

June 2 III! 390 ogo 
1120 415 ogo 
1123 415 ogo 
1152 550 ogo 
1210 39° ogo 
1243 430 ogo 
1246 360 ogo 
1031 290 ogo 
1049 37° ogo 
1125 360 ogo 
1130 410 ogo 
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0938 410 033 
095! 54° 045 
IOIr 345 045 
1042 100 058 
1043 340 058 
1121 290 ogo 
1201 350 105 
1007 370 ogo 
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1029 420 ogo 
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1051 37° ogo 
1102 390 ogo 
1106 365 ogo 
1106 240 ogo 
1120 330 ogo 
1126 390 ogo 
1133 300 ogo 
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1209 350 090 
1217 240 ogo 
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Tasie II 


Range Velocity 
(km.) (km./sec.) 
255 ogo 40°5+1°3 Arietid 
305 ogo 32°7+1'5 Arietid 
385 ogo 40°9-+2°8 Arietid 
235 ogo 53°30°8 
295 ogo 37°3+0°7 Arietid 
300 ogo 31°2-+1°6 
193 ogo 45'O+2°2 
333 120 37°I+1°3 Arietid 
358 120 39°7+2°0 Arietid 
320 120 33°8+2°7 Arietid 
385 ogo 43°1I+2°1 Arietid 
420 ogo 42°0+2°4 Arietid 
312 ogo 44°7+4'5 Arietid 
382 ogo 34°9+2'5 Arietid 
240 ogo 29°5+2°2 Arietid 
285 ogo 45°1+2'1 Arietid 
218 ogo 41°7+3'8 Arietid 
300 ogo 40°5+1°1 Arietid 
223 ogo 41°0+3'5 Arietid 
225 ogo 34'9+4'5 Arietid 
300 ogo 31°8+1°2 
472 ogo 26:0-+0°6 
184 ogo 36°8+2°3 
357 ogo 36°2+0'9 
335 060 65°3+4'8 
277 o60 41°8+2°4 Arietid 
225 ogo 38:7+2°3 Arietid 
240 ogo 37°8-+1°3 Arietid 
195 ogo 35°7+2°4 Arietid 
240 ogo 34°0+3°8 Arietid 


Azimuth (°) Shower 


TABLE III 


Mean Standard Number of 
Shower velocity deviation determina- 
(km./sec.) (km./sec.) tions 
o-Cetid (1950) 36°7 4°2 37 
Arietid (1950) 

(i) May 31 to June 11 37°8 4°0 
(ii) June 12 to June 16 37°6 4°5 
(iii) May 31 to June 16 37°7 4°3 

t-Perseid (1950) 28°8 3°2 
f-Taurid (1950) 31°4 4°1 
Arietid (1949) 38°5 4°0 


~ 
nN -— N NSIT 
nv Or N Dh 


Jodrell Bank Experimental Station, 
University of Manchester : 
1950 December. 


References 


(1) C. D. Ellyett and J. G. Davies, Nature, 161, 596, 1948. 

(2) J. G. Davies and C. D. Ellyett, Phil. Mag., 40, 614, 1949. 

(3) J. A. Clegg, V. A. Hughes and A. C. B. Lovell, M.N., 107, 369, 1947. 
(4) A. Aspinall, J. A. Clegg and A. C. B. Lovell, M.N., 109, 352, 1949. 
(5) C. D. Ellyett, M.N., 109, 359, 1949. 

(6) A. C. B. Lovell, C. J. Banwell and J. A. Clegg, M.N., 107, 164, 1947. 





THE SUMMER DAYTIME METEOR STREAMS OF 1949 AND 1950 
III. COMPUTATION OF THE ORBITS 


Mary Almond 
(Communicated by A. C. B. Lovell) 
(Received 1950 December 8) 


Summary 


The orbital elements of four daytime meteor streams have been computed 
from the observational data given in Parts I and II. The orbits are found 
to be of short period and to lie within the orbit of Jupiter. The origin of 
these meteor streams is discussed. ‘Two are of low inclination and are 
shown to be associated with well-known visual streams produced by the 
passage of the Earth through the orbit before perihelion. The main 
stream of the daytime sequence is shown to have a probable association with 
the visual 5-Aquarid stream. One of the orbits appears to be coincident 
with that of Encke’s comet, but there are no prominent objects which can 
be associated with the others. 





1. Introduction.—The orbits of many of the major meteor showers recognized 
by visual observation have been inferred on the assumption of a parabolic velocity. 
In a few cases—for example the Bielids, Giacobinids and Leonids—the orbital 
elements have been computed from their periodic returns. For the Taurids and 
the Geminids exact orbits have been computed from actual measurements of the 
radiant positions and velocities (x, 2). Of these streams Whipple (1) has associated 
the orbit of the Taurids with that of Encke’s comet, but the period of the Geminid 
stream (1-7 years) is much shorter than that of any known comet. 

The computation of the orbits of the recently discovered major streams active 
in the summer daytime is therefore a matter of considerable interest. The results 
described in Parts I and II of this series of papers have now given reliable data for 
the radiant coordinates and velocities of four streams. In this paper the orbits of 
these streams are computed and the results are discussed in relation to their 
probable origin. 

2. Observational data and computation of orbital elements.—The 1950 observa- 
tions described in Parts I and II refer to the o-Cetid, Arietid, f-Perseid and 
f-Taurid meteor snowers. ‘The observational data for these showers which were 
obtained with an aerial system at latitude A= + 53° 14’ and longitude ¢ = 2° 18’ W. 
are given in Table I. 

TABLE | 
Radiant Positions and Velocities of Four Daytime Meteor Streams 
o-Cetids Arietids ¢-Perseids B-'Taurids 
Mean date 1950 May 19 June 8 June 8 June 30 
Radiant « (degrees)  29°4+2 44°3+1 61-642 86-242 
Position 0 mi - 2°7+2 + 22°641°5 +238+2 +18°7+42 
Velocity V (km./sec.) 36°7+.4°2 37°744°3 28°8+4.3°2 3144411 


The radiant positions were corrected for zenith attraction and diurnal aberra- 
tion, the necessary corrections being in no case more than one degree, and the 
resulting positions were used to compute the orbital elements given in Table IT. 
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Tasie II 
Fundamental Data and Orbital Elements of Four Daytime Meteor Streams 
o-Cetids Arietids ¢-Perseids B-Taurids 
Corrected « (degrees) 29°0 +2 43°8 +1 61°0 +2 85°7 +2 
Radiant 56 za 3°97 +2 21°9 +1°5 +228 +2 +179 +2 
A i 26-0 +2 480 +1 63°6 +2 85-9 +2 
2 


” 136 +2 5°7 +1°5 + 3°O 4°6 +2 
Heliocentric velocity 


Vq (km./sec.) 32°5 +2°7 34°0 +2°9 34°8 +2°0 ‘7 £28 
a (A.U.) 1*3 +04 1°5 +0'7 1°7 +0°5 =a +9 
—o'8 
e o'91+0°04 0'94+0°04 0°79-+0°'08 o°85+0°10 
q (A.U.) O*11-+0°03 0°10+0°'03 0°35+0°04 0°34+0°04 
i (degrees) 34 +7 18 +5 4 +2 6 43 
wo » 211 +3 29 +2 61 +4 244 +4 
7 oe 89 43 106 +2 138 +4 162 +4 
Q be 238°0 77°O 77°0 278°1 
AQ, 9 18 16 8 

The duration AQ of the showers varied from 8 to 18 days, and although the 
radiant position was found to shift from day to day, the movement would not be 
sufficient to produce a detectable change in the meteor velocities. Orbits were 
therefore computed from the mean radiant positions found over these periods, 
together with the corresponding most probable velocity values. 

The dates of observations, which determine the nodes of the orbits, were taken 
to be the mean days of the o-Cetid, {-Perseid and £-Taurid showers, although the 
date of maximum activity was used for the Arietids. 

3. Probable spread in the orbits.—The spread in the orbital elements given in 
Table II was estimated by considering the effects of a movement of the radiant 
within the given limits, and also of the probable spread in velocity. 

The values of the spread quoted for the elements a, e and g, which give the size 
and shape of the orbits, were determined by the error in velocity. The effect of a 
movement of the radiant was negligible in comparison. 

The inclination 7 was influenced by the variation in the velocity and in the 
declination 5 of the radiant. The variation in the right ascension « of the 
radiant had the greatest effect on the longitude of the perihelion w. ‘The value of 
w would be altered significantly by taking a different position of the node Q, but 
this effect has not been considered in computing the mean orbits of the streams. 

m is the sum of w and Q, and so has the same limits as w. 

4. Discussion.—The mean orbits of the four daytime meteor streams projected 
on to the plane of the ecliptic are shown in Fig. 1 and are discussed individually 
below. 

(i) Arietids.—The radiant position of the Arietid stream shows a steady change 
in right ascension and declination throughout the shower, and two orbits, with the 
same semi-major axes as the mean orbit, have been computed using the smoothed 
radiant positions for the first and last dates of observation (Table III). Although 
the shape of the orbit remains practically unchanged, it can be seen that 7 and 7 
show considerable variation with position in the shower, such that the line of 
apsides swings gradually forward, and the main centre of the stream rises farther 
away from the ecliptic day by day. 

The portion of the stream producing a shower at the beginning of June has 
particularly small inclination, and so could be seen at a point of the orbit other than 
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the node, therefore a meteor shower should occur when the stream cuts the Earth’s 
path again after perihelion passage. From the computed orbit it can be deduced 
that the radiant position for this return shower would be « = 336°,5=—11° on 





C-PERSEIDS 


ARIETIDS 


JUPITER 


Fic. 1.—Orbits of summer daytime meteor streams projected on to the plane of the ecliptic. 
Tase III 
Limiting Orbits of the Arietids 
Ist day Mean day Last day 
Corrected « (degrees) 38°4 44°3 51°6 
Radiant 5 5 +15°5 +22°6 +31°8 
Date June 1 June 8 June 18 
a (A.U.) rs $s 15 
e 0°94 0°94 o'92 
q (A.U.) 0'09 o'10 o'12 
i (degrees) 3 18 34 
w ms 28 29 33 
7 "4 98 106 120 
Q Pr 70°3 77°0 86-6 
This radiant coincides with that of the visual 5-Aquarid stream (Table IV) 
which produces only a minor shower in the northern hemisphere, but forms one of 
the most prominent southern meteor showers having a duration of about 20 days. 
Various values have been given for the velocity of the 5-Aquarid meteors. 
McIntosh (7) estimated from the zenith attraction a hyperbolic value of 
46°51 km./sec., but Davidson (8) later showed that this method of calculating 
velocities was incorrect. Prentice (g) quoted a much lower observed value of 
29km./sec., that is a heliocentric velocity of 32 km./sec., while Hoffmeister (5) has 
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determined the heliocentric velocity by an indirect method, and obtained the 
value 34°6km./sec., which is very close to the Arietid velocity of 34-0km./sec. 
Hoffmeister computes an orbit similar to that of the Arietids, but having different 
w and an inclination of 23°+7 (Table V, Fig. 2). Such an orbit would not approach 
the Earth again in June. 


TABLE IV 
Radiant Positions of 5-Aquarids 
Return 5-Aquarids 
of ~- A ’ 
Arietids Norton (3) McIntosh (4) Hoffmeister (g§) Opik (6) 
a(degrees) 336° 339° 345° 344° 336° 
8 a —11 —t31 —17°0 —16 —II 
Date July 28 July 25-30 July 28 August 3 July 29-30 
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Fic. 2.—Orbits of Arietid and 8-Aquarid meteor streams projected on to the plane of the eckptic. 
TABLE V 
Orbits of Daytime Arietids and 5-Aquarids 
Daytime Arietids 5-Aquarids 

Date June 8 August 3 
Vay (km./sec.) 34°0 34°64 
a (A.U.) I°5 1°5995 
e 0°94 0°9264 
q (A.U.) o'18 o'118 
i (degrees) 18 23°7 
w - 29 147°0 
7 a 106 97°0 
Q = 77°0 310°0 
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However, the 5-Aquarid radiant is known to cover an appreciable area (10), and 
some of the meteors travel much closer to the ecliptic than those in the orbit given 
above. Although the main streams of the Arietids and 5-Aquarids are now 
different, it seems probable that they may have had a common origin in the past. 

(ii) ¢-Perseids and B-Taurids.—The mean orbits of the {-Perseid and 8-Taurid 
streams have inclinations of less than 6°, and therefore are sufficiently close to the 
ecliptic to produce another shower after passing perihelion. The expected 
radiant positions are found to agree with the night-time showers known to visual 
observers as the Arietids and Taurids (Table VI a and 5). 


TABLE VI 
(a) Radiant Positions of the Night-time Arietids 
Return Arietids 
of A . 
£-Perseids Norton (3) Whipple (11) Prentice (12) 
North South 
a (degrees) 33 42 42 41°6 38 
C) ie +10 +21 L419 +10°3 +12 
Mean date Oct.13 Oct. 12-23 Nov. 10 Oct. 20 Oct. 12 





(b) Radiant Positions of the Night-time Taurids 
Return Taurids 
of r a — 
B-Taurids Norton (3) Whipple (11) Hoffmeister (g) Prentice (12) 
North South 
sew 55 56°5 55°0 55 54 





” + 24 rts +22°5 14'S +24 +20 
Mean date Nov.g Nov. 3-15 Nov. 11 Nov. 9 Nov. 12 Nov. 6 


EARTH 





‘ 

\ SOUTHERN 
+ ARIETIDS 
*  Wo642 


C-PEASEIDS 
Fic. 3.—Orbits of (-Perseid meteor stream and an individual S. Arietid meteor. 


Whipple (13) has determined the velocities of meteors of these showers by 
photographic methods, and his individual orbits agree well with the mean orbits 
of the {-Perseid and 8-Taurid streams computed above (Table VII, Figs. 3 and 4). 

Whipple (11) has distinguished two radiants on opposite sides of the ecliptic 
for each of the visual showers, which he has called the northern and southern 
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streams respectively. ‘The radiant positions to be expected for the return showers 
of the mean orbits agree better with the S. Arietids and N. Taurids, although the 
radio-echo observations are not sufficiently sensitive to separate two such showers 
in the summer daytime, and the mean radiants may be a combination of positions 
above and below the ecliptic. 





a 


HOFFMEISTERS, 
TAURIDS : 


Fic. 4.—Orbits of B-Taurid and Taurid meteors (Nos. 778, 789 and 1009 are Whipple’s 
individual orbits for three Taurid meteors). 


TABLE VII 
(a) Orbits of €-Perseids and S. Arietids 
f-Perseids SS. Arietids No. 642 
Vy (km./sec.) 34°'8 36°3 
a (A.U.) z°9 1°9gI 
e 0°79 0'845 
q (A.v.) 0°35 0°296 
i (degrees) 4 6-1 
61 122°3 
138 150'2 
77°° 27°91 
(6) Orbits of B-Taurids and Taurids 
B-Taurids N. Taurids, Whipple (11) | Taurids, Hoffmeister (5) 
No. 778 No. 789 No. 1009 
V gq (km./sec.) 36: 37°1 36-9 37°1 35°11 
a (A.U.) ° 2°19 06 2°16 1°5738 
e "Bs 0891 ‘799 0858 0°7294 
q (A.U.) . 0239 "414 0°307 0°426 
i (degrees) 3°9 , 2°5 o'7 
307°7 ; 299'6 290°5 
160°1 : 160°5 160°5 
212°41 ‘ 220°92 230 
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As long ago as 1940 Whipple (1) predicted the return of the Taurids in the 
summer daytime, and he pointed out the similarity of the orbit to that of Encke’s 
comet. 

Hoffmeister (5) included the Taurids in his system of planetary streams and 
computed the orbit given in Table VII 6, and shown in Fig. 4. He estimated its 
return date to be June 22. 

(ili) o-Cetids.—T he o-Cetids stream observed in May has an inclination of 37°, 
and longitude of perihelion of 211°, and so does not approach the Earth at any 
other time of the year. Hoffmeister (§) suggested that the night-time Piscids 
should return on May 3, but these have not been observed as a permanent stream. 





~ ~ 


MINOR PLANET ADONIS ¢-PERSEIDS 


Fic. 5.—Orbits of ¢-Perseid meteor stream and minor planet Adonis. 


5. Conclusion.—The results described in this paper show that the major 
summer daytime streams all move in short-period orbits well inside that of 
Jupiter. 

Apart from Encke’s comet, which has been associated by Whipple (1) with the 
Taurid and f-Taurid streams, there are no similarly short-period cometary orbits. 
In fact, in size, shape and inclination the orbits resemble more closely some of the 
minor planet orbits, as can be seen in Fig. 5. It is interesting to compare these 
results with the theory proposed by Hoffmeister (5) of a system of planetary 
meteor streams in addition to the cometary streams. 
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THE ORIGIN OF GALACTIC RADIO-FREQUENCY RADIATION 
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Summary 


The results of observations of the intensity and distribution of 
radio-frequency radiation from the Galaxy at frequencies from 9°5 to 
3000 Mc./s. have been collected. Some of these data are used to determine 
spectrum curves of the radiation from chosen regions of the Galaxy. 

Using the equations defining the propagation of radio waves in an 
ionized gas, the forms of spectrum curves from distributions of gas are 
determined. It is shown that the power flux of radio energy from a 
thermally radiating gas cloud increases with frequency except under special 
conditions which are defined. 

From a comparison of the observed and the theoretical spectrum curves 
it is shown that galactic radiation at radio frequencies probably originates 
partly in hot ionized interstellar gas and partly in stellar atmospheres. The 
ionized gas provides most of the radiation at the higher frequencies and 
evidences itself by absorption at the lower. The properties of the stellar 
sources are discussed. 





1. Introduction.—Since Jansky (1, 2) discovered radio-frequency emission 
from the Galaxy many writers have speculated on its origin. Early theories 
relating to interstellar dust (3, 4) have been discarded in favour of theories of 
origin in interstellar gas or in stars. 

Emission by ionized interstellar gas has been discussed by Reber (5), Henyey 
and Keenan (6), van de Hulst (7), ‘Townes (8) and others. ‘The mechanism 
envisaged is that of free-free transition of electrons in the fields of positive ions, and 
treatments by classical and quantum-mechanical methods give consistent results. 
However, to explain experimental data, electron temperatures of the improbably 
high value of 10° deg. K. would be necessary. Furthermore, the distribution of 
galactic radio emission does not fit the accepted distribution of interstellar gas and, 
finally, such a theory would fail to explain radiation from “‘ discrete sources”’ 
found by Hey, Parsons and Phillips (g), Bolton and Stanley (xo) and others. 

The possibility of the radiation originating in stars has been discussed by Hey, 
Parsons and Phillips (11), Unséld (12) and others. The difficulty is to account 
for the very high intensity of observed radiation ; according to Unsdld this could 
by explained by an abundance of stars whose radio-frequency emission is about 
101 times that of the quiet Sun. 

In the present communication fresh observational data are included in an 
analysis covering the widest possible range of frequencies. ‘The radio-frequency 
spectra from different regions of the Galaxy are used to show that radiation 
originates partly in interstellar gas and partly in stellar atmospheres. The presence 


* Received in original form 1950 May 30. 
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of gas is demonstrated by both emission and absorption processes, the former being 
thermal in character. Emission from stellar sources is probably by a process other 
than thermal. 

2. The determination of spectra of galactic radiation.—The results of a number 
of measurements of the intensity of galactic radiation are now available. Some of 
these consist of surveys of the intensity of galactic radiation over large regions of 
the Galaxy. Such surveys have been made by Reber (13, 14) at frequencies of 
480 Me.'s. and 160 Me./s., by Allen and Gum (15) at 200 Mc./s., by Bolton and 
Westfold (16) at 100 Mc./s., by Hey, Parsons and Phillips (11) at 64 Mc./s. and 
by Shain (17) at 18:3 Mc./s. ‘T'wosurveys (at 100 Mc./s. and 64 Mc./s.) have been 
made by measuring absolute values of power received as the aerial beam moved over 
as much of the Galaxy as possible. From these data sets of contours of absolute 
intensity of the radiation were constructed by a process of trial and error. At 
frequencies above 100 Mc./s. the radiation from the “ cold’ parts of the Galaxy 
was too low to be measured. ‘Thus, the 200 Mc./s. contours indicate relative 
intensities above an unknown base level. A method of estimating this level is 
given in the Appendix. ‘The two sets of contours at 480 Mc./s. and 160 Mc.:s. 
have not been converted to intensities in a given direction ; they represent 
average intensities within the solid angle of the aerial beam. Once again they are 
only relative values above an unknown base level. In the case of the 18-3 Mc.:s. 
survey the data were limited to a series of measurements around the strip of the 
sky observed with the aerial beam pointed vertically (declination — 34°). The 
contours were then drawn by assuming that they would be similar in general form 
to the 100 Mc.,s. contours and by varying the contour intervals until the results 
were consistent with the experimental data. ‘The measurements were, however, 
absolute values. 

The results of the four surveys which consist of equal-intensity contours 
(allowance having been made for the aerial beam-width) are shown in a slightly 
modified form in Fig. 1. A general similarity between the four will be noticed. 
A similarity between the 18-3 Mc./s. and 100 Mc./s. sets may have been anticipated 
from the method of derivation; however, it should be remembered that unless the 
sets were generally very similar (at least in the regions explored at 18-3 Mc./s.) it 
would have been impossible to fit the experimental data to the assumed contour 
shapes. 

Of particular interest in Fig. 1 is the fact that the distance between successive 
contours is smaller at the higher frequencies. This will be discussed in more 
detail later. 

In order to investigate possible sources of galactic radiation, not only the 
contour sets of Fig. 1 but also detailed emission spectra for representative parts of 
the Galaxy are required. If contour charts of absolute intensities existed for all 
frequencies at which measurements have been made, such spectra could easily be 
prepared. Unfortunately, as indicated above, some of the data require careful 
interpretation. In the first place the base level of intensity in the ‘‘ cold” parts 
of the Galaxy must be determined for frequencies above 100 Mc./s. Secondly, in 
order to make them comparable on different frequencies all data must be reduced 
either to intensity in a given direction or to power received by an aerial of given 
beam-width, pointed in a given direction. ‘The data used by ‘Townes (8) in his 
determination of a galactic emission spectrum were not reduced to correspond to a 
uniform beam-width and for that reason his results differ somewhat from our own. 
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It was decided to determine the spectra of three regions A, B and C, in terms of 
the power received by an aerial of beam-width + 3° to half-power*, pointed in the 
three directions : 

A—Galactic latitude — 2°, galactic longitude 330°. 
B—Galactic latitude 0°, galactic longitude 180°. 
C—Galactic latitude — 30°, galactic longitude 200°. 


Region A was chosen because it lies near the centre of the Galaxy, which is the 
only place where measurements at the highest frequencies have been made. 
Region B is in the Milky Way but in a region where radiation is much less intense 
than at A. Region C is well away from the Milky Way in an area where the con- 
tours are widely spaced, so that the actual size and position of the region is not 
critical. ‘The positions of A, B and C are shown in Fig. 1 (18-3 Mc./s. contours). 

At the low frequencies concerned the Rayleigh-Jeans law provides a linear 
relationship between the intensity of radiation emitted by a black body and its 
temperature. ‘The intensity of radiation from a given direction may therefore 
be simply specified in terms of the temperature of a black body. Similarly the 
power received by a given aerial may be specified in terms of the temperature of a 
black body which entirely fills the aerial beam and gives the same power input. 
This temperature is called the (average) equivalent temperature of the region 
covered by the aerial beam. ‘The equivalent temperatures of the regions A, B and 
C will be referred to as T',, 7, and T;,, respectively. 


TABLE I 


Frequency Equivalent ra i 
Observer f temperature deg. K. (Te=5000deg. KK.) ¢ 


(Mc./s.) aa Tp T 

Friis and Feldman (18) 9°5 = 2"4 = 10° fe 
Shain (17) 18-3. 2°0* 10° 75000 50000 
Moxon (19) 40 67000 11900 8500 
Hey, Parsons and 

Phillips (11) 64 21000 3100 2200 
Moxon (19) go 7700 
Bolton and Westfold 

(16) 100 6000 
Reber (14) (modified) 160 2180 ; a : 4°16 
Allen and Gum (15) 200 1190 . 3°62 
Reber (13) (modified) 480 145 os : 2°72 
Piddington and 

Minnett (20) 1200 17°9 one we . 2°26 
Piddington and 

Minnett (20) 3000 «ee _— 105 2°36 


The method of finding 7,, 7, and 7;, from the various data available is 
described in the Appendix. Values are listed in Table I together with the reference 
source. The figures listed in the last three columns will be explained later. It 
should be noted that all data for 200 Mc./s. and above have been adjusted by adding 
an estimated “‘ base level”. ‘This correction is about 5 per cent of T, and the 
accuracy of determination is discussed in the Appendix. Experimental data at 
18-3 Mc./s. exist for regions A and C but not for region B, where conditions could 
only be estimated from the assumed similarity of contours discussed above. 

Although some results are given to three significant figures it is unlikely that 
the accuracy is anywhere better than 10 per cent: it will usually be nearer 
20 per cent. However, the best determination is quoted in each case. 


* The shapes of this and other aerial beams are discussed in the Appendix. 
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’ 

Values of 7’, and 7’, are plotted against frequency in Fig. 2 and the two curves 
marked 7, and 7, are drawn throughthem. ‘The points corresponding to Reber’s 
160 and 480 Mc./s. data and those of Moxon at 40 Mc./s and go Me./s., and Friis 
and Feldman at 9-5 Mc./s., are marked as circles (all*others being disks) since they 
are considered less reliable, for reasons given in the Appendix. A plot of 7, is 
omitted because of its similarity to 7),. The slope of the spectrum curves is of 
principal interest, and is expressed by 7'«f-*, where T corresponds to any portion 
of the 7, or 7, curves, f is the frequency and « a constant for the portion of the 
curve in question. « hasa value of 2-73 for the central portion of the 7’, curve and 
a smaller value (2°51) for the same part of the 7’, curve. At frequencies above 
1000 Me. s. « falls to approximately 2-0, while at the low-frequency end of the 
7, curve x becomes nearly zero. 












































Equivalent temperatures ef regions A and B plotted against frequency. 


Ordinates: Equivalent temperature (deg. K.). 
Abscissae: Frequency (Mc./s.). 


It is of interest to compare these conclusions with those of Herbstreit and Johler 
(21), who measured the relative amounts of radiation received by wide-beamed 
aerials at 110 Me./s. and 25 Mc./s. ‘They found a value of « = 2-41, which may be 
compared with values of « over the same frequency range of 2-00 and 257 for regions 
A and B respectively. As the effect of the 7), type spectrum would be relatively 
more important the agreement is satisfactory. 

Three further sets of galactic radiation measurements have been omitted 
from Table I but should be mentioned. Sander (22) used a rather wide beam at 
60 Me./s. His data may be compared with those of Hey et al. (11) by the graphical 
integration method given in the Appendix. ‘The results are about ro per cent 
smaller than Hey’s instead of being slightly higher (because of the lower frequency). 
Franz (23), observing at 30 Mc./s. with a beam of +15°, found an equivalent 
temperature near the galactic centre of about 1-2 x 10°deg.K. This value is 
very close to 7’, given by the curve of Fig. 2, which would be expected if the contour 
lines were widely spaced as for 18-3 Mc./s. (Fig. 1). Observations by Jansky 
(2,24) at 18 Mc./s. with a wide beam have been analysed by ‘Townes (8), who found 
an equivalent temperature of 9:2 x 10* deg. K., which is slightly lower than the 
value of 7, given in Table]. Each of these results shows fair agreement but 
was omitted because other results on nearby frequencies were considered more 
reliable. 
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3. Theoretical radio-frequency emission spectra.—\n considering the theoretical 
emission spectra of various components of the Galaxy the following theorem 
is important. Consider a distribution of ionized gas in which the temperature, 
pressure, degree of ionization and chemical constitution have any desired distri- 
bution in space. ‘Then, unless one of the special conditions listed below obtains, 
the radio spectrum of the gas is such that the power flux of radiation never decreases 
with increasing frequency. In general, the power flux increases with frequency 
but in the limiting case of an optically very thin radiator it remains constant. 
‘The special conditions are : 

(a) A magnetic field of sufficient intensity exists within the gas. 

(6) ‘The dielectric constant of the medium (at the frequencies concerned) 
departs appreciably from unity. 

(c) ‘The emission is non-thermal. 

A fourth possibility exists : that emission is by the atomic or molecular processes 
discussed by Saha (25) and not by collision of electrons with positive ions, atoms or 
molecules. Such emission would tend to provide irregularities in the spectrum 
and not to determine its general form over a wide frequency range. It is probably 
not important in the present discussion and will not be considered further. 

In a distribution of ionized gas let us consider radiation issuing axially from a 
cylindrical element having unit area, thickness ds, optical thickness dz and situated 
at optical depth r. ‘The intensity of the radiation is proportional to dr and also, 
according to the Rayleigh-Jeans formula, to the square of the (radio) frequency /. 
Absorption is suffered as the radiation leaves the gas cloud so that the contribution 
to the power flux at any point outside the cloud ts given by 


dW xf*dre’, (1) 
the constant of proportionality being independent of frequency. 
If the medium is free from magnetic fields then radio waves are propagated 
according to the formula (see, for example, Booker (26) ) 


( ick \? Ne* 
b _ os I a a * 

27 , iv 

f amf{?{1— — 
; anf 

where , is the dielectric constant which we take as unity, « is the absorption 
coefficient, c the velocity of light and N, «, mand v the electron density, charge, mass 
and collision frequency respectively. Expanding, collecting imaginary terms and 
writing dz =«ds (by definition) 


It is evident, from an examination of (2), that the values of both f* dr and e-’ 
will increase as f increases. Thus, from (1), dW will increase with frequency and 
hence the power flux from the whole gas cloud will increase with frequency ; which 
proves the theorem. 

A special case of particular interest is that of an optically thin cloud of gas. 
This condition implies (unless the cloud has negligible dimensions) that v<f. 
Equations (1) and (2) then reduce to 

sa 
Z2rme 
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so that the power flux is independent of frequency. This case, which is the 
important one of interstellar gas (at some frequencies), has already been considered 
by Townes (8). 

The value of the absorption coefficient of an electron-proton gas may be derived 
by several methods, which have been summarized by the author (27) and found to 
show good agreement. We have 

CN* 


— 


PT 2?’ 


where 7, is the electron temperature and ¢ a slow variable which differs slightly for 
different methods of deriving x. If the value of « is found from Gaunt’s (28) 
formula for free-free emission of an electron in the field of a proton, then 


€=9'70 x 1073 log, (=F) ; 


(4) 


where k and A are Boltzmann’s and Planck’s constants and units are C.g.s. 
Comparing equations (2), (3) and (4) it is evident that v is not a constant but a 
slowly varying function of frequency. ‘This may appear strange, but is due to the 
fact that equation (2) does not take account of the complex mechanism of collision 
between a proton and electron. 
The equivalent temperature of an optically thin cloud of gas is given by 

T =f Txds, 
integration being along the line of sight. If « is given by equation (3) then 

nt NP 

T= | —, ds. 

f°) fe" 

‘The expression under the integral sign is independent of frequency so that the 
spectrum curve is defined by 


Tf?t-1 =const. 


The above argument which referred to an electron-proton gas is applicable to any 
ionized gas, the slope of the 7f*¢~! curve being zero. 

We conclude, therefore, that if the value of 7f?¢-! for radiation from any distri- 
bution of gas is plotted against f, then the slope is, in general, positive. A limiting 
case for optically very thin gas shows a zero slope while a negative slope indicates 
that one of the special conditions listed above must hold. Some illustrative 
spectra are plottedin Fig. 3. The first, curve (a), corresponds to emission by an 
optically thin volume of gas. Curve (4) is a typical spectrum of thermal radiation 
from a gas in which magnetic fields play no part and the dielectric constant does 
not depart appreciably from unity. It is the spectrum curve of radiation from the 
quiet Sun. It was derived from experimental data by the author (27), who showed 
that, at least for frequencies above 600 Mc./s., the refractive index did not differ 
appreciably from unity.* An example of a radio spectrum curve in which Tf?¢-! 
decreases with increasing frequency is shown as curve (c) of Fig. 3. This represents 
emission from the disturbed Sun in the range 30 Mc./s. to 200 Mc./s. ona particular 
occasion, and is plotted from a curve given by Appleton and Hey (29) obtained from 


* The ordinate scale is arbitrary and differs for curves (b) and (c) so that they may be displayed 
together. 
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experimental results. It is not known to what extent the source size influenced 
this apparent spectrum, but there is little doubt that on occasions radiation from 
the disturbed Sun does have a spectrum curve with negative slope. 

It is desirable to consider the special conditions listed above for which the 
spectrum curve has a negative slope. ‘These are : 

(a) ‘Thermal emission in the presence of a magnetic field. 

An examination of Appleton’s magneto-ionic equations indicates that to affect 
propagation appreciably the strength of the field would have to approach the gyro 
value, which for a frequency of 100 Mc. s. and electrons is 36 gauss. Under these 
conditions the spectrum might take almost any form. 

(6) ‘Thermal emission when « departs appreciably from unity. 

Small changes in ,« may cause bending of the rays leaving the source and may 
result in a change in the amount of radiation leaving in a given direction. Since 
galactic radiation probably has a large number of randomly orientated sources the 
averaging effect should overcome this bending effect. However, if u varies by large 
amounts, any radio spectrum may result. For example, consider radiation from 
the ionosphere in the absence of galactic radiation. Waves of certain frequencies 
emitted from the Earth are reflected and partially absorbed. Hence by Kirchhoff’s 
law thermal radiation is emitted by the ionosphere at these frequencies. At higher 
frequencies the ionosphere becomes practically transparent and hence does not 
emit, so that the slope of the 7f*f-' curve is negative. ‘The electron density 
necessary to cause these effects is considerable, of the order 10° em. * for a frequency 
of 100 Me. s. 
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Fic. 3.—Some theoretical radio-frequency emisston spectra. 


Tf? 


Ordinates: ; ak 


Abscissae: Frequency (Mc.|s.). 


(c) Non-thermal emission. 


Emission by particles moving in an ordered manner has been discussed by 
Martyn (30), Haeff (31), Bailey (32) and others. Because of the complexity of the 
problem almost any spectrum might be emitted. 

4. A review of the experimental data.—’The experimental data to be reviewed 
consist of the four sets of equal-intensity contours shown in Fig. 1 and the spectra 
for regions A and B shown in Fig. 2. 

The various sets of equal-intensity contours published have been drawn with 
contour intervals chosen arbitrarily. ‘Thus successive contour lines of Hey et al. 
(11) are marked (in their units) 2°5, 3-0, 3:5, 4°0, 4°5, 5°O, 7°5, 10, 12°5, 15, 17°5, 
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20, 22:5 and 25. ‘Those of Bolton and Westfold (16) are marked (in their units) 
5, 6°25, 7°5, 10, 12°5, 15, 20, 30, 40, 50 and 60. Since intervals between successive 
contour lines vary arbitrarily it is impossible to make a fair visual comparison 
between different sets of contours. In Fig. 1 this difficulty is overcome as follows : 
for each set of contours a minimum intensity of radiation is determined, corre- 
sponding approximately to the intensity from region C. This intensity determines 
the first contour line (of lowest intensity) and also the interval between every 
successive pair of contour lines. ‘These equal-interval contour lines are drawn by 
interpolating between the published lines. In addition to equalizing contour 
intervals the galactic longitude scale is made identical for all four sets of contours ; 
this allows easy comparison of contour shapes. 

In comparing the four sets of contours, minor (and for the present purposes, 
unimportant) differences may be noted. ‘Thus the 200 Mc./s, chart has amaximum 
at latitude 0°, longitude 230°. ‘This is missing from the 100 Me. s. chart. 
The 64 Mc./s. chart has a large central piece missing because this was not visible 
from the northern hemisphere where the observations were made. In addition it 
has several small maxima around the point latitude — 3°, longitude 340 , and no 
maximum at latitude 0°, longitude 150°. At the latitude at which the observations 
were made the regions concerned always lay near or below the horizon, making 
observations difficult; the effects may be experimental errors resulting from such 
difficulties, 

The features about Fig. 1 which are striking and important are, first, that there 
is a marked similarity between all four sets of contours* and, second, that the 
distance between contours, particularly near the galactic centre, decreases steadily 
as the frequency rises. ‘Thus it may be seen that the total number of contour 
lines at the four frequencies 18-5 Mc./s., 64 Mc./s., 100 Me./s. and 200 Me. s. are 
4, 9, 12 and 23 respectively. ‘These features are also found in the sets of contours 
published by Reber (13, 14) for frequencies of 480 Mc./s. and 160 Me.:s. 

Let us now examine the features of the spectrum curves for regions A and B 
shown in Fig. 2. In order to compare them with the theoretical curves of Fig. 3 it 
is desirable to replot them as 7f*{~! against frequency. ‘The values of ¢ (from 
equation (4)) and of 7, f?¢-' and 7), f?¢—" are listed in the last three columns of 
Table I and the spectra are plotted in Fig. 4. 

Comparing these spectra with the theoretical curves of Fig. 2, the following 
features may be noted : 

(a) In the frequency range 3000 Mc.)s. to a few hundred megacycles per 
second the spectrum of region A is constant within the uncertainty of observation, 
and hence consistent with an origin in optically thin, thermally radiating gas. 
It should be realized, however, that experimental errors, which may be as great 
as about 20 per cent, may necessitate a revision of this part of the spectrum. 

(6) Inthe frequency range 400 Mc./s. to 40 Mc./s. the spectrum curves of both 
regions have negative slopes similar to that of curve (c), Fig. 2, and are therefore 
inconsistent with a thermal origin in a gas of low density (such as interstellar gas). 
They are, however, consistent with a non-thermal origin or a thermal origin in the 
presence of a strong magnetic field or high electron density. Because of the 
unanimity of results, it is unlikely that this conclusion could be invalidated. 


* The method of derivation of the 18-3 Mc./s. contours, discussed in Section 2, somewhat 
reduces their value as evidence in this respect. 
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(c) Over most of the frequency range 400 Mc./s. to 40 Mc./s. the relativerate of 
change of 7),f*¢-! is greater than that of T, f?¢-'. That is, for regions where the 
equivalent temperature T is low, 

1|dT 

T | df 
is greater than for regions where 7 is high. An inspection of Fig. 1 shows that this 
is not a peculiarity of regions A and B but applies in general; it is evident from the 
fact that the intensity ratio (or number of contour lines) between two given regions 
increases with frequency. 
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Fic. 4.—Spectra of regions A and B plotted as Tf*¢ 
Ordinates: Tf?! for regions A and B. 
Abscissae: Frequency (Mc./s.). 


(d) Inthe range 40 Mc./s. to 18-3 Mc./s. the slope of the 7’, {*¢-! curve changes 
rapidly to a positive value while that for 7), f°¢ ' remains negative. 

5. The sources of galactic radiation.—All the above peculiarities of galactic 
radio-frequency spectra find a ready explanation in terms of emission by stellar 
atmospheres and both emission and absorption by interstellar gas. 
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‘The manner in which the spectrum varies according to the distance from the 


galactic plane as indicated by Fig. 1 might, at first sight, be explained in one of 
the following three ways :— 


(a) Emission is from a single type of source whose properties vary according to 
the distance from the galactic plane. 


(6) Radiation emitted from a single type of source is subsequently partially 
absorbed by an amount which varies with galactic latitude. 


(c) Emission is from two or more types of source having different distributions 
in galactic latitude. 


The first possibility may be ruled out as most improbable, particularly as the 
various phenomena have simple alternative explanations. 

Turning to the second possibility, it may be noted that the radiating medium 
could not be interstellar gas radiating thermally. In the first place the electron 
temperatures necessary to explain the radio data would have to be at least as high 
as the highest observed equivalent black-body temperat: ¢, namely 2 x 10° deg. K. 
Spectroscopic evidence indicates that interstellar gas, e. ca in the vicinity of very 
hot stars, is maintained at a temperature of about 10'deg. K. The shape of the 
radio spectrum provides additional evidence that the radiation cannot originate in 
interstellar gas. As shown in Section 3 above, a negative slope of the 7f?¢-! curve 
indicates a strong magnetic field, a high electron density or non-thermal emission. 
Itis unlikely that interstellar gas could exhibit any of these properties ; we conclude, 
therefore, that if radiation is from a single type of source, then the source cannot be 
interstellar gas and is probably stellar atmospheres. ‘The absorption effect, on 
the other hand, could only be caused by ionized interstellar gas since this is the only 
medium which is extensive enough. Further, if ionized gas absorbs, it will also 
emit radiation, and since its temperature is about 104 deg. K. such emission might 
be of appreciable intensity. A subsequent quantitative treatment shows that 
this is so. ‘Thus the second possibility listed above really implies radiation 
from two sources, one of which is interstellar gas. 

Let us now consider the third possible combination of sources listed above. 
Stellar atmospheres and interstellar gas may be the two types of source or, 
alternatively, emission may be from two or more types of stellar atmosphere 
existing in different proportions throughout the Galaxy. While the second 
possibility cannot be ruled out, it appears improbable in view of the simple 
explanation of all the spectral peculiarities in terms of radiation by stellar 
atmospheres and interstellar gas and absorption by the latter. In a communica- 
tion in preparation, further evidence is provided by a quantitative discussion of the 
spectra. We conclude that the observed galactic radio-frequency radiation 
originates principally in interstellar gas and the atmospheres of stars. 

The spectrum curves of Fig. 4 may now be explained in terms of these sources. 
In the frequency range above about 500 Mc./s. the spectrum curve of region A is 
similar to that of optically thin gas. Evidently the stellar radiation which pre- 
dominates at lower frequencies but which decreases rapidly at higher frequencies is 
negligible here. At frequencies below about 100 Mc./s. the slope of the Tf?¢-! 
curves indicates radiation from stellar atmospheres. The fact that the slope is 
steeper for regions away from the galactic plane indicates that absorption by inter- 
stellar gas (which is highly concentrated near the galactic plane) is increasing in 
effect. Finally the marked change of slope for region A and frequencies below 
about 60 Mc./s. indicates strong absorption by gas. 
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6. Interstellar gas and stellar sources.—In the previous section we have seen 
that the spectra of galactic radiation provide evidence of the sources of such 
radiation. It is possible to proceed a step further and make quantitative deter- 
minations of the electron temperature and density and other properties of the 
radiating gas. For example, the shape of the spectrum curve for region A (Fig. 4) 
indicates that radiation at 3000 Mc./s. is probably thermal emission from optically 
thin gas. If zis very small, then equation (3) gives 

T= T,7, 
where 7 is 2-77 deg. K. and 7, is probably about 10‘ deg. K., so that 


rT=28xi107%, 


Again at 64 Mc.,s., since 7 <f-* the optical depth would be 0-61. ‘If the radiating 
(and absorbing) gas were uniformly distributed along the line of sight, this large 
value of optical depth would cause much more absorption at the low frequencies 
than is indicated by the spectrum curves. ‘The conclusion is that the gas exists 
in the form of relatively small clouds. ‘The properties of these clouds and of the 
stellar atmospheres are the subject of a further communication in which a model 
radio galaxy is determined. 

We have seen that the emission from stellar atmospheres may be thermal and 
controlled by magnetic fields or high electron densities, or it may be non-thermal. 
Ryle (33) has suggested a thermal origin, but in view of the very high equivalent 
temperatures observed at the lower frequencies this appears most improbable. 
Thus Shain’s (17) results at 18-3 Mc./s. indicate an average equivalent temperature 
over the celestial sphere of about 6 x 104 deg. K., and if Dunham’s (34) estimate of 
10 for the dilution of all stars (the proportion of the celestial sphere filled by 
stellar disks) is adopted, then an electron temperature of about 6 x 10! deg. K. 
would be necessary. Not only are there great difficulties in accounting for the 
continued existence of such an enormous temperature, but the gas would not 
radiate appreciable radio energy because of the very small collision cross-section of 
protons (or other nuclei) for electrons. ‘The conclusion is, therefore, that emission 
is by a much more efficient (non-thermal) process. 

It has been assumed that all radiation discussed above was galactic in origin. 
At least two of the “‘ discrete sources ’’ of radiation mentioned earlier have been 
identified (35) with what may be extra-galactic nebulae. However, the manner in 
which the contours of Fig. 1 follow the general form of the Galaxy makes it highly 
improbable that any but a small proportion of the total radiation is extra-galactic in 
origin. Confirmatory evidence is provided by a comparison between the radio 
data and the distribution of observed extra-galactic nebulae (36). No corre- 
spondence is observable. 

7. Acknowledgments.—The author would like to thank Dr J. L. Pawsey and 
Mr B. Y. Mills for useful discussions during the progress of this work. 


APPENDIX 
In determining 7, 7, and 7,. from the various published data two processes 
may be necessary. ; 
(i) Conversion to equivalent aerial temperatures of a standard + 3° beam. 
The standard is shown in Fig. 54 ; it is the average (in shape) of two beams 
used by Piddington and Minnett (20) and should not differ appreciably in shape 
from those used by other observers when all are scaled down to + 3”. 





Jj. H. Piddington Vol. 111 


When the observational data are presented in the form of intensity contours, 
T,, Ty, and 7,,may be found by a straightforward method of graphical integration, 
account being taken of the aerial gain and the intensity of radiation in every 
direction. Where data are given as observed power flux of radiation per unit 
solid angle, these are converted to equivalent temperatures by means of well- 
known formulae. 

At some frequencies the only available data are values of equivalent aerial 
temperatures in a limited number of directions. Where a set of contours is 
available on a nearby frequency this is considered to apply (with all values of 
intensity altered in a fixed ratio) for the data in question. For example, the values 
of 7,, 7), and T,, at 90 Mc./s. are found from Moxon’s data by using the 100 Me. s. 
contours. Where contours are available at higher and lower frequencies two 
determinations of 7, 7, and 7;, are made and the average used. For example, 
the 18-3 Mc./s. and 64 Mc.’s. contours are used to find 7,, 7;, and 7,. from 
Moxon’s 40 Mc./s. data. 2 
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At frequencies above 200 Mc./s. the correction for beam-width is more difficult, 
since there are no sets of intensity contours available on nearby frequencies. 
However, as will be shown below, the form of the contours, even near the galactic 
centre, does not change greatly above 200 Mc..s., so that the results of Allen and 
Gum may be used to determine 7',, 7’, and 7, at higher frequencies. Given the 
equivalent aerial temperature for any beam pointed at the centre of region A the 
value of 7, may be found directly from Fig. 5 6, which gives relative responses for 
a wide range of beam-widths. It was derived from the 200 Mc./s. contours by a 
series of graphical integrations.* 

* It may be noted that the response indicated by Fig. 5 6 for Allen and Gum’s beam of 
+12°5° is higher than the value actually measured. The apparent discrepancy is explained by 


the fact that only about 75 per cent (a normal value) of the total aerial response is in the main beam, 
the remainder being dissipated in side lobes. 
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The values of 7, and 7,, may be found similarly, the integration process being 
much simpler. As there probably is some decrease in contour spacing with 
frequency the values so found may have small errors. 

(ii) Base-level determinations. 

For each frequency above 100 Mc./s. the data available are intensities above an 
unknown base level. The base-level intensity, corresponding to regions well 
away from the galactic plane, is only a small proportion of 7',, but an estimate is 
desirable and may be found from the following considerations. 

An examination of the various published galactic surveys suggests that near 
the galactic plane the contours become more crowded as the frequency increases. 
However, lack of knowledge of base-level intensities at some frequencies warrants a 
closer examination, and for this purpose Fig. 6 is drawn, showing relative intensities 
at the various frequencies plotted against displacement from the centre of region A 
along the line of declination —27°. Each intensity value is the mean of the two 
values on either side of region A. ‘The data are obtained from the sources listed 
in ‘Table I. 
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Plots of relative equivalent temperatures along a line (dec. — 27° through the 
galactic centre. 


Ordinates: Relative equivalent temperatures. 
Abscissae: Angular displacement from centre of region A (degrees). 
Note: Base levels omitted from dashed curves. 


Before comparing these curves, account must be taken of two factors. First, 
the curves for 160 Mc.)s., 480 Mc./s. and 1200 Mc./s. show relative equivalent 
aerial temperatures for the aerials used. If allowance is made for beam-width then 
the curves would be steeper (compare, for example, the 100 Mc. s. contours (16) 
before and after correction for beam-width). ‘The effect on the 160 Mc./s. curve, 
which applies for a beam of about + 6°, would be considerable. ‘The curves for 
480 Mc./s. and 1200 Mc./s., which apply for beams of about +2°-5 and +1°-7 
respectively, would be little changed. 

Secondly, the base level of radiation must be added to the curves for frequencies 
above 100 Mc./s. This level has fallen from 25 per cent of the maximum intensity 
at 18-3 Mc./s. (17) to about 8 per cent at 100 Mc./s. and is hardly likely to rise again 
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for higher frequencies. Even if the proportion remains unchanged for higher 
frequencies, there is obviously a marked increase in slope of the curves from 
18-3 Mc./s.* to 200 Mc./s. This result, illustrated in a different way in Fig. 1, 
is vital to the whole of the theory developed in this paper. Its implications have 
been discussed in Section 5. 

In contrast to the marked increase in slope between 18-3 Mc./s. and 200 Mc./s., 
there appears to be only a small increase for higher frequencies (even when 
allowance is made for aerial beam-widths). ‘This justifies the approximate method 
of beam-width correction used above. Also, an inspection of the various sets of 
contours indicates that the spacing of contours in the cold parts of the Galaxy 
changes very little throughout the whole frequency range. It is concluded that 
the change for two reasonably adjacent frequencies (say 100 Mc./s. and 200 Mc/s.) 
is negligible; this provides a method of determining the base level of intensity at 
frequencies above 100 Mc./s. 

Determination of 7',, 7, and 7, for frequencies from 9-5 Mc./s. to 3000 Mc.’s. 
is made as follows: 

(a) Friis and Feldman (18) have measured galactic radiation at 9:5 Mc./s. 
which ‘Townes (8) has interpreted as indicating an equivalent temperature of 
1-2 x 10°deg. K. Reber and Greenstein (37) have estimated the direction of the 
beam as being in the constellation of Cygnus, so that 7, will be somewhat higher. 
From Shain’s contours an approximate value of 2-4 x 10°deg. K. is deduced, 
although not much reliance should be placed on this figure, as the “ discrete 
source’ in Cygnus (10) might have influenced the result, and also the spacing 
between contours may have decreased between 18-3 Mc./s. and 9:5 Mc./s. 

(6) The derivation of Shain’s (17) 18-3 Mc./s. contours from rather limited 
data has been discussed. The values of T,, 7, and T,,can be made directly from 
the contours. 

(c) Using a rather wide aerial beam, Moxon (19) measured increases in aerial 
temperature at 40 Mc./s. as the beam crossed the Milky Way. He then attempted 
to reduce these to a plot of excess equivalent temperature, 7, along the galactic 
equator. ‘The average value of T in the region of galactic longitude 140°-220° is 
3°4 x 10°deg. K. From the contours of Hey et al., or Shain, it is evident that the 
average equivalent temperature in this region is approximately equalto 7. Thus, 
if T,. is taken as the minimum equivalent temperature (well away from the galactic 
plane), we have 

T, — Te = 3°4 x 108 deg. K. 
Also the ratio 7),/ 7; at 40 Mc./s. should lie between the ratios at 64 Mc./s. and 
18-3 Mc./s. which are I-41 and 1-39 respectively. Putting 7/7, = 1-40 we find 


T, = 11900 deg. K., 
T,. = 8500 deg. K. 


In the region of galactic longitude 320°-340° Moxon finds a value of T of 
about 3:1 x 104 deg. K. In this region the contours are closely spaced, and as the 
results apply to an “ energy width ” of about 35°, a beam-width correction factor 
must be found by graphical integration. Factors of 2-2 and 1-2 were found for 

* The derivation of the 18-3 Mc./s. contours depended on an assumed general similarity with 
those at 100 Mc./s. However, the observations were made along a line of constant declination, 


only a few degrees from the line for which Fig. 6 is drawn. Thus the curve shown appears to be 
supported directly by the experimental evidence. 
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frequencies of 64 Mc./s. and 18-2 Mc./s. respectively. Using a mean value of 1-7 
and adding 7, we find 

T, =6°7 x 104 deg. K. 
Not much reliance can be placed on the accuracy of the above determinations. 

(d) Values of T,, T, and 7;, at 64 Mc./s. may be found directly from the 
contours of Hey, Parsons and Phillips (11). 

(e) Moxon’s (19) observations at 90 Mc./s. with a +17:5° beam show an 
equivalent temperature near the galactic centre of 3500 deg. K. Using Bolton and 
Westfold’s contours to correct to a + 3° beam the factor is approximately 2-2, so 
that 7, =7700 deg. K. 

(f) The contours of Bolton and Westfold (16) at 100 Mc.:s. allow values of 
T,, T), and T,, to be determined directly. 

(g) Reber’s results at 160 Mc./s. indicate an equivalent temperature of 
5900 deg. K. for a beam of width + 6° (estimated from published aerial dimen- 
sions and frequency). ‘The correction factor to a standard aerial beam is 1-4 
at 200 Mc./s. and 1-3 at 100 Mc./s.; a value of 1-35 is used. The base-level cor- 
rection factor is about 1-06 for the 200 Mc./s. contours (as determined below) 
and 1-08 for the 100 Mc./s. contours; a value of 1-07 is used, leading to a value for 
T, of 8500 deg. K. However, ‘Townes (8) and Thomas and Burgess (38) have 
shown that Reber’s results are too high by a factor of about 3-9. The very high 
values of solar emission published by Reber indicate a similar factor. ‘The cor- 
rected value of 7’, is 2180 deg. K., but little weight should be given to this figure ; 
values of 7, and 7;, are not estimated. 

(h) Before Allen and Gum’s (1§) intensity contours can be used it is necessary 
to determine a base level. ‘They mention a value of roo deg. K., suggested by 
Bolton, who extrapolated data from lower frequencies. Since an absolute 
determination of intensity may be in error by as much as 20 per cent, the extra- 
polation method may result in an error of 40 per cent or more. Fortunately, as 
indicated above, the base level may be found by comparing the 100 Mce.’s. 
(absolute) and the 200 Mc./s. (relative) contours on the assumption that the distri- 
butions in the cold parts of the sky are similar. In principle we may compare the 
effective temperatures, at the two frequencies of the region C, 7.100 and 7;.200, 
with those of an adjacent region D. ‘Then from our assumption 


T,.100 =T,.200 





Ty100 = T2007 


Since 7,,100, T)100 and (7,200 — T,.200) are known, the values of the unknown 
base level at 200 Mc./s. may be derived. If our assumption is not quite accurate it 
is likely, as shown above, that 

T,.100 ~=—T,,.200 


‘sy al ’ 
T,100 ~=7))200 





so that the derived value of the base level is likely to be too high rather than too 
low. 

The use of single points C and D introduces the possibility of errors due to 
local irregularities in one or other of the sets of contours. In practice two complete 
contours at 100 Mc./s. were compared with two at 200 Mc./s. Suitable contours at 
200 Mc./s. were the minimum level x (in units of 10~*! watts m.-*(c./s.)~! sterad.~} 
as used by Allen and Gum)and the 1-0 unit contour. The corresponding 100 Mc./s. 
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contours are 500 deg. K. and 1072 deg. K., the latter being found by interpolating 
between published 100 Mc./s. contours for a large number of points on the 
I-O unit 200 Mc./s. contour line. ‘The value of x is found to be 0-87 units. 
A similar determination, using the 200 Mc./s. 0-4 contour line, gives x =0-85. 
‘The mean value of 0-86, corresponding to an equivalent temperature of 71 deg. K., 
is adopted. ‘The values of 7',, 7, and 7, may then be found directly from Allen 
and Gum’s survey, allowance for the beam-width being made when determining 7°,. 

The accuracy of the above determination of base level may be estimated as 
follows: Bolton and Westfold quote an accuracy for the cold parts of the Galaxy of 
4 per cent, so that the error in the left-hand side of the equation should be less than 
8 per cent, and probably much less because a large number of points are averaged. 
Allen and Gum quote an accuracy for most of their diagram of 20 per cent but the 
process of integration around the whole of a contour would lead to a much more 
accurate (relative) value; probably better than 10 per cent. Substituting these 
errors in the equation an uncertainty in x of 18 per cent is found. ‘The close 
correspondence of the two determinations, however, suggests a greater certainty. 

Moxon’s measurements of galactic radiation at 200 Mc./s. cannot be interpreted, 
as the aerial beam-width is not specified. 

(z) At 480 Mc./s. Reber’s (13) contours indicate a value of 7, of 250 deg. K. 
To this should be added a zero-level correction which might be taken (following 
Allen and Gum) as about 5 per cent of 7,. ‘The corrected value of 7, would 
be 262deg.K. ‘There appears, however, to be an error in the estimated 
intensities which may be traced as follows. In an earlier publication Reber (39) 
measured solar radiation over a period of several months and found it to be 
5 x 10 *° watts cm.-*(c./s.)-', which is equivalent to a black-body temperature 
of 1-05 x 10° deg. K. for a disk subtending the same diameter as the Sun. ‘The 
observed fluctuations were + 74 per cent, so that minimum temperature measured 
was 9°7 x 10°deg. K. A recent investigation by the author (27), however, has 
shown that the anticipated minimum disk temperature at this frequency is about 
5°5 x 10°deg. KK. ‘Therefore, unless the discrepancy can be accounted for by 
an error in the aerial gain, it would seem that Reber’s solar, and therefore cosmic, 
results should be reduced by a factor of about 1-8. Although both sets of measure- 
ments were part of thesame series and should, therefore, have been self-consistent, 
little weight will be given to the corrected result. Including the correction, 


T, =145 deg. K. 
Similarly, 
T), = 16°60. 

(j) At 2100 Me./s. Piddington and Minnett (20), using a + 2°-85 beam found a 
difference in equivalent temperature of 15-3 deg. K. between two regions centred 
at 1=330°, b= —2° and ]1=340°, b=—17°. The beam-width correction 
(Fig. 54) is 0-98 and the correction for radiation from the cooler region (from 
Allen and Gum's contours) is I-19, so that T, = 17-9 deg. K. 

Using a +1°-7 beam Piddington and Minnett also determined the curve of 
intensity along a line of constant (— 27°) declination as shown in Fig. 6. 

(k) At 3000Mc./s. and using a +1:7° beam, Piddington and Minnett 
made a similar measurement to that described above, and found a difference in 
equivalent temperature of 2-6 deg. K. The beam-width correction is 0-89 and the 
zero-level correction 1°19 (as before) so that 7, =2-77 deg. K. If the contours 
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at these highest frequencies are more closely spaced than at 200 Mc.,s., then both 
estimates of 7, should be reduced, the value for 3000 Mc./s. perhaps more than 
that for 1200 Mc./s. Judging from the curves of Fig. 6 the correction should be 
small and should not atfect the discussion of Fig. 4. 


Division of Radiophysics, 
CSE, 
Chippendale, N.S.W., 
Australia: 
1951 January 15. 
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THE OCCULTATION OF ANTARES OF 1950 JUNE 27-28 


David S. Evans* 
(Received 1950 September 11) 


Summary 


The variation of light intensity of Antares was measured during the 
period of extinction at the occultation of 1950 June 27-28. The results 
cannot be fitted within the estimated limits of error by the assumption of a 
single symmetrical star with any law of darkening to the limb. The results 
are explicable by the hypothesis that Antares A is itself a double star, the 
components being a star of diameter 0”-027 and a companion 3 magnitudes 
fainter of diameter o”’:016. This solution, although mathematically 
satisfactory, is physically improbable, since the temperature of the companion 
would have to be very low. It is suggested as a preferable alternative that 
Antares A is pulsating and has an ellipsoidal form. 





Introduction.—\t was pointed out more than forty years ago by MacMahon (1) 
that it should be possible to estimate the angular diameters of stars by timing the 
period of extinction of their light when they were occulted by the Moon. Eddington 
(2) pointed out that difficulties would be introduced by the diffraction pattern 
produced by the light from the star passing the edge of the Moon. He estimated 
that, as a result of this effect, all stars should show a spurious angular diameter of 
0”-008. However, just before the war, J. D. Williams (3) computed the diffraction 
patterns in detail for stars of various assumed angular diameters, and showed that 
it should be possible to compute the true angular diameters of stars from the 
observed degree of smudging of the diffraction bands. At the same time, 
Whitford (4) reported observations of 8 Capricorni and v Aquarii. In these cases 
the diffraction bands expected for stars of zero angular diameter were actually 
observed. 

At intervals of about twenty years, the Moon comes sufficiently far south to 
occult Antares, and such a series of occultations takes place in 1950-1952. In 
view of the large estimated angular diameter of Antares (about 0”-040) (§) it was 
obviously desirable that attempts should be made to observe the phenomena with 
the Radcliffe reflector, and the organization of the work was entrusted to the 
writer. 

An attempt to observe the reappearance of + Sagittarii on 1950 April 9 was 
clouded out, as was the reappearance of Antares on 1950 May 4. 

At the Cape, on May 4, the weather cleared just at the time of disappearance 
and Guelke and Cousins (6) obtained an observation somewhat reduced in value 
by the presence of a 50-cycle ripple on their tracing. 

For the observation of the occultation of Antares in the early morning of June 28 
(Greenwich date, June 27) we were able to improve our arrangements. In this 
case Antares disappeared at the dark limb of the Moon; asa result, not only was the 

* Report of a project undertaken by David S. Evans and A. D. Thackeray (Radcliffe Observatory, 


Pretoria), J. C. R. Heydenrych, J. D. N. van Wyk, R. T. Jamieson and J. H. N. Loubser (National 
Physical Laboratory, S.A. Council for Scientific and Industrial Research, Pretoria). 
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period during which guiding would be necessary much reduced, but also we were 
no longer confined in our choice of guiding stars to those immediately north and 
south of the track swept by the Moon. We were therefore able to select a sixth 
magnitude star which itself was occulted later. ‘This star remained visible in the 
guiding eyepiece of the double slide plateholder throughout the necessary period. 
and a successful observation was obtained in a sky clear except for a little scattered 
high cirrus cloud. 

The equipment.—The equipment was provided by the Electronics Division of 
the National Physical Laboratory, C.S.1.R., Pretoria, working under the direction 
of Dr J.C. R. Heydenrych. It is impossible to speak too highly of the enthusiasm 
and skill which these workers brought to the task, and the success of the work is 
due entirely to their efforts. ‘The installation of the equipment and its operation 
during the rehearsals was undertaken by Dr Heydenrych and Mr van Wyk. 
Unfortunately, Dr Heydenrych then left on a visit to the U.S.A., and, at the 
successful observation, the work was done by Messrs van Wyk and Jamieson and 
Dr Loubser. 
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3. 1.— Electrical circuits. 


A 931A photomultiplier tube was used. ‘This was fitted in a housing adapted to 
one of the eyepiece holders which fitted onto the Newtonian double slide plate- 
holder. A diaphragm about 0-8 mm. in diameter was placed in front of this. The 
diaphragm was in the focal plane, and the cone of rays then expanded to a diameter 
of 3-4mm. before impinging on the sensitive surface of the cell. ‘The maximum 
response of the cell was at 4000+ 5004. Operated with 100 volts per dynode 
stage, and with 100 volts between dynode No. 9 and anode, the cell gave a maximum 
anode dark current of 0-25 microamperes. Operating characteristics were: 
Average sensitivity : 9300 microamps. per microwatt, or 10 amps. per lumen; 
Current amplification : 10°; Luminous detectivity : 10 !! lumens. 

The differential balanced amplifier consisted of a double triode 6 SN 7 (Fig. 1). 
The signal from the photomultiplier tube was fed directly onto the grid of one 

5 
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section. Since about go volts positive anode supply is required for the final 
stage of the photomultiplier tube, the triode had to be biased to — go volts. ‘This 
was done by means of a battery in the cathode circuit. A small potentiometer of 
1000 ohms resistance was provided for cathode bias to permit all irregularities in 
the tube sections to be balanced out so that slight changes in plate supply did not 
affect output signal. 

The grid of the other triode was connected to earth through a fairly high 
resistance, variable through coarse and fine controls. ‘This enabled the changes 
in internal impedance of the photomultiplier tube, due to background light, to be 
compensated over a very large range of background light intensities. ‘The 
100 kilohm potentiometer in the grid circuit of the first triode was used as a 
voltage divider, so that a potential of + go volts could always be applied to the final 
stage of the photomultiplier tube. 

‘Two microammeters were connected in the grid circuits to ensure that, when a 
balance had been obtained, the grids were not positive and drawing grid current. 
Had this occurred there would have been a reduction of sensitivity and stability. 
The most sensitive setting was secured when the grids were only slightly negative. 
The signal output was taken between the plates of the triodes and fed, via 
a sensitive potentiometer, to the Hathaway recorder. 

For the dynode stages of the photomultiplier tube a regulated power supply of 
1000 volts was used. ‘This was stable to within 100 millivolts over long periods, 
and to within 10 millivolts over short periods. ‘lhe positive side was earthed, so 
that the cathode of the photomultiplier tube was 1000 volts below earth. ‘The 
differential amplifier was fed from a regulated power supply operated at 280 volts. 
The filaments were fed from a 6 volt storage battery. 

The screened cables from the photomultiplier tube at the Newtonian focus of 
the telescope to the power pack and recorder on the observing floor were about 
60 feet long. ‘The computed time constant of these lines is 1-65 milliseconds. In 
addition a line was run up the telescope to a small pea lamp in the photomultiplier 
tube housing. ‘This was turned on for a time before using the equipment in order 
to dry out the cell thoroughly. This had been found to be necessary in the early 
stages owing to trouble with humidity. An intercommunication system was also 
installed, using Sound-Power handsets boosted with a small speech amplifier on 
the Newtonian observing carriage. ‘This proved extremely useful owing to the 
difficult acoustics of the turret. Drs Evans and Loubser on the observing carriage 
were able to keep in easy conversational communication with Dr Thackeray at the 
finder and clocks, and Messrs van Wyk and Jamieson at the recording apparatus. 

The computed response of the latter, on the ranges used, was completely 
linear. The final recording was on a Hathaway Recorder, registering on bromide 
paper 10 inches wide, driven at 10 inches per second. The response on Antares 
(with a reserve sensitivity factor of something over 100 times) was about 50 mm. 
The galvanometer used had an undamped frequency of 580 cycles second and 
a guaranteed accuracy of 2 per cent at a frequency of 200 cycles/second. ‘The 
computed sensitivity for the optical path-length employed was 209 mm. /milliamp. 
‘Transverse time markers were put on by means of an internal standard, while, 
in addition, a second galvanometer recorded the ordinary 50-cycle mains 
frequency. A third galvanometer ran idle, and provided a useful straight trace 
which served as a datum for measurement. 
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The observation.—The circumstances of the occultation for Pretoria were 
kindly supplied by Miss F. McBain of H.M. Nautical Almanac Office. These 
were: 

1950 June 27, D_= 22" 27™ 068 G.M.T., P =87°. 

R23" 41m 278 G.M.T., P = 285°. 

Age of Moon: 12-3 days. 
‘hese circumstances mean that the angles of the normal to the Moon’s limb to 
its direction of motion at the points of disappearance and reappearance were 
small. In what follows no correction for this inclination has been made, nor 
for the extension of the time scale which would result if the occultation had 
taken place behind a steep mountain slope on the Moon. 

The equipment was tested and tuned up well in advance of the occultation, 
and the 6-2 magnitude star H.D. 148760 picked up for guiding. Balances were 
obtained and centring on the diaphragm checked by moving the plateholder 
back and forth with the screw settings in the two coordinates. ‘The motor driving 
the bromide paper was started 25 seconds before the predicted time, and stopped 
10 seconds after disappearance had been called by Dr Thackeray. The observed 
disappearance was 0-8 seconds early, i.e. at 22" 27™ 058-2 G.M.T. 

The thirty foot length of tracing was then developed and found to include a 
cut-off of the expected type lasting a total time of 100 milliseconds, which would, 
at the adopted rate for the motion of the Moon (0"-437 per second), correspond 
to a diameter of 0”-044 for Antares A. However, a more detailed scrutiny 
soon showed a number of complications. 

Analysis of the results.—We had expected to see on the trace two steps, one 
corresponding to Antares B (distant 3” in P.A. 274°) about 7 seconds before the 
predicted time of occultation: then, at the moment of occultation we expected 
to see an oscillatory pattern corresponding to the diffraction bands, followed 
by a gradual cut-off. In fact we found the following: the traces were dead 
straight with only minute ripples (of amplitude less than 0-3 mm.) except at two 
points. About 20 seconds before occultation there was a small sharp kick, 
occupying about 5 milliseconds, causing a zero shift of about 15 mm. Interpreted 
as due to variation of incoming light this would have indicated an increase. After 
the zero shift the trace remained perfectly straight. A visiting astronomer, 
Dr John B. Irwin of the Goethe Link Observatory, University of Indiana, who 
was watching the sky at the time, reported that it was perfectly clear. In any case, 
the change was far too rapid to be due to cloud. The equipment appears to be 
subject to very occasional random zero shifts of this kind. Nothing else of the 
kind was detected on the whole of the rest of the tracing, but, earlier in the evening, 
a slight jump had been reported when the electric motor tracking the observing 
carriage had been operated. The probability that anything of the kind occurred 
during the critical 100 milliseconds is very small, and would in any case almost 
certainly have shown up as a vertical break in the curve. Care was taken 
not to operate any of the turret motors during the recording period. Apart from 
this, random zero drift appears to be less than I mm. over the whole period of 
operation (35 seconds). 

There was nothing on the tracing which could be definitely identified as the 
occultation of Antares B, nor were there any oscillations indicating the presence 
of diffraction bands. 
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The curve presumed to be due to Antares A also showed a number of 
unexpected features (Fig. 2, Plate 2). The drift rates before and after the 
occultation showed a slight difference—about 0-055 mm./millisecond before and 
0-010 mm./millisecond afterwards. At the top, at the beginning of occultation, 
the curve showed a smooth rapid turnover, but the lower toe tailed off slowly, 
making the whole curve asymmetrical. We considered a number of possible 
explanations for this. ‘The time constant computed for the circuits is rather high 
(1°65 milliseconds), but this is nowhere near adequate to account for the final 
slow change as due to instrument lag. We were also able to check on the value 
of the time constant by examining a number of traces in which, with the paper 
running at slow speed, the roller blind shutter at the Newtonian focus had been 
snapped on and off through the (out of focus) star image. ‘The shutter is rather 
slow, but it was possible to assign a maximum possible time constant of 
5 milliseconds, much of which must be due to the lack of focus and to the shutter 
action. What would be needed to account for the observed effect as lag in 
recording would be an averaging over at least 15 milliseconds. Ifthe time constant 
were as large as this the initial turnover and the rate of change on the straight 
line portion could hardly be as rapid as they are observed to be. 

We then considered the possibility that the star might not have been fully on 
the diaphragm. The hole was about 0-8 mm. in diameter and the separation of 
Antares A from Antares B about 0-14mm. This possibility was also rejected, 
since, during practice, a number of tracings had been taken with the stars just 
at the edge of the diaphragm. When this is the case there are large rapid 
fluctuations due to seeing. The straightness of the traces seems to rule this out. 

Computations were then made, combining the spectral sensitivity curve of 
the photomultiplier tube, the colours and magnitudes of the stars, atmospheric 
transmission, and the reflectivity of silver to determine (a) the expected response 
on Antares B (4 magnitudes fainter than A) and (4) the effective colour response 
for Antares A. It was concluded that the response on Antares B should be in the 
region of 1 per cent or less of the response for Antares A. ‘The effective wave-length 
for Antares A was 5100 A. and the response curve is narrow, the 50 per cent 
values being at 4300 A. and 5800 A. 

The first result cleared up the question of the absence of any definite response 
from Antares B, and the analysis of the trace was embarked upon, taking the 
observed curve to be a true representation of the variation for Antares A. 

Using the idle galvanometer trace as a datum line the deflections were 
measured and plotted against the time. Here a difficulty was encountered, since 
it soon became clear that the internal time marker of the Hathaway recorder had 
an error of several per cent. This is probably to be attributed to the fact that 
all the apparatus had had to be transported by car a distance of several miles over 
gravel roads in very poor condition, and in spite of the greatest care it had been 
badly jolted several times. ‘The time basis adopted was therefore the 50-cycle 
trace from the mains supply. 

The Pretoria power station authorities reported that, at the relevant time, this 
station had been running in parallel with the Electricity Supply Commission 
generating station at Germiston. ‘The Chief Control Engineer there was good 
enough to report that, at the relevant time, the metered frequency was 
50°05 +0°05 cycles per second. Fortunately, too, there is an arrangement 
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Fic. 2.—Observed trace (slightly enlarged). 
Time increases to right, intensity upwards: 50-cycle mains trace at bottom. 
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between the Union Observatory, Johannesburg, and the Electricity Supply 
Commission whereby a differential clock reading between the mains frequency 
and the quartz-crystal clock of the Union Observatory is exhibited on the 
generating station control panel. The reading on this clock at midnight was 
2 seconds and was unchanged at I a.m. South African time. It seems clear then 
that if we adopt a frequency of 50-05 cycles per second we shall have a time standard 
more than adequate for the purpose of reduction. 

The rate of travel of the bromide paper was slightly variable about an average 
of 10°5inches per second. It was possible to identify the times corresponding 
to the transverse marks in milliseconds from an arbitrary datum. Finally a curve 
was derived showing the pre-occultation readings as 100 units and the 
post-occultation as zero, with the abscissae in milliseconds. 

As has been remarked by Eddington, Williams and others, the analysis of such 
a curve is beset with diffraction difficulties. Asa preliminary a series of computa- 
tions was made for uniform disks of various sizes using the data and formulae 
supplied by Williams, which are particularly convenient for calculation. But 
Williams adopted a mean wave-length of 43004. and a mean lunar distance of 
3-8 « 10" metres, whereas in the actual case the mean wave-length is 5100 A. 
and the lunar distance 3-6 10° metres. Williams works in terms of v, the 
argument of tabulation of the Fresnel integrals, and his figures are so chosen that a 
change of 0-1 in v would correspond to a shift in position of a point source through 
an angle of 0”-0005, whereas in our case the shift would be 0”-00055, i.e. 10 percent 
greater. However, in the event it was found that the rate of change of the 
function {I (diffraction) —I (geometrical) } is very small over the range of diameters 
from about 0”:024 to 0-040 in which we are interested. ‘Throughout this range 
the effect is to raise the intensity computed from geometrical optics by about 
I°5 per cent near extinction, to lower it by about 1°5 per cent throughout the 
duration of eclipse and to raise it by between I per cent and 2-5 per cent (the 
larger figure corresponding to the smaller diameters) in the region just before the 
beginning of eclipse. ‘This result cleared up the difficulty that no oscillations of 
brightness were observed just before the occultation. All the fringes have been 
almost completely smudged out, but the effect of diffraction is not negligible. ‘The 
geometrical curve is modified by its effects, and an error will be committed if it is 
neglected. In particular it became clear that merely to measure the time interval 
between the points where the darkening curve deviates from the best straight lines 
which can be fitted to the traces just before and just after the occultation will 
result in a value for the diameter which is more than 10 per cent too great (Fig. 3). 
Moreover, diffraction will produce what looks like a differential rate of zero drift 
just before and just after the occultation, since, just before the occultation, the 
brightness will be rising slowly. We attribute the small differential drift already 
noted to this cause. ‘The two drift rates should not be quite the same, and to draw 
lines which make them equal will lead to error. 

The situation which now presented itself was somewhat confusing. Clearly 
the most accurate fit would be obtained by using the slope of the extinction curve, 
but this slope is affected not only by the assumed diameter of the star but also by 
the law of darkening over the disk, and these two together will determine the 
diffraction which will, in its turn, affect the slope of the observed curve. Further, 
the effect of diffraction will modify the curve from symmetry about its centre in the 
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sense actually observed. At the top it will turn down fairly sharply, while at the 
bottom, where the intensity is low, it will drag out the curve into the type of “ toe”’ 
shape actually observed. 

The procedure adopted was then the following: assuming a disk of a certain 
diameter the diffraction was computed by Williams’ formulae. As has been 
pointed out, the diffraction deviations are not at all closely dependent on the dia- 
meter assumed so that, within limits, one can apply results computed for one 
diameter to another, and one need not bother with changing to the true values of 
lunar distance and effective wave-length. Further, since the wave-length-response 
curve for Antares A is so narrow, it will be sufficient to make the computations for a 
single wave-length. This, of course, would not be the case for very much 
smaller diameters where the diffraction produces rapid changes of intensity, and 
where the curve alters appreciably with changing diameter. 


DURATION OF 
EOMETRICAL  —f 
OCCULTATION “a 








50 
Fic. 3. 
Observations. 
Computations, darkened disk, 0”-027, including diffraction. 
Computations, uniform disk, 0-024, including diffraction. 

These simplifications are fortunate since the computations are very laborious. 
By using Williams’ values one can read off the values of the Fresnel integrals 
from tables in forming the table of J (the intensity distribution for a point source), 
whereas, if fractional values of the argument had been involved, it would have been 
necessary in deriving each of the 150 entries needed in the table of J, to interpolate 
or to compute each value of these integrals afresh. As it is, in the computation of 
each point on the intensity curve for a disk of diameter 0”-040 the calculation of an 
integral expressed as the sum of 80 products is involved. 

When the diffraction had been computed, the slope of the curve was fitted to 
the observed slope to derive a value for the diameter of the star. If this agreed 
with the diameter originally adopted, the solution so found was compared with the 
observations (Fig. 3). It will be seen that, while a reasonable fit can be obtained 
for the upper part of the curve, there is a large discrepancy at the toe, which, we 
believe, is far outside the possible range of instrumental error. 

In the hope that diffraction alone might be capable of accounting for the 
observed discrepancies, some consideration was given to other models. For 
example, the possibility that Antares A might consist of a bright central star 
surrounded by a fainter atmosphere was considered and rejected. In the end, 
the following appreciation of the situation was reached: the total period during 
which the brightness is observed to be changing is too long for the slope of the 
straight line part of the curve. Any symmetrical distribution of brightness 
variation over the stellar disk must produce a geometrical extinction curve having 
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symmetry of the following type: namely, that if the moment of mid-eclipse is taken 
as the zero of time (when the intensity is 50 per cent), and if two moments are 
selected at times +¢ and —t at which the intensities are J and J’, then J+J must 
be constant. ‘That is, if the geometrical curve were rotated through 180° about 
its centre, it would be unchanged. The effect of diffraction will be to raise the 
intensity at the beginning and end of occultation (i.e. where the geometrical 
values are 100 per cent and zero respectively) by something of the order of 1-2 
per cent, and to lower the intensity at mid-occultation by about the same amount. 
We have not been able to devise any model with central symmetry in which the 
diffraction will raise the intensity at the end of occultation from the geometrical 
value of zero to the observed value of 6 per cent. 

It is also worth noting that this argument does not depend on circular symmetry 
in the star. Central symmetry, e.g. an elliptical disk having symmetrical 
illumination, will produce a geometrical curve having the same symmetry property. 

At this impasse a search of the literature was made, from which it appears that, 
on spectroscopic grounds, duplicity or deformation by pulsations of Antares A 
has been suggested by Spencer Jones (7). But for this suggestion we should have 
hesitated before proposing such a solution to the problem, since the introduction 
of a second star, or distortion of a single star, provides so many disposable variables 
that almost any curve might be fitted. 

Assuming first that Antares A was double, we proceeded on the following 
lines. First, the best geometrical fit was obtained in the two cases (a) uniform 
surface brightness, (6) the extreme case of darkening to the limb represented by 
1=I, cos #. ‘The computed results were subtracted from the observations and a 
curve resembling that for the occultation of a second star obtained. ‘This was 
then analysed and the computed values subtracted from the observations. ‘These 
values were then reanalysed, taking diffraction into account, to give definitive 
results for the larger star. The new curve of residuals was then analysed again 
to give a definitive result for the smaller star. 

The best results are: 


Larger star: Diameter 0"-027, Darkening / = J, cos 6, Luminosity 0°939. 
Smaller star: Diameter 0”-016, Uniform disk, Luminosity 0-061. 
Separation between centres projected on to line of Moon’s motion: 0”-0182. 


‘The accuracy of fit is to about 1 percent. For a uniform disk the diameter of 
the larger star would be 0”-024 but this solution is less satisfactory than that given 
above for the following reason. For a uniform disk of this diameter the diffraction 
produces a rise in intensity just before occultation of 2-4 per cent (Fig. 3). This 
has not been observed, and we beiieve that it would be observable. ; 

It will be noticed that the projected separation between the centres of the two 
stars is slightly less than the sum of the radii. This would not necessarily mean 
that the stars partly eclipsed each other, since all that is determined is the projec- 
tion of the separation on the direction of the Moon’s motion. ‘The second star 
may be centred at any point on a line perpendicular to this direction. In any case, 
even if there were an eclipse, the intensity reduction would probably be so small as 
to be unobservable. 

Discussion.—Now let us consider this solution critically. The curious result 
obtained makes it a matter for great regret that a successful observation was not 
secured at Pretoria on May 3-4, and that the Cape result is such that it is impossible 
to determine the shape of the curve precisely. There is no opportunity for a 
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further observation during the present year, and in 1951 only one occultation is 
visible at Pretoria, and that a reappearance, which we cannot be entirely confident 
of observing successfully. 

What we have been after is a description of the shape and intensity distribution 
of a patch of light. We have found it impossible to make that patch circular in 
shape, the shape to which general ideas lead us first. But we can get concordance 
if we assume that there are two circular patches. If these patches are stars, how 
do they accord with our knowledge of stars gained from other sources? The 
projected separation is small—assuming a parallax of 0”-0095 it is about 2 astro- 
nomical units. It is inherently suspicious that the light from one star has to take 
over immediately that the other has finished its contribution to the light variation. 
Physically there is a discontinuity: observationally there is none, and we have to 
assume that the two stars conspire to produce a smooth curve with no plateau. 
We have no evidence one way or the other, but the result looks odd. Spencer 
Jones in his analysis found a value of a sin 1=67-4 x 10° km. =0-45 astronomical 
units, and was most careful to insist that the analysis in terms of a spectro- 
scopic binary model did not necessarily imply duplicity, but might indicate that 
Antares A was pulsating. Others have been equally dubious as to the inter- 
pretation in terms of a true binary (8). 

According to our solution the larger star would have 2-85 times the surface 
of the smaller, and would emit 15-7 times as much light. Hence each square 
centimetre of the surface of the smaller star would emit 0-181 of |» radiation per 
unit area of the larger. ‘This is, of course, in terms of light to which the photocell 
responds. If these were bolometric data the conclusion would be that the tempera- 
ture of the smaller star was about 1400 deg. K., which seems impossibly low. 
The result is not very different if the calculation is worked through in terms of 
radiation of 5100 A. emitted by black bodies. 

We might defend our binary solution on the ground that interferometer 
measures of such a close double would be approximately those which have 
actually been found: we might argue further that the diameters computed from 
radiation laws involve the parallax, and that quoted parallaxes for Antares show 
a range of three to one. But we cannot evade the conclusion that the assumption 
of two separate stars leads to a very unsatisfactory temperature for the second 
component. If we assume that Antares is pulsating and is a single star of distorted 
form, we find the difficulty that we are without any kind of model to go on. 

Suppose that Antares has the shape of a prolate spheroid with the major axis 
at some angle to the direction of view of the observer (Fig. 4). As observed, 
this spheroid will have a variable surface brightness which will be greatest near 
the region where the line of sight is normal to the surface (point P). ‘The iso- 
photes round this point will be roughly circular, and the distribution of intensity 
will probably approximate to that of a spherical star. The more distant, “ tail 
end’’, of the spheroid will provide the light from the “ second component ” 
of the ‘‘binary”’, and this will be of low intensity, since here the line of sight 
makes a large angle with the normal to the surface. At points above and below 
the plane of the paper (Fig. 4 (4)) this angle will be greater still. The situation 
will be much as illustrated in the diagram. The general picture is supported 
by the fact that for the “ brighter component’”’ we found the best accordance 
on the hypothesis of complete limb darkening. 

In order to make the picture more precise we must draw an ellipse enclosing 
the “ binary ’’ already deduced. Ifthe projection of the major axis of the spheroid 
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Fic. 4. 

(a) Binary star solution. 

(b) Ellipsoid (plan). 

(c) Ellipsoid (elevation) showing approximate isophotes. 
on the sky is parallel to the projection of the motion of the Moon, we find a projected 
length for the major axis of 0”-040 and a projected length for the minor axis of 0”-027. 

It should be noted that this solution should be in accordance both with the 

interferometer measures and the radiation computations. Further, it seems likely 
that some difficulty may be encountered in determining the parallax if the centre 
of brightness shifts on the star’s surface in the course of pulsations. A shift of 
only 10 per cent of the estimated major axis is almost 50 per cent of the measured 
parallax. 
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In carrying out this project it was hoped that we should solve a problem, 
instead of which we have merely created a number of new ones. ‘There seems no 
doubt that occultation observations should be repeated at every opportunity, 
since the expectation is that, if Antares changes shape by pulsation, the results 
will differ appreciably among themselves. It would also seem desirable to try to 
make interferometer observations in different position angles. ‘There is obviously, 
too, a considerable scope for detailed theoretical discussion of these results. 

In the course of this analysis I have been much helped by discussions with 
Dr John B. Irwin of the University of Indiana, a specialist on eclipsing binaries, 
who has been visiting this observatory. I also wish to record my thanks for the 
very illuminating and profitable discussions which I have had with Dr A. D. 
Thackeray. Warm acknowledgment is also due to the staff of the electronics 
division of the S.A. National Physical Laboratory for their expert work, and to 
the Director of the Laboratory, Professor S. M. Naudé. Drs Heydenrych 
and Loubser and Messrs van Wyk and Jamieson wish to acknowledge the 
permission of the South African Council for Scientific and Industrial Research 
to participate in this investigation. 

I am indebted to Dr C. B. Watts of the U.S. Naval Observatory, Washington, 
for a trace of the relevant portion of the Moon’s limb. It is unusually smooth, 
and the region of the occultation is a smooth flat plateau. Naturally, however, 
detail shown on this trace is larger than that which would be relevant in the present 
case, since, as was pointed out in Monthly Notes of the Astronomical Society of 
South Africa, 9, No. 8, 1950 August, the projected diameter of Antares at the 
distance of the Moon corresponds to a distance of only about 75 metres. It will 
be appreciated that the effect of an irregularity in the Moon’s limb will be greatest 
when that irregularity is of almost the same dimensions as the projected diameter 
of the star, and when it passes centrally across the disk. ‘The following cases had 
been computed for a uniform disk : (1) straight edge; (2) straight edge with slope 
joining it an at angle of 45°, the intersection passing centrally over the disk; 
(3) central passage of a hemispherical depression having half the radius of the star : 
(4) same as (3) but depression has same radius as star. ‘The results are: the general 
slope of the curve is unchanged in all cases. Cases (2) and (3) increase the total 
duration by 25 percent. ‘The degree of asymmetry introduced can be measured 
by the displacement of the point of 50 per cent intensity from the mid-point of 
duration, measured in terms of the duration. The values of this displacement are : 
Case (2) I per cent, Case (3) 2 per cent, Case (4) 10 percent. ‘The required value 
is more than 20 per cent. In addition, Case (4) does not resemble the observed 
curve in other respects. 


Radcliffe Observatory, 
Pretoria : 
1950 August 18. 
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ERROR BALANCING IN FAST SCHMID'T CAMERAS 


E. H. Linfoot 
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Summary 


An approximate solution is obtained to the problem of designing the 
Schmidt camera so as to obtain best performance, in an agreed sense, over 
the field surface taken as a whole. The approximation consists in neglecting 
small quantities of the seventh and higher orders and basing the discussion 
on the fifth-order aberrations of the system. ‘The optimization covers 
arbitrary variations of the plate profile and of the adopted field surface ; 
it also takes account of the spectral characteristics of the light in which the 
camera is to work. In all cases, the best field surface turns out to be 
spherical and the best corrector plate profile such as would give axial 
stigmatism to fifth-order accuracy when the refractive index of the corrector 
plate is assigned a suitable fictitious value. In a system intended to work in 
monochromatic light, the best position of the neutral zone is at the edge of the 
plate aperture and the root mean square image radius is then } of that of a 
classical Schmidt camera. When the light is not monochromatic, the best 
radius of the neutral zone always lies between H4/§# and H, where H is the 
radius of the aperture stop. Its value is given by a simple expression 
depending only on the angular diameter of the field and the effective 
index-spread of the corrector plate for the given spectral distribution. 

In the last section of the paper it is shown that a slightly better root 
mean square performance can be obtained by the use of the modified 
Schmidt system described in a previous paper (x), and the relation between 
the properties of the classical and modified Schmidt systems is analysed. 





1. /ntroduction.— The problem of the design of aspheric systems is essentially 
one in the calculus of variations. ‘lhe profiles of the different aspheric surfaces 
are described by arbitrary functions which have to be chosen so as to optimize, 
in some agreed sense, the performance of the system over a suitably chosen 
field surface. If this field surface has to be flat, its position is another free 
variable which the process of optimization must take in; if a spherically curved 
field surface is allowed, its description involves two free variables; while if 
the field surface may be of arbitrary form, the determination of its best profile 
requires that another arbitrary function shall be covered by the optimizing 
process, in addition to those which describe the profiles of the optical surfaces. 
In systems designed to work at a finite distance, variations of the obiect surface 
also come up for consideration. 

The analytical difficulties of the general problem place it wel! beyond the 
reach of present-day geometrical optics, but in one special case of practical 
importance I have been able to find an approximate solution. This is the case 
where the system effectively contains only two surfaces, one of them a reflecting 
sphere and the other a near-flat aspheric refracting surface, where the object is 
at infinity, and where the assessment of performance is based, in a manner more 
fully described in Section 4, on the notion of mean square image radius over the 
whole field. 

The solution, obtained in Sections 2—6 of the present paper, provides an 
approximate answer to the question: How can the errors of a fast Schmidt 
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camera best be balanced (in the sense of Section 4) over a given field and what 
departures from the classical design of plate profile and of field surface are needed 
to secure this balancing? In Section 7 it is shown that the modified Schmidt 
described in an earlier paper (1) can give a better mean square performance than 
the “ best’? Schmidt camera with spherical mirror, but that the improvement 
in mean square image radius is not large enough to be of much practical 
consequence. 

Some progress towards an analytical theory of error balancing in Schmidt 
systems was made in a recent paper (2) by E. Wolf and the writer, where a solution 
was obtained to the problem of optimizing, in the above sense, the errors of a 
Schmidt camera over a given angular field in the “‘customary”’ field surface 
(defined in Section 2) by modifying the profile of the corrector plate. This joint 
work formed the starting point of the present investigation. 

2. Notation; the aberrations of the ordinary Schmidt camera.—Fig. 1 shows 
the optical system of the ordinary Schmidt camera. For the present we suppose 
that the aspheric surface of the corrector plate faces the spherical mirror and 
passes through its centre of curvature C. x, y, z are rectangular Cartesian 
coordinates with C as origin and with CA as z-axis. We take the aperture stop 
of the system to be the circle x? + y? = H® in the plane z=0, i.e. in the tangent 
plane to the aspheric surface. In an actual system of focal ratio f/2°5 or longer, 
the aspheric surface lies everywhere within a few thousandths of an inch of this 
plane. 


Fic. 1.—The Schmidt camera. 


It is convenient to define the parameter 
w=H/R, (1) 


where R is the radius of curvature of the mirror, and to introduce new variables 
u, v in the plane of the aperture stop by means of the equations 


x=pRu, y=pRv; (2) 


then the boundary of the aperture stop is the circle u°+v*=1 and the equation 
of the aspheric surface can be written * 


(n—1)3 = pRpA(r* — ror’) + RyA(3r8 — 374 — Fr*) + O( Rp’), (3) 


16 
where n is the refractive index of the plate in light for which the system is axially 
stigmatic, r* =u? +? and ry is the radius in (u, v)-coordinates of the neutral zone 
on the plate. In a practical case, the value of rp is usually about 0-8. 


* See (1), equation (37). 
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If terms of order Ru® may be neglected, (3) can be written in the simpler form 
(n—1)z = } Ru4(r4 — ar*) + O(Rp*), (4) 


where, to a sufficient approximation, @ = 27;.t 

Rays entering the system parallel to the axis CA, and distant ry from it, 
meet the aspheric surface of the plate on the neutral zone and are therefore not 
deviated by the plate. After reflection by the mirror, they meet in a point I on CA, 
which we call the axial image point of the camera. ‘The sphere centred at C and 
passing through I is customarily taken as the field surface of the camera; its 
radius CI, which is also the focal length of the system, is given by the equation 


f=Cl=}RU-?n)'*. (5) 
It follows from (5) that 


1 — per)? 


I o » 1/2 
4F“ —p roy”, (0) 


where F is the focal ratio of the system. Ina system of focal ratio f/2°5 or larger, 
p< i003 in an f/1 system p?<4. Thus we are justified in treating pz? as small in 
approximations. 

If the system works over a field of angular diameter 24, in the above customary 
field surface, the linear diameter of the field surface is 2CI sin 4, and the central 
obstruction ratio 


c=2F sin dy = $0 {1 + O(u*))]. 
2h 


In Fig. 1, 1, represents a point at the edge of the field surface in the plane xCs ; 
I’ represents an intermediate field-point in the same plane, at an angular distance 4 
from the centre of the field. ‘The parameter V’, defined by the equation 


, sing 
sien sin dy’ (6) 


runs from zero at the centre of the field to 1 at its edge. 

The point I is, in third-order (Seidel) approximation, the image corresponding 
to an object point at infinity in the uz-plane, at an angular distance ¢ below the 
axis of the system. ‘That is to say, the rays of a parallel pencil from this object 
point will eventually meet the field surface in points whose distances from I’ 
are O( Rp’). 

More precisely, the angular displacement components from I’, measured in 
seconds of arc, of a ray which enters through the point u, wv of the aperture stop 
are respectively 

mee” ae — es 
bX,5Y=Kp° (55° x) @*(u, v, VV; n')+O( Kz‘), (9) 


where K =648,000/7, 


O*(u, v, V; 1c?V? (: + 5.) tans | (04 — 20%), (10) 


+ We can define a= 2r' if the power series expansion of the error term in (4) is allowed to contain 
a term in r? ; otherwise (3) and (4) give a=2r§+ O(p*). 
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n’ is the index of refraction of the plate for light of the wave-length under con- 
sideration, m that for the light in which the system is axially stigmatic and 
n'—n 
Res, 4c(n—1)” (11) 
‘The value of &,,,,,, seldom exceeds 2 in a practical case, and the inaccuracy 
represented by the error term in equation (9) usually amounts to only a few per 
cent of the mainterm. It follows that the “ fifth-order aberration components” 
oD* ao* , 
du’ Ov ) aes 
describe the geometrical images with an accuracy amply sufficient for the needs of 
the optical designer. 

When the plate, of thickness O( Ry), is turned with its aspheric side away from 
the mirror, all the above results continue to hold, provided that, in measuring the 
distance of the aspheric surface from the pole A of the mirror, that part of the 
distance which lies in glass is taken with a factor 1/n. 

In practice the system usually works, not in monochromatic light, but over an 
effective range which varies in different applications, and the value of n’ may vary 
by a few per cent over this range. Such a variation may be said to be O(”) in the 
case of a fast Schmidt camera. 

‘The formulae (9) remain valid if the definition (10) of ®* is replaced by the 
new definition 


bX*, 5Y*= Ku ( 


Sead ef 0 - a 
*(n, v, V;n')= [ ser 2 (= tus) 4 Ae | (r4 — 2ar-r’), (13) 


in which np is the refractive index of the corrector plate at a selected point of the 


effective spectral range, k\).,, is given by the equation 


n'—n 
ee: en 
.=— I 
a mn } 4p7(ny — 1) ( 4) 
and it is assumed that 
ny =n + O(u"*). (15) 
(13) is more convenient than (g) for discussions which include chromatism because, 
cs : ») 
in (13), anda enter only through the quantity k, 
The new expression 


| wn akitge [MED au 4 #8 4 in Bt 
5X*+16Y* =} Kp*¢* [ =" 5 Th * et Oe (r* — ar*) 


saxon (2 +i) (rt—ar’) (16) 
for the fifth-order angular deviations in the customary focal surface is effectively 
independent of mp, in the sense that it is altered only by an amount O( Ky’) 
if the choice of m, is varied in accordance with(15). It is permissible to use a in 
place of 275 in defining these aberrations, since the corresponding change in 
their values is only O( Ky‘). 

By inspection of (16) we see that, to the order of numerical approximation 
represented by the fifth-order theory, the colour-error of the system is a sphero- 
chromatic aberration, uniform all over the field and of the same “‘ shape” r4 — ar? as 
the corrector-plate profile. 
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‘The parameter p, defined by the equation 


/ 
n—n 


pd = (17) 


’ 
Ny—1 


provides a convenient means of measuring this aberration. By (7) and (14), 


nn’ 


on = he [1+ O(y2)]. (18) 


Evidently p increases as the wave-length shortens, and p =o corresponds to the 
wave-length for which the system is axially stigmatic. 

The formalism of the general discussion can be simplified by introducing 
normalized deviations t U*, V*, connected with 6X*, 5Y * by the equations 


* = a 7% "* 4 : "* ( 
l Kop DX *, | Kage? ; (19) 


7” I :, 
by writing « = — and by introducing the new parameter 
; ; ) 
0 


,_ ¢ _ sing : 7 
O= 5 = Fag E+ OW) (20) 


(© is for all practical purposes the same as the parameter V). ‘Then (16) becomes 


U*+i1V* =0? | + 2 +hu a ried + hu Al (r*—ar*) 


0 70 4 2 
+e(5. HiZ Jo —ar-) 


ou" ov OQudv 


(U*, V*) =2(u, v)[ p+ (1 + «)0*|(2r? — a) + (u, v)ad® 
+(4u®—au, — 4v° + av)O?. 
The effect of replacing p by p — py in this equation is to add a term 
—}Kp* pods (5. 7 ix) (r4 — ar”) 
to the displacement components (16) and a term 
— pot*(r* — ar’) (23) 


to the aberration function (10). This represents a pushing back of the wave- 
fronts converging towards the images in all parts of the field by an amount 

Spgc*p8(r' — ar*) + O( Ry’). (24) 
When fp, >0, it means that a smal! amount of spherical under-correction has been 
introduced into the system. { 

On the other hand, the effect of replacing p by p— fp in (21) is to move the 
“stigmatic wave-length”’ of the system from that corresponding to p=o to that 
corresponding to p=py. It follows that, in the approximation represented by the 

+ This improvement in notation was suggested to me by Dr E. Wolf. 

t Although the effect of the term in ar* is to change the position of the best focal surface, it 
seems better not to separate off this term as focus shift, because the case a=o does not correspond 
to a zero shift of the best focal surface as py varies. Instead of doing this, we use the notion 


“spherical aberration of type r*—ar®’’, which corresponds to the “ doughnut” of the practical 
mirror-maker. 
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fifth-order theory, the introduction of primary spherical under-correction into 
the system by change of plate profile is equivalent to a shift of the ‘‘ stigmatic 
wave-length’’ towards the blue end of the spectrum. Of course, this would no 
longer be the case if the seventh-order terms were included in the analysis. 

3. Variation of field surface—In the field surface so far adopted, the 
fifth-order angular deviation components, measured in seconds of arc, are 


om* oao* 
ou’ dv 


where @* =@*(u,v, V; n’) is given by (13). A moving forward (against the light) 
of the focal surface by }efutc? adds Kep*c?(u,v) to these angular deviations; it 
can therefore be expressed by adding }ec*(u? + v*) to the aberration function ®*.f 
To the normalized displacements (U*, V*) it adds e . (u,v). 

More generally, a moving forward of the axial point of the adopted focal 
surface by fc*u‘c,, combined with a lengthening of its radius of curvature by 
| fu?d,, results in a moving forward by 


Spice + dO") + O(Rp*) (25) 
of the part of the surface at an angular distance ¢ from the axis. This adds 
2(c, + d,O*) (u, v) (26) 


to the normalized displacements (U*, V*). 

We call c,, d, the parameters of the new field surface with respect to the 
original field surface as reference surface. An arbitrary spherical field surface, 
symmetrical about the axis of the system and lying everywhere within O( fu‘) 
of the original field surface, is specified by giving c,, d, suitable “‘finite’’ values. 

4. Performance over a given field; best field surface in a given wave-length.— 
The notion of optimizing the performance of an optical system over a given field 
can be interpreted in various ways. ‘The most primitive interpretation, and for 
some purposes the best, is to take the greatest image diameter for all object points 
in the given field and for all wave-lengths in the adopted spectral range, and to 
say that the performance is optimized when this diameter is made as small as 
possible. But even if diffraction effects are ignored, this procedure does not 
correspond very well to the requirements of the physical situation when, for 
example, the effective intensity of the light rises to a peak at a certain point in the 
spectrum and tails off to very small values at either end. 

Moreover, in cases where the size of the image is comparable with the resolving 
power of the photographic film, the performance is affected by photographic 
image-spread and there is no very close correspondence, even in monochromatic 
light, between image diameter as defined above and the diameter of the 
star images actually formed on a photographic plate. 

In such cases it seems more appropriate physically, as well as more convenient 
analytically, to define the effective radius of a single monochromatic image as the 
square root of the expression 


(8X*, 5Y*)=Kp? ( 


*|{ (8X)? + (5Y)*] dude, (27) 


ul+o<l 





+ This corresponds to focusing by trial. Deviation terms of order Ky’ are of course neglected 
in the present approximation. 
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and the effective monochromatic image-radius for the system working over the 
region @< ¢y of a given field surface as the square root of 
. me tx e\ r\9 r\9 
E=—, | db dd | | (8X)? + (SY)*] du dv. (28) 
m7$5 J 0 7/ui+u'sl 
‘This definition is therefore adopted in the present discussion. It is equivalent 
to defining the effective radius of each stellar image as the radius of gyration of 
its ray-density distribution about its centre (6X, 5Y)=(0,0), and the effective 
monochromatic image-radius over the whole field as the mean square average of 
these effective radii over the field-area ¢ ~ dy. 
When the expressions (27) and (25) are respectively replaced by 
- | | [(8X*)? + (8Y*)?]du dv 
ur+v'<1 


7 


I _ 
= — Kp%4! U®+4V%) dudo 2 
107 pi _—, r ) ” ( )) 


E* = — r b dd | | (OX *) + (6¥*)?] du dv 
( 1 


) 
75, 2. 


) ~/ui+v 


= Kuss | Odo (U*+V) dudv (30) 
Sr J0 J/7 @ 


+v®— 1 
l £72,644, 
= 7, R*e poe* T, (31) 
the changes in value are negligible in the present approximation, since they can 
be expressed by an error-factor 1 + O(*). 
The square root of the quantity 


I 
e* = | 
T 


20 dO || (U"+V*) dudv (32) 


v 1 


~0 


will be called the normalized image-radius of the system over the given field. 
The etfective image-radius in seconds of arc, viz. 


VE*=1K3¢2Ve*, (33) 


is different for different wave-lengths, and in assessing the performance of the 
system it is usually necessary to take account of the fact that the images are not 
monochromatic. 

To determine the ‘“‘best”’ field surface for the system working in 
monochromatic light of given wave-length, we impose on the customary field 
surface an arbitrary symmetrical deformation, which moves forward the field 
surface by an amount $fc*‘e(@) at an angular distance ¢ =©¢, from the centre 
of the field. ‘he even function e(@) may be taken to be everywhere O (1) without 
real loss of generality, and it is then immaterial in the fifth-order approximation 
whether the field-displacements are taken to be normal to the original field surface 
or parallel to the axis. 

In the new field surface, the normalized deviations of the system are 


U*, VS = Ut +ue(O), Vi +ve(0) 





+ This e* is twice the e* of (2). 
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and the value of e* is given by the equation 
1 * 
e* = | 20 dO | | {[US +ue(O)}* +[V* + ve(O)]}?} du du 
0 ui+v'<1 
1 - 
=e) + =| 2e(O).20d0 | | (uU* +0V#) dude 
0 Jd 


1 an 
Py. | €(0)20 do | | (u? + 0%) dude, 
Tig ’ 


where Uf, V3, ef denote the values of U*, V*, e* in the customary field surface. 
Thus 
“1 1 
e* =e* + | (0) #(0).20d0+ | (0)0d0, (35) 
0 Jo 
where 
H(O) = : | | (UU? +vV*) dudv. (36) 
~/v+e<l 
An infinitesimal variation 5e(©) in the function e(@) changes e* by an amount 
“1 “1 Pa 
| 5e(0)#(0)20 dO + | 2€(O)5e(O)O dO = | de(O).[ (0) + €(O)].20d0. 
0 0 Jo 
When the surface is a best field surface, this must vanish for every choice of 
5e(Q). Therefore in the best field surface 
«(Q) = —- #(0) 


; | | (uU} +vVF) du dv, (37) 


J ut+or< 


by (36), 


=2{(a—$)[p+(«+ 4)0"] —30%}, (38) 
on substituting for U*, V* from (2 


that the best field surface is, in fifth-order approximation, a spherical surface 
of parameters 


2) and evaluating the integral. It follows 


¢,=(a—$)p,  d,=(a—$)(a+4)- (39) 


with respect to the customary field surface. That is to say, its radius of curvature is 


S(t + dyPdy) =ERi1 + (a — 3)(1 +a)? ] + O(Rp*) (40) 


and its axial point is at a distance 


fe*utp(a— 3) + O(Rp®) (41) 
in front of that of the customary field surface. 

That the best field surface is spherical appears here as a consequence of the 
fact that the parameter © enters into the expressions (21) only as ©?, other powers 
of © being absent. 

5. Best monochromatic Schmidt camera.—Equations (39)-(41) specify the 
best field surface of a monochromatic Schmidt camera of any given focal ratio 
and angular field in terms of the plate-profile parameters a and p, which remain 
at our disposal We next show how these parameters can be chosen to minimize 
the effective »mage-radius over the given angular field. 
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In the best field surface (39) we have from (22), (34) and (38) 


(Ut, V¥)=4[p + (1 +«)0?](r? — 3)(u, v) + O7(4u8 — au, — 4v* + av); 


=| | (Ut? + V#2) du dv =p? + 2p@%X1 +2) + Or + 2] 
JJwttor<i 
+ 40a 4a+5); 


“1 
et = | 20 do (U#? 4 V#) du do 
tT) 9 j 


7 ut+y?<1 


=§lp? + p(I +) +4(1 + %)*] +4(a? — 4a+5) 
-i| (p+ a) +4(1 +a) | +4[(a—2)? +1]. 


2 


The last equation shows that the minimum value of e*, namely 
s(t +a)? +4 (45) 
is taken when 
I+a@ 
po » a= 2; (46) 


? 
“« 


that is to say, when the neutral zone of the corrector plate is at the edge of the 
aperture and when the axial image is spherically under-corrected by a wave-front 
distortion 
I ie aw 9 9 °° 9 
f.—— p8(r4 — 2r*) =} Re*p8(r* — 27°) nearly. 

In this case the field surface has radius of curvature }R(1 + 4,*) and its axial 
point is set at a distance 4fc?u* behind the marginal focus of the axial pencil. 
The effective image-radius is then 


a 2 x4)" 1/2 
EK sg3Ve8= 8 pags | OES 5 =) (47) 
seconds of arc. ‘This is 0-6 of the value for the classical Schmidt (in which a = #) 
of the same focal ratio and angular field, working in the customary field surface. 
In the case of an f/ I system working over a 15° field it amounts to 6-0 seconds of arc. 

6. The “best mean square’? Schmidt camera.—To extend the analysis to 
systems working over a whole spectrum, we recall that each value of the 
parameter p corresponds to a particular wave-length A by means of the equations 

pd = - thee : n’ =n'(d). 
Ny —1 

The first of these is sirnply the definition of p; the second expresses the refractive 
index n’ of the corrector plate in terms of the wave-length A. The value p=o 
corresponds to the “stigmatic wave-length”’ of the system, at which n’ =n and 
the angular deviations of the axial image are O( Ky’). 

The equation 


1 -_- 
e*(a, p,n; €)= * | 2040 | (US? 4+. Vi") dude (48) 


uv+rcl 


then defines the normalized mean square effective image-radius in “ p-light’’ of 
the system working over the region ¢< ¢, of an adjustable field-surface specified 
by the deformation function e(@). 

6* 
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To obtain the normalized mean square effective image-radius when the 
system is working in integrated light comprising a whole spectrum of wave- 
lengths, e* must be replaced by its weighted p-mean over the spectrum, that is 
to say by 


e*= | e*(a,p,n; cop) dp p(p) dp, (49) 
where the weighting function p(p) expresses the effective spectral density in terms 


of the parameter p. The effective image-radius in seconds of arc for the system 
is then 


Kp die*!?. (50) 
Evidently p(p) depends on the properties of the photographic emulsion used, 
as well as on the spectral characteristics of the light entering the system. 
If «(A) is the effective spectral density in terms of the wave-length A, then 


p(p) =(% — 1) di, (a) w(A(n')) 


(ny 1)45 (Fp) w(Aln+ (mp— 1)45P)) 


where n’ =n'(A) is the refractive index of the corrector plate at wave-length A 
and A(n’) is the inverse function of n'(A). p(p) is effectively independent of the 
choice of the free parameter mp, since this choice is restricted to a range of values 
of length O(,*). 

6.1. Best field surface.—Set 


L- | p(p)dp, p= =| pp(p) 4p 


and, generally, 
7 a . 
f= Z| APe(p) ap. 


Then, by (49) and (35), 


e* Z| e*(a,p,n; «)p(p) dp 


“ oe hee i 
ey + | €(O)A (O).20d0 4 | €*°(O)0 do, 
/0 /0 
where Jf (©) is defined by (36).. It follows, just as in Section 4, that the best 
field surface for given a, p, is now that for which 


«(O) = —# (0) 


7 | {(a—s) [p+ (a+ 30] 30"}p(p) dp 


2{(a—4)[p+ («+ $)07]— 307}. (54) 


That is to say, the best field surface is, in fifth-order approximation, a spherical 
surface of parameters 


C,=(a—+)p, d, =(a—)(x+4)—3 (55) 
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with respect to the customary field surface. In the special case where the light 
is monochromatic, (55) reduces to (39). 

(55) shows that the best field surfaces for a given system for different 
wave-lengths all have the same curvature. ‘hey differ from each other only 
by a focus shift and the best compromise focus depends only on the mean centre 
of the spectral intensity distribution with respect to the parameter p. When 
a=, the best field surfaces in the different wave-lengths all coincide and the 
best compromise field surface is therefore independent of the character of the 
light. 


2. Calculation of e* in the best field surface.—By (22), (34) and (54), the 
normalized deviations in the best field surface (55) are 


(US, VS) =4 p+ (1 + 2)07](7? — §)(u, v) 
+ ©?(4u® —au, —4v*+av) 
—2(p—p)(a— })(u, v) 


-(U*, VT) —2(p—py(a—4)(u, 2). 
Whence 


= | | (0 3? + VI) du dv 
: l(t m4VP*) dudv— + (p-p)a-$) || (wu +uV)dudv 


+ 4 (p—p)(a- a)? || (u? + v7) dudv 


Alp? + 2p07(1 +a) + O11 +a)?] + }O4%(a* — 4a + 5) 
+2(p—p)*(a—3), 
since, by (34) and (37), 


Z2({ i ‘ 2 ( , ‘ 

=| | (uU* +0V$) dudv = - | | (uU* +o0V*) dudv+e(O) =o. 
a). 7). 

Therefore the value of e* in the best field surface (55) is 


“ *F 
es(a, p,n) = * } 29 de) | | (US* + VI") du dv 


=i[p? + p(it+a)+h(1+«)*]+}(a—4a+5) 
+2(p—p)(a—3)*, 


and that of e* is 
e*(a,n) =+[p? + p(1+a)+ 4(1+)?]+)(a*—4a+5) 
—— \(a—$) 
+2 b(t +2)’ y | +8 (a—2)°+1] 


2 
4 
+2r2|(a—*)? +2]; 


P-()? =(p—Py 
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is the square of the radius of gyration of the spectral weighting-distribution, 
while p is its ‘centre of mass”. r, depends on the angular field radius 4, 
(it is proportional to 1/42), but not on the focal ratio of the system. 

6.3. Minimization of e*(a, n).—A design which may be expected to give 
something near to the best practical performance is obtained on minimizing 
e*(a, n) by choice of the plate-convexity constant a and the stigmatic index n; 
we shall refer to this design as the “‘ best mean square’’ Schmidt system. By (47), 
a change in n is equivalent to a change in p-origin; and equation (60) shows that 
the “‘ best” choice of p-origin is that which makes p= —(1+«)/2. In other words, 
the stigmatic index m should exceed by (m)—1)¢2(1 +«)/2 that value of n’ which 
corresponds to the p-centroid of the spectral density-distribution. Because n’ 
is very nearly a linear function of A, it is of little practical consequence whether 
the p-centroid or the A-centroid is used. We can therefore say, to a sufficient 
approximation, that the stigmatic index m of the system should exceed by 

(ty — 13" -4¢; nearly (61) 


that index m) which corresponds to light from the mean centre of the spectrum. 

(60) also shows that the “ best”’ choice of a is that which minimizes 

2 oe $\2 , 2 
[(a— 2)? +1]+2r[(a—5)* +35}. 
For this expression to be stationary, a must satisfy the equation 
\(a—2) + 4r(a— 4) =0, 
which gives 
I 
a=}+ }— . 
I+12ro 

‘This is therefore the “best” value of a. It always lies between 4 and 2; ina 
monochromat it is equal to 2, since then r, =o. 

When a, are given the above values, e#(a, 7) attains its minimum value 


* o it fo 
afi + a)" + } t ale ——s ’ 
27 : 2r- 


(62) 
which reduces to (45) when 7, =0. 

6.4. Systems with free plate profile.—Variation of the plate-convexity constant a 
and of the “‘ stigmatic index”’ » is obviously not the most general type of variation 
which can be given to the corrector plate profile, since it is always restricted by 
the condition that there is a wave-length at which the system is stigmatic in fifth- 
order approximation. ‘The question therefore arises whether a further improve- 


ment in the value of e* can be obtained by introducing quite general variations 
of the plate profile, together with general deformation of the field surface. 
We shall prove the rather unexpected result that no further improvement can 


be obtained in this way; no plate profile and no choice of field surface can give e* 
an appreciably smaller value than do the plate profile 


» I 
a=i+= 


_—_ = ~ 1 : H2 
3 I + 12r? ’ n n, + Adi, ( 3) 


and the spherical field surface (55). 

In order to be capable of optimizing the performance of the system, the general 
variation introduced into the plate profile must be of the form (#)—1) 'e*y®G(r), 
where G(r) is an even function and G(r), G’(r) are both O(1) for o<r=1. 
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The effect of such a variation in plate profile is, to a sufficient approximation, 
the addition of a term 


raph? 2 O\ rr wes _\ G(r) 
Key (5, +15) G0) —Kewru- io) (64) 


to the angular aberration components (12). The corresponding addition to the 
normalized deviations U*+iV* is 4/Ky°¢> times this, namely 
G'(r 
oo (u+1v). 


In the “free plate’ system, therefore, the equations (34) for the normalized 
deviations in an arbitrary field surface are replaced by 


U*, V* = US=G'(r) +ue(), Ves ~G'(r) + ve(9). (65) 


It follows, by the same argument as before but with U¥* + = G'(r), V+ - G'(r) 
in place of U#*, V#, that the best field surface of the new system 


e(O) = — ; | H {O)p(p) ap, (66) 


where 
KH (O)= 2 | | (uU* +0V*) dudv + . | rG'(r) dudv 
mJ. mJ) J 
“1 
= —2{(a—4)[p+(a+4)07]—307}+4 rG'(r) dr. (67) 

“0 

(66) and (67) give 

€(O) =2{(a—5)[p + («+ 3)O*] — 30°} — 41(G), (68) 


where 
“I 

KG)= | PG (dr. (69) 
“0 


Comparing (68) with (54), we see that the only etfect of the free plate variation 
on the best field surface is to move it bodily along the axis and that the values 
Us, V# of U*, V* in the new best field surface are obtained on adding 


| = G7) 416) | 0) (70) 
to the expressions (56). They are therefore given by the equation 
(Us, V3)= (UF, VI) +(u,0) | =G')—4N1G)—2(p-Pa-$)]- (72) 
To evaluate e* in the new best field surface, we first calculate 


‘ | | (Us? + Vs?) dudv 


rer 


- : | | (Ut+ VP) dude 


+ ={{ wor +vVP) [- G'(r) — 41(G) — 2(p—p)(a— | dude 


+ *|{ E G'(r) — 4I(G)—2(p—p)(a -)] (u2 + v*) du dv. 
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The first term is simply (43). The evaluation of the second term is simplified 
by observing that, by (34) and (37), 

| [ (uur +vV¥)dudv=o. (73) 
Therefore the second term can be written 


P “1 
= | [= Gut +0r9) du dv = 16[ p+ @?(1 tay) } *—8)G'(r)r dr, (74) 


7 


and (72) gives 
1 [ (U4 VS) dude 
= 5 p* + 2p0°(1 +a) + ON(1 +x)?] + $O%(a? — 4a +5) 


+ 16[p+07(1+«)] [ (r? —§)G'(r)r° dr 
0 


‘1 


+2 | E G'(r) — 41(G) —2(p—p)(a— >| P dr. 


Jo 
Hence 


» If _ , 
ef =— | 20 dO | | (Us? + V**) dudv 
0 


7 


“1 
=e +16 | p+ = J l. (7? —2)G'(r)r* dr 


$2 i; [; G'(r) -41(6) | r3 dy 


- “1 I 
~—8(p—p)(a—}) , | Fo -416) rdr 
by (59), and we obtain 
fae +16] p+ 24) |" -newra 
“0 


+2 [ GQnrdr—4 (| cme ar) }. (75) 


The problem is now to minimize e* by choice of a, m and the arbitrary even 
function G(r). 

The discussion can be shortened by noticing that p may be set equal to 
—(I+«)/2 without loss of generality provided that at the same time the arbitrary 
even function G(r) is replaced by the arbitrary even function 


h(r) =G(r)+ (p+ ~<*) (r*~ar*). (76) 


For these two changes together leave the expressions on the right of (71) unaltered 
in value. 
On making the changes, (75) becomes 


eae +24 [ wy dr—4 ( | hon ar) }. (77) 
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Now from Schwarz’s inequality 
“1 “I cl 2 
| Fi(r) dr | F3(r) dr (| F(r) F,(7) ar) 
~0 “0 0 
with 
Fi(r)=r'"h'(r), Flr) 


it follows that 


“] “1 2 
| h'*(r)rdr 4 (| h'(r)r* ar) : 
~0 “0 


with equality only if 
I / 
- h'(r) =constant. 
r 


It therefore follows from (77) that 


with equality only if 
I sa T m / 9 
* G(r) =2B- ( 4-20), 
r 2 


G(r)=A+ Br — (2 + - >) (r4—ar*). (Sr) 
(80) shows that no choice of “free’’ plate profile and no choice of field surface 
can bring e* below the value (62), while (81) shows that, if this value is to be 
attained, terms in r®, r*, ... must be excluded from G(r); that is to say, the 
plate profile variation must be of the special type considered in Sections 6.1-6. 3. 

In particular, the (a,m)-optimal solution obtained in Section 5 remains the 
“best”? monochromatic Schmidt in the mean square sense when free variation 
of the plate profile is admitted. 

7. The modified Schmidt.—\n this section we consider whether a further 
reduction in e* can be obtained by using a modified Schmidt system of the type 
discussed in (1). In this system a small amount O(R,°) of the asphericity on 
the plate is transferred to the mirrort, while the corrector plate is moved a 
small distance O( Ry?) along the axis of the system towards the mirror, so as to 
keep the coma coefficient in the fifth-order aberrations equal to zero. 

It is easy to show by the methods of (x) that the effect of these changes is to 
add a primary astigmatism term 2a(—u, v)©* to the expressions (22) for the 
normalized deviations, giving them, in the customary field surface of the original 
system, the values 


(Uf, VE) =2(u, v)[p + (1 + %)O?](2r* — a) + (u, v)a* 
+(4u® —au, — 4v° + av)O* + 2a( —u, v)O?. (82) 
The design-parameters at our disposal are now a, @ and the p-origin. Ona field 
surface which at the point © is displaced }fc?4e(@) forward from the above 
surface, 
(U*, V*) =2(u, v)[ p+ (1 +2)0?|(2r? — a) + (u, v)a@* + 2( —u, v)a0? 
+(u, v)e(O) + (4u®—au, — 4v* + av)O” 
=(U}, V3) + 2a( —u, v)O* + (u, v)e(O) (83) 





+ With a factor $(n,—1), to allow for the difference between refraction and reflection. 
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-j “A 
~| 2040 | | {[U* — 2au0? + ue(@)? 


0 Yuv'+vol 


+ [Vs + 2a? + ve(O)]?} du dv 


+ 2)" 2e(@).20d0 | { (ult + eV 9) dude 


— 


“sg 


rep 


4a dO | | (—uU* +0V*) dudv 


j, 
: 


| 84205 dO {| (u? + v*) du dv 


rl rer 
e(Q).20d0 | (u* + v*) du dv 


I 
TlO 


-l 
=et + | «(0)2(0).20d0—3a(x~$) +48 
“0 


=< 


“1 
| <(0).0d0, 
“0 

where 


# (CQ) = > | | (wU* +0V*) dudv. 
Therefore 7" 


; | e*p(p) dp 


ae a a aXe 
e* 4 | (0) A (O). 20 dO) — Fa (: _ “) + 5a" 


1 
+ | (A)Odo, 
0 
where 


ee I P 
WH (©) T | H ()p( p) dp 
= —2{(a—*)[ p+ («+ 4)0*]— 307}. (87) 
By the same argument as in Section 4, the best field surface is now that for which 
(0) = — #(0) =2{(a—$)[p + (x + 0] - 30%}; (88) 


it therefore coincides with (54) and is spherical to fifth-order accuracy. On 
this surface 


( #(O)}2@ dO — 2a ( 


me 2)" + (1 +2)°| + 2r°{( 


+}a+2a—2)?+4. 
It follows that the best mean square modified Schmidt is that for which 


a=}, a=}, (90) 
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while the p-origin is chosen to make 


By (88), (89) the best field surface of the system is 


e(@) = — 40? (92) 
and on this surface 


e*=1(r+a)?+h4+4r. (93) 


9 @ 


Comparing (93) with (62), we see that the use of a modified Schmidt does 
give an improvement in the effective image-radius 1 Ky3d2(e*)"2. However, 
the improvement is not large enough to be of very much practical consequence 
in the majority of cases. For example, in an f/1 system with c=} and r5=4 


the respective values of e* are 0-440 and 0-381; in an f,2 system with c=} and 
r- 6-02 they are 2-98 and 2-91. 

Fig. 2 shows twelve diagrams from which the performance of the optimized 
system, working in light of given spectral distribution, can be estimated at any 
given focal ratio and field size. 

The diagrams are computed at four selected off-axis distances and at three 
selected wave-lengths. In each diagram the dots show the intersections with 
the image surface of initially parallel rays entering the system through the lattice 


OSu, U=0, OI, 0-2, ...; ur=+v2 <1 
of points (u,v) filling the positive quadrant of the aperture stop. ‘The curved 
line in each case represents the image of the curved boundary of the quadrant 
on the image surface by rays of the same pencil. Because of their double 
symmetry, the structure of the ray-theoretic images can be shown more clearly 
by these ‘‘ quadrant diagrams’’ than by spot diagrams corresponding to a lattice 
of (u,v)-points filling the whole of the aperture stop. To obtain such spot 
diagrams we have merely to reflect each quadrant diagram in both coordinate 
axes. 
The parameter p’, defined by the equations 


I+a@ I 
p=- 1= — No 


2 2M,’ 
specifies the wave-length in which the image is formed. By (91), the value p’ =o 
corresponds to the wave-length for which the system is a best mean square 
monochromat. A change of 1 in the value of p’ corresponds to a change in n’ of 
id, nearly. In practice, mn’ may change by about + between the spectral lines 
Cand H. Thus p’ may range from about —1 to +1 in an f I system with a 22 
field. 

One unit of length in each quadrant diagram corresponds, by (19), to an 
angular distance }Ky°d; seconds of arc on the field surface. 

From Fig. 2 we see that, in the best mean square modified Schmidt, the 
images are smallest near the centre of the field and we infer that no considerable 
reduction in their effective radii near the edge could be obtained by rebalancing 
so as to increase the effective image-radii near the centre. 
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8. Classical and modified Schmidt cameras.—The relations between the 
classical and the modified Schmidt systems, so far as mean-square image 
optimization is concerned, can be summed up as follows: 

In the fifth-order aberration theory of the classical Schmidt, the effect of the 
central bulge is threefold. First, it shifts the best field surface towards the 
mirror; secondly it introduces primary off-axis astigmatism; thirdly it alters 
the chromatism. 

V*+ 











© =1 


. 1 Fo 


M 








Fic. 2.—Images of the best mean square modified: Schmidt at selected points of the field. The 
dots correspond to rays entering through a lattice of points filling one quadrant of the aperture 
stop; the curved lines correspond to the part of the circular boundary of the entering pencil which 
lies in this quadrant. Unit length in each diagram corresponds to a length }y°¢3f in the field 
surface. The parameter p’ specifies the wave-length, as explained in the text. 


In the result, there is a different best field surface in each wave-length; the 
best field surfaces are spherical in fifth-order approximation and all have the 
same curvature. When a= they all coincide; otherwise they differ in focal 
position and the best field surface of the system working in integrated light'is 
a compromise between them. 
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In the best field surface for any given wave-length, the mean square 
performance is best when a=2, but this value of a is not the best value for 
a system working in integrated light because colour error then spreads the 
different best field surfaces too far apart; the best value of a is always between 
3 and 2. 

In the modified Schmidt of (1), a primary astigmatism term, controlled by 
the free parameter @, is added to the fifth-order aberrations, which are otherwise 
identical with those of the classical Schmidt working in the same field surface. 
Colour error and astigmatism can now be controlled independently by choice 
of a and a. In each wave-length the best field surface is spherical and the 
fifth-order aberrations in that field surface agree with those of a classical Schmidt 
of plate-bulge parameter value a + 2a, working in its own best field surface for that 
wave-length. By choosing 

a=}, a+2a=2, (94) 
we make the different best field surfaces of the modified Schmidt all coincide 
with each other and with those of the classical Schmidt, and at the same time 
optimize the mean square performance in each of these surfaces. 

For the modified Schmidt working in the above common best field surface 
we have, by (45), 

e*=H(1 sah th+s(p4 +) 


in p-light and consequently 


arta des (B- 


in integrated light. In particular, 


e*=T(I+a)y+h+ire (97) 
if the system is designed with the amount of axial under-correction indicated by 
the equation 
I+o ° 
- {xX 
i (99) 
The Observatories, 
Cambridge : 


ea Cortor 


IQ5c 
References 
E. H. Linfoot, M.N., 109, 279-297, 1949. 


E. H. Linfoot and E. Wolf, ¥. Opt. Soc. Amer., 39, 7: 
E. H. Linfoot and P. A. Wayman, M.N., 109, 535 





INTEGRATION OF THE COWLING STELLAR MODEL 
J. G. Gardiner 
(Communicated by W. H. McCrea) 
(Received 1950 September 28) 


Summary 
A numerical integration of Cowling’s stellar model is performed to give a 
solution of known accuracy. ‘The method is described and the solution 
tabulated. 





1. Introduction.—The stellar model introduced by Cowling (1) has become 
a basic model for stars of the main sequence. The integration of the equations 
for the model were performed by Cowling (1) and have been repeated by other 
workers (2, 3). However, these computers did not claim any considerable 
accuracy and made no estimate of the accuracy attained. There is now a demand 
for a definite integration of known accuracy. In response to requests reaching 
him, Mr D. H. Sadler, Superintendent of H.M. Nautical Almanac Office, 
suggested to the writer that she should attempt to provide this. Mr Sadler also 
suggested the particular choice of dependent variables which have actually 
been used. 

The purpose of the rather high degree of accuracy which it is here sought 
to achieve is principally : 


(a) to discover whether previous solutions have been appreciably affected 


” 


by the mathematical “‘instability’’ of the equations of stellar structure; 

(6) to estimate the degree of accuracy attained by the various methods of 
integration used by other workers, especially those of C. M. and H. Bondi (3). 
These methods are not intended to yield high accuracy but it is useful if the 
results obtained can be compared with those of a method of known accuracy so 
that for other similar cases the degree of approximation can be inferred. 

It is not, of course, implied that the accuracy arrived at in the present work 
is of direct value for comparison with observational results. Neither the accuracy 
of the observations nor that of the assumptions underlying the equations would 
warrant this. 

2. Problem.—The Cowling stellar model consists of two regions: a central 
convective core and an outer region in radiative equilibrium. All the sources 
are assumed to lie within the convective zone and the polytropic index m =(y—1)>! 
has there its least value 1-5. 

The usual physical variables are defined as follows: 

r= distance from centre, a=factor of proportionality in 
M =mass contained within a radius r, Kramers’ law, 
p=density of matter at distance r a= Stefan’s constant, 
from centre, c=speed of light, 
?=total pressure at distance r from G =constant of gravitation, 
centre, A =gas constant, 
T = temperature, # =molecular weight. 
= luminosity. 
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The equations of hydrostatic equilibrium and conservation of mass are 
dP = GMp 


- 
dM 


yy =47pr-. 


In the convective zone 
pa T3 2 
which can be written p= AT?32, (3) 


whereas in the outer region in radiative equilibrium, using Kramers’ opacity law, 
d —aLp* 
= ine emia 
dr (ja7") 4ncr*T72° (4) 
Changing 7, 7, p and M to the dimensionless quantities x, y, z and w, where 
r=Xx, T=Yy, p=Zz, M= Wa, equations (1) to (4) become 
d(yz) 


dx 


dw _ 
dx 


y*?, in convective zone, 
dy kz" 


ones xylan » in outer region, 


5 
and the constants X, Y, Z and W are given by 
X=)", Y=u"*V4, Z=Au*V*, W=4rAu*yV*, 


where u=(54/8rpGA)"", V=(3xLA?/167ack)"", and k is a numerical 
constant (see below). 

The problem is to find the solution of these equations such that: at the 
centre, x=0, y and z are finite and w=0; at the interface of the two regions, 
x=X, y, %, w and their first derivatives with respect to x are continuous; and 
at the outer boundary, x=R, y and z tend to zero together (i.e. J and p are 
zero together at the boundary of the star). For any x, the continuity conditions 
of dy/dx at x= Xy give a unique numerical value of k. For the solution satisfying 
all the required conditions there is a unique x, and hence a unique value of k. 

The solution of equations (5), (6) and (7) is the solution (4) of the Emden 
equation for m=1-5. Hence a solution satisfying the required conditions is one 
that is the same as the solution of Emden’s equation up to x=xy, and which 
for x >X, satisfies equations (5), (6) and (8), and has y and z zero together for 
some x= R( >»). 

Near the outer boundary x= R, a series solution of equations (5), (6) and (8) 
can be obtained. Outer boundary values of the variables are denoted by the 
suffix s. 

Put €=R/x—1, so that =o at the boundary. Then assuming that y can be 
expanded in a power series of the form y=y, + y€+ terms in higher powers of €, 
expansions of z and w as power series in terms of € can be obtained by using 
equations (5), (6) and (8). 
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The purpose of deriving a series solution was to calculate boundary values 
of the variables (see next paragraph) and it was found that for the required 
accuracy it was sufficient to take the first two terms in y, i. e. y=y,+ Yé, and 
the terms dependent on this expression in the expansions for z and w. 

The resulting series are 


y=Yé, y,=0, 
g= RYVARWE+ 13/4, 
woe, ALM +8 — BE + Be 
where w,=7RY, 
A, =% Ri? YMA i, 


w 


It may be noted that to this degree of approximation the function zy"'* is of 
constant value. 

3. Method of integration.—The numerical integration was carried out 
using the variables x, y, w and 6, where 6=zy™4, With these variables 
equations (5), (6) and (8) become 


dy 
dx 
dw 


= x26ylsis 
’ 
dx 


dé 50 
—— GS ae 1-7k6? — WwW). 
dx aay | " ) 
The boundary condition for 6, to the same degree of approximation as 
in (9), is 


6=0,=(RY/R)!2, 


i.e. 8 tends to a constant non-zero value as x->R and yo. 

This is a unique solution separating two sets of solutions for which either 
4-+0 as y->0 or 6-> 0 as y->0, according as xX», the radius of the convective core, 
is smaller or larger than the required radius. (See Fig. 1.) 

An approximate value of the radius of the convective core was found using 
the method of integration suggested by C. M. and H. Bondi (3). This gave a 
value of x) between 1-1920 and I-1924. 

‘The procedure used for the integration was as follows: 

(a) for some x=Xxy calculate yy, 2, w, and their derivatives from Emden’s 
solutions for n=1'5 (4); 

(b) from equation (8) find & at x=xp using Vg, 3, wo, etc. ; 

(c) integrate equations (10), (11), (12) by the usual methods described in 

. . aa . . dy 
Interpolation and Allied Tables (5), i.e. assume values of y and 2h Te ate known 
7 ax 


for Xp, X1, %, (where x,=X%)+A, x,=X y+2h) and thus the differences 51(y), 


dy 
Sy), 5! [ 2h dy , o- (| 2h s , etc., and then 
wn dx dx 


; 1 dy 
(i) extrapolate a value of 6° (24 i) at x2, 
dx 
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(ii) calculate a value y, (say) for y, at x, by the equation 


dy of, : 
Ye=Nit (2n 2), + $5 (2n z), + higher terms, 


where suffixes I and 2 denote values of functions at x,, %2, 


(iii) use yy, (also wy and @, which are similarly extrapolated) to calculate 


(ene from equations (10), (11) and (12), 
dx 3 


, dy ’ dy 
(iv) difference (2n) at x, to give a revised value of 8 (ane . 
ax 7 dx)» 
: dy 
and, if necessary, to correct yp, and 2h= . 
x 


(v) calculate 


ie al ( . dy dy I . dy . dy \ 
50) =ye— n= 24 (2h), +(2h2) - rit (2n z). ~§ (20%) )| 


which gives y, (also w, and @,) ready for the next step in the integration. 

Initial values of y, w and 0, etc. found from (a), together with two other 
sets of values calculated from series solutions of equations (10), (11) and (12) 
were sufficient to start the integration. ‘The first few steps of the integration 
were performed using an interval of 0-01 in x; this was changed to 0-05 and 
later to 0-10. The integration was continued until it became clear from the 
behaviour of @ and dé/dx that the solution was deviating from the required one. 
(See Fig. 1.) 
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Fic. 1.—Diagram to show the extreme sensitivity of the solution for 0 according to the choice of Xe 


A new integration was then started using different values of x and k and the 
procedure was repeated. From two integrations on either side of the required 
solution an intermediate integration could be started. It was not necessary 
to recalculate all the initial values of y, w, 4, etc., since these could easily be 
interpolated from the initial values of the other integrations. 


7 
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By this method it was possible to find the numerical solution as far as x = 6-5. 
To finish the integration, appeal was made to the series solution (9). Using the 
known values of y, w and @ at x = 6-4 and 6-5, and assuming that 6=6, = constant, 
the values of R and Y were calculated from the expansion y= Yé. 

Then using the series for w and the values of R, Y and k, the values of w, 
and A,, were calculated. The value of 6,=(RY)!*k-'* was recalculated, but 
was found to differ from the assumed value of 6, by only a few units in the sixth 
decimal place. This discrepancy made no difference to the calculated solution 
however. 

The integration was carried out using six decimal places, but owing to the 
cumulative effect of rounding-off errors, etc., the sixth decimal is unreliable 
and is not given in the table. In general, the fifth decimal should be reliable 
especially in w and 6 for which the solution is very stable. 

4. Results.—The principal results obtained by this method of integration 
are as follows: 


Radius of convective core = 1-192296, 

k = =0:197409, 
Radius of star = 701168, 
Mass of star = 3:12610 Emden units, 
Mass of convective core =0°45759 Emden units, 
Mean density/central density = 0-02749. 


The number of decimals given corresponds to five decimals in the solutions 
and is therefore expected to be reasonably correct. 

5. Comparison.—The results are compared with those obtained by other 
authors in the following table: 


TABLE | 
Comparison with previous solutions 


R (Emden units) Teore/R = Moeore'Mstar Pmean’ Pcentrat 
Present paper 7'O117 01700 01464 00275 
T. G. Cowling (1) 7°027 0'169 O°144 0'027 
F. Hoyle and R. A. Lyttleton (2) 6-967 O'r71 0°146 0028 
C. M. and H. Bondi (3) 7026 0°170 ° 
° ° 


147 0'027 
C. M. Bondi (6) 7662 "163 


‘180 0020 

6. Boundary conditions.—While it is irrelevant to the performance of the 
computations, it is of interest to see the physical meaning of the so-called 
“instability of the solution”’. ‘This instability of the solution is, in fact, a result 
of the stability of the model with regard to its boundary conditions. 
A slight departure from the required solution when integrating outwards 
leads to a solution of the equations differing greatly from the former near 
the boundary. Conversely, this means that a very large difference in the 
boundary conditions is needed to give an appreciable difference in the interior 
of the model. ‘Tables II, III and IV illustrate the effect on the boundary 
conditions of a slight departure of @ from its true value. Tables II and 
III give the solutions for two models A and B which start to deviate from 
the Cowling model at x=6-60 and 6°75 respectively. Table IV gives a 
comparison of values of the temperature and radius for models A and B and 
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TABLE II 
Boundary values for model A 
y 2 
00164 4°78 X 10 
"0143 3°06 X IO 
‘0122 1°84 x 10 
‘O102 I1'OI X 10 
‘0082 4°97 X 10 
“0065 1°99 x 10 
"0057 "74 X10 
‘0056 *37 X10 
‘0056 ‘29 X 10~” 
"0056 's <8 
‘0056 <10 
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TaBie III 


Boundary values for model B 


y 
a 


4°99 X 107% 
"03 X 10 
‘96 x 10 
*43 X10 

x10 
6 XI0 


TABLE IV 


Comparison of y and x for a common density z= 5:0 x 107"! for Cowling model and 
models A and B 


Cowling 

model 

y (Temp.) ©*0004 005 0'0027 
x (Radius) 7°00  «. 7°04 


Model A Model B 


TABLE V 
Numerical solutions of the Cowling stellar model 
y z 

1°O 1'O 

0°99533 0°99751 
99335 “99005 
“93510 97773 
97305 96074 0208: These 
"95910 *93929 ‘ values 
‘94159 91367 "06823 - not used 
"92125 “88424 ‘ 2 in 
“89828 $5136 "155 integration 
87285 ‘31547 
"84517 ‘77699 
"81547 "73040 
*78398 “69415 "40528 "53092 
°75224 64965 "5702! ‘03877 
"72124 "60300 ‘ "74404 
"69097 55526 ‘80608 ‘84609 
"66144 "50742 0°93304 ‘94434 
*63267 *46036 1°06526 "03831 
*60469 "41484 I*19913 *12764 

0°57752 0°37143 1°3335! ‘21202 
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"46679 
59753 
*72450 
"84667 
‘96323 
"07358 
17732 
‘27417 
"36407 
“44704 
*§2321 
"59280 
"65609 
*71338 
"76504 
“81143 
"85291 
“88988 
‘92268 
95169 
‘97721 
‘99962 
‘O1QIQ 
03622 
°05097 
06368 
"07459 
"08389 
‘09179 
"09844 
10401 
10864 
11246 


*29128 
"30534 
"43416 
"49782 
55644 
‘O1017 
"65925 
*70388 
‘74434 
78089 
"81380 
"84333 
"86975 
"89332 
"91429 
‘93287 
94931 
"96380 
‘97653 
‘98768 
‘99742 
‘00590 
01324 
‘01960 
"02506 
‘02974 
03372 
03710 
03995 
04232 
04433 
"04598 
04733 
04842 
04930 
"05000 
*05055 
"05008 
‘05131 
*O5155 
05173 
"05185 
‘05194 
"05199 
"05202 
"05204 
"05205 
"05206 
"05206 
05206 
"05206 
‘05206 


o'55118 0°33058 
"52571 *29263 
“50110 "25769 
‘47738 "22588 
"45454 "19713 
"43258 ‘17135 
"41149 "14540 
*39127 “12810 
*37189 ‘11023 
35334 09458 
"33558 08093 
*31860 ‘06908 
"30236 ‘05882 
-28683 "04997 
‘27200 04235 
‘25781 "03581 
"24426 "03022 
‘23130 "02543 
‘21891 "02136 
"20707 "01789 
‘19573 ‘01495 
*18489 "01246 
“17450 ‘01035 
-16456 "570 X 107° 
"15502 ‘O71 X 10° 
-14589 ‘814 x 107% 
*13712 "759 x 10°73 
"12871 ‘879 X 10 
*12063 "145 X10 
"11286 "$35 Xx 107° 
*10540 ‘031 xX 107° 
‘09822 ‘616 10 
‘09131 275 * 
08466 ‘978 x 10-4 
07825 "729 Xx 107* 
‘07207 ‘916 x 1074 
06610 467 x 10-4 
06035 *324 x 10~* 
‘05478 ‘427 x 10~* 
"04942 ‘737 x 10-* 
"04422 210 x t0~* 
"03920 180 x 107° 
03433 5°315 x 10% 

‘02963 3°2904 x 10° 

‘02507 1914 xX 107° 
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Cowling’s model at a common density of z=5:0 x 10". For the two models A 
and B, the temperature remains approximately constant after x=7-0, while 
the density continues to approach zero. 

It should be noted that these differences in boundary conditions are 
independent of other boundary changes such as a change of opacity coefficient. 
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Summary 
A model of a red giant star is constructed having a definite composition 
and the correct boundary conditions. It is shown that a pure hydrogen 
star having the characteristics of this red giant cannot exist, but that a model 
having an outer region composed of hydrogen and an inner core of hydrogen 
and helium is self-consistent. 





1. Introduction.—The problem of the red giant stars has been considered 
recently by a number of workers including Hoyle and Lyttleton (1), Li Hen 
and Schwarzschild (2), and C, M. and H. Bondi (3); the work of these authors 
was based on a suggestion by Hoyle and Lyttleton (1). Their methods of tackling 
the problem have been to construct stellar models based on certain physical 
assumptions and to compare these models with actual observed stars. In this 
paper, a model is built up starting with observed data of mass, radius, luminosity 
and effective temperature of a typical gK star. An integration is performed 


working inwards from the boundary of the star, assuming a definite composition 


‘ 


for the ‘‘atmosphere”’ and using the accurate absorption coefficient for this 
composition. In doing this, the “atmosphere” of the star is dealt with fully 
and the star has the correct boundary conditions at the surface. ‘The 
composition for the ‘atmosphere’? was assumed to be hydrogen with an 
almost negligible quantity of metals, the latter having no effect except to give 
negative hydrogen ions in the outermost parts of the “atmosphere” and thus 
make the absorption be due to H~ there. 

As a result of integrating inwards, it became apparent that a red giant star 
composed of hydrogen only cannot exist, since the mass m(r) within a radius r 
is found not to tend to zero with 7, i.e. there is an excess mass. However, when 
allowing for a change in composition from hydrogen to a mixture of hydrogen 
and helium at some level inside the star this excess mass is accounted for and 
a self-consistent model is obtained. ‘hus, in accordance with the views of 
the previously mentioned workers, there must be a change of composition 
between the inner and outer regions of a giant star. 

The outer regions of the atmosphere, as calculated for this model, will 
probably possess convectively unstable zones. However, in order to give a 
definite model, this has been ignored, and the calculations given are based on 
the assumption of radiative equilibrium in this part of the star. 

The model is, then, assumed to consist of three regions : 

(1) a radiative envelope composed of hydrogen, 

(2) an intermediate zone in radiative equilibrium, with a composition of 
hydrogen and helium, 

(3) a central convective core, with the same composition as region (2). 


‘ 


’ 
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The outer region (1) was considered in some detail, and conditions in the 
atmosphere were investigated in the manner similar to that used by Strémgren 
in his investigation of stellar atmospheres (4). 

2. Opacity coefficient.—The value of the opacity coefficient used for the 
integration is the Rosseland mean « of the continuous absorption coefficient «,, 


IH (*r wu*e* hv 


= Toe ts U= Tp (1) 


I 
K - 4r4 , 0 K, (1 =¢ us , 


This is evaluated for hydrogen in the way suggested by Strémgren, with 
extensions of his tables where necessary. 

Following Strémgren, we write @ = °°; then, for @>0-go and in the presence 
of metal atoms, « is due to the absorption coefficient of the negative hydrogen 
ion; for 0-06<@<0-9, « is due to absorption by neutral hydrogen and negative 
hydrogen ions; and for @<0-06, « is due to absorption by neutral hydrogen. 

(a) 6>0-9. 

Let «,-— be the absorption coefficient of the negative hydrogen ion (4), then 


log Ky- = 7°21 — 0°12 + 0-700 — 3 log T'+ log p,+ log (1—xy), (2) 


where p, is the electron pressure, and x, is the degree of ionization of hydrogen, 
given by 
YH 


———, = = 048 — 13°S 5 ‘—log p.. 
(ites) O'4 13°530+ 3 log T—log p, (3) 


log 


The excess electrons are supplied by the metals. ‘These elements (1 metal atom to 
10’ H atoms) have no other effect upon the calculations than to allow the 
H opacity to operate. If these other elements are not present, the opacity is 
so low for the lowest range of temperature that it leads to impossible results for 
this type of star for pure hydrogen. 

(b) 0-06 <@<0-0. 

« is calculated according to Strémgren (4) with extensions to his tables where 
necessary. 

(c) 6<0-00. 

An approximation is used to give a simple expression for « in terms of 
temperature and density. 

x,(H) the neutral hydrogen absorption is given by (5) 


6474 me} Z° hy rT oS et e's . . 
x,(H) iy oT ag —~ 7s 1-3} = - eae Z=A= 1 for H. 
343 \ch®myk A ree te an © 2u, | 


For small 6 (large 7), u....u,; are negligible in comparison with u,, and hence 
we may write as an approximation 


6.4774me!® 


t ew 3 ith e 
—=je= —— with = —L 
u? fis 4 3ch3?myk®’ 


I 15 “oa T?u3 ute u 
so that -=— —- | —— au. 
x 4m), \ C )aQ-e” 


The lower limit in this integral is taken to be zero, since as long as x, is not 
exactly zero in the interval 0 to u,, the contribution from this part of the range 
makes very little difference to the integral, 


-!) 


6s Tf u'e~ 
= | du 
o ( 


I 
x 4ge'C),(1-e)) 
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This integral is evaluated by Unsdld (5). 


Then log «=8-865 +3 log @ per gram neutral hydrogen, 
log x =8-865 +3 log 6+ log (I —x,) per gram, 


log (I—xy) = -- log Ge + :) = log ot for highly ionized hydrogen, 


N,=number of neutral atoms in ground state, 


Ny* =number of positive ions, 
log (I — xq)== — 0°48 — 13°530 + § log T— log p,, 

and log p,=7°916+ log T+ log p. 
Hence, log x=+28-37 + 13°530 —2 log 7'+ log p, (4) 
or write log «= 28-37 — log T+ log p for small 0. (5) 

This value of « for the photoelectric opacity was taken into account 
throughout the interior of the star. As the integration proceeded inwards, the 
contribution to the opacity from another source became increasingly important. 
‘This was the electron scattering «,, which for hydrogen is«,=0°39. The Rosseland 
mean of these two opacities should have been calculated, but in order to make 
the integration simpler, a simple addition of the two factors was taken. The 
contribution to the opacity from the free-free transitions was also considered 
but was found to be negligible in comparison with the other factors. 

3. Equations of radiative equilibrium.—The outer region of the star is 


assumed to be in a state of radiative equilibrium and the conditions in the star 
are governed by the following equations: 


ede : Pai 
hydrostatic equilibrium ect Pe) =— = P, (6) 


P 1M 
conservation of mass = =4nr"p, (7) 


energy transport =f = “pF, (8) 
€ 


where p, =radiation pressure, 
p,=gas pressure, 
F = flux of radiation, where 47 F =acT! = L/r*, 
L=luminosity, 
T, = effective temperature, 
M = mass contained within a radius r. 
Define the optical depth r by 
dr= — Kpdr. 
Then, dp,=}aTidr, 
and since p,=4aT’, 
it follows that 4aT* = }aT4> + constant, 
i. €. T* = T3(3 + 37), 
the standard equation for radiative equilibrium (8). 
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4. Integration of the equations.—{a) Through any sufficiently small depth 


J. . 
of the atmosphere g = =z May be considered to be constant. 


Equation (6) gives d(p,+ r+ Po) a 
pdr 


dp, a. Bi 


dr K 


dr . Pg 
d(logp,) loge(g«—(a/4)T3)- 





or 


This is written as 
I '  p,d(log p,) 
t (log —--- | 544 II 
(log py) = gloge! (1k—(a 4g)T4)* (11) 
The integration of this equation was performed by the standard central difference 
formula 
ra+nh 


; | f(x)dx = $f(a)+fla+h)+ ... +fla+(n—1)h) +43 } f(a+nh) ++ 50 f(a) 


“a 


~— 48’f(a+nh) 


and the interval h was taken small enough for the second differences to be 
negligible. 
The calculations were penne for a large proportion of the wees 


since, for the most part, the term — —T! was negligible in comparison with = = 


This meant that the radiation pressure had become important in ceunguitinn 
with the gas pressure, with the further result that the density decreased instead 
of increasing with decreasing radius. ‘This term was taken into account in the 
integration, but further investigation showed that if it had been neglected the 
difference would have been very slight. 

(5) The above method of integration was suitable until r began to vary too 
rapidly for small increases in log p,, and g was no longer constant due to changes 
in r?. At this stage, s=1000, different variables were introduced to make the 
numerical work easier. From equation (6) 


Upyt+P,) ___ GMp 


dr ry 


Since, for r> 1000, p, and —~ Pe are negligible, write p, =p, 
dr 


GM, 
then ? =— f. 


dr r° 


Writing dr=—kp 2 dr (r,=radius of star), 


also d= - 2 dé6, 
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dp .|GM 


then from (12) and (13) Toe aoe, 
dr on 
16GMé@é 1 

cnt = = 40, 
16G0 M(r) 
3 1, KO p, 

1664G loge ¢ M(r) 
log p= log p, - — — —— 
gp Z Po 3 2 | Kp, 


d (log p) = do, 


dé. 


dr yp dp 
tr,  GMp 
Tdp ; n kp 
=— kGMnp +X), wp a’ 


I I ‘my(l+xy) T 
Yo i fF k M(r) d (log p) . (15) 


Also M(r)= | 4nr°pdr. (16) 


Equations (14) and (15) were integrated step by step by the central difference 
method. At each step of the integration of (14), an extrapolated value of g was 
used, then the value of r was calculated from (15) and was used to recalculate g 
and to repeat the step in the integration if this was necessary. The value of M 
calculated from (16) varies very slowly and for the main part of this integration 
was found to be constant. 

(c) This integration was used up to 6=0-05, when a different method of 
integration was adopted, this being an adaption of the Bondis’ method of 
integration (6). 

From equation (8) 


dp, _ 


adr c 





I 
Hence ~ 
F 


kp F 


aed T ol 
ot la 
dT 0 Lk y + Kg) 
dr tonacr?T® ° 
_ aT Be 
dr PT PTS 
Using the Bondi variables 
S= ce al <a a i d(logp) 
GM? d(log T) 


and equation T= 


aT = 
dr A\pdr_ p? dr}’ 
dT dp (: p z) 


a" ha Se 
; ‘ t-. ae. 
using (18) GMs N= 77 —* 
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th vf af Y f 

Writing MN= e + S , B= (*)' ” BG 1 C= (*)' CG", 
He fe 

dM , dN _ a 

mM t Ww 74+55e)— -G+45e) 7, 


then 


C’grs —I 
where w= (1+ Fa) : 
SdN _ S—(4N—1)—0(5°5N—1) 
N dS — 1+2S—40 


w satisfies the relation 
w® _ 4nB fp \® i 
(i1-—w)*®= CB (5) Oso.” 
Equations (6) and (7) become 





Hence 





SdQ_ S+N-O 

OdS 1+2S—40° 
Equations (19) and (20) were integrated numerically, starting with values of 
S, Q and N given by the physical variables p, p, etc. at a convenient point in 
the integration (at S=o0-0015). At each stage in the integration a value of w 
was extrapolated and then checked; w varied very slowly and it was rarely 
necessary to repeat an integration. 

To satisfy the normal physical conditions pertaining to the centre of a star, 

the variables S, O and N should satisfy the conditions (6), (7). 


(20) 


(a) O>a, S>wasr->o. 
(6) OS} as r-o. 


(c) N increases from N =0-21 to 0-40, this latter point indicating where the 
convective zone would start. 
(At S=0-0015, N is given by N= = ( a) fdr, “(=) pl) 
G \ ph Mp® | C\ ps edt 
Conditions (a), (6) and (c) were not satisfied, however, since after an initial 
increase S and N reached maximum values and then started decreasing again. 
Physically this showed that although ro, M did not tend to zero with r and there 
was, in effect, a surplus mass. 

5. Change of composition.—The star, as considered so far, has a composition 
of hydrogen. If, at some point just before the maximum values of S and N are 
reached, the composition is assumed to change to a mixture of hydrogen and 
helium, there would be a discontinuity in QO and N to values Q, and \N, which 
would then yield the required solution. At the discontinuity, i.e. near the 
maximum values of S and N, the opacity was due entirely to electron scattering. 
Therefore the solution for the central region, required to fit onto the outer 
hydrogen region, was the electron opacity solution given by C. M. Bondi (7). 
The procedure, then, was to calculate values of QO and NV from this solution at 
some value of S. The discontinuity in QO gave the discontinuity in the molecular 
weight and the new composition of X per cent H: Y per cent He, since, if suffixes 
I and 2 denote values of the variables for the outer and inner regions (1), 


OQ: Pe ea tee (m=bme= 
O. py Me 2/3+5X°\U2 7"? 345X/7 
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From the tables of electron opacity solution, the values of p, p, M, T and r for 
the edge of the convective core, and the values of p, p and 7 for the centre of 
the star were calculated. ‘The central temperature depended on the position 
of the discontinuity, so, by adjusting the level of discontinuity, central 
temperatures suitable for the energy-producing cycles could be obtained. Three 
models with differing central temperatures are given: one with a central 
temperature of 2-5 x 10’ deg.; one with 7'=1-g x 10? deg., this being a recent 
estimate of the central temperature of a star (g) ; and the last gives J = 1-4 x 10’ deg., 
a central temperature suggested by H. Bondi (in press). 

Providing that the underlying physical assumptions are correct, it has been 
shown that the giant K star consists of a central core, comparable to a main- 
sequence star, and a highly distended atmosphere of low density. ‘The 
composition of the central region is accurately given by the discontinuity of Q 
and is found to be approximately the same wherever the discontinuity is 
assumed to occur, i.e. a central ratio of hydrogen: helium=1:3. Probably the 
presence of other elements such as carbon, nitrogen and oxygen will slightly 
alter the values of 7, p, p, etc. at the centre of the star, but it does not seem 
probable from the numerical calculations that it would have a large effect. 

The author wishes to express her indebtedness to Professor W. H. McCrea 
for suggesting this work, and for his invaluable advice and criticism. Her 
thanks are also due to the Department of Scientific and Industrial Research 
for the grant under which this work was done. 


A Model of a gK Star 
Mass =4:90 x 10% g, 
Luminosity =7:21 xo* ergs./sec. 
Radius =2:-77 x10! cm. 
Effective temperature = 3580 deg. K. 


‘TABLE I 
Model1 Model2 Model 3 
Hydrogen zone (1) 
Inner radius r, 3°73 x 10!® 5:00 x 10!" 6°67 x 10!® cm 
‘Temperature at 7, 8-0 x10® 6:0 x10% 4:4 x 10° deg 


Mass of zone (1) 2°46 x 10" 2-29x 10*7 2:19 x 10" g. 


Hydrogen-helium zone (2) (Radiative equilibrium) 
Molecular weight 4. 0:96 93 o'"9! 
Ratio X: Y I : 3°29 22°85 Ea°S7 


Inner radius r, 2°15 x 10! 2*84x 10! 3:67x 10!9 cm. 


Temperature at r. 16 x107 1I'2 x10? og x10? deg. 


Mass of zone (2) 1°52x 10" 1°65 x10" 1°72x 10% g. 
Hydrogen-helium core (3) (convective core) 


Central temperature 2°§ x10" ‘9 x10" 1°4 x10’ deg. 


Mass of core 0:98 x 10" 1-02xK 10" 1:05 x10" g. 
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TABLE II 


T, r, p and p for the atmosphere of a gK star (M constant) 


Method (a) 

T log p log p- logk (1—Xyg) logp 
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Method (6) 
bp, X7=1, M=4'96 x 10* g. 
log p log K log fe) 
‘51 2°11 8-95 
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Method (c) 
Using Bondi’s method of integration. p.—4p, xy 

log p log p 
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A NEW ECLIPSING BINARY OF VERY LONG PERIOD 
D. F. K. O'Connell, S.7. 
(Received 1950 September 11) 


Summary 
The light curve of a recently discovered eclipsing binary, based on about 
1100 photographic observations at Riverview, is discussed. The eclipse 
lasts more than 17 years. The period is probably of the order of 200 years. 





‘The variable star C.P.D. —60° 3278 was discovered by the writer in 1949 
on Riverview plates. 

Comparison stars.—The comparison stars are listed in Table I, with 
their C.P.D. numbers, positions for 1900, the spectra given in Potsdam 
Publication, No. 89 (P) and in Harvard Annals, 112 (H), and their photographic 
magnitudes. ‘These magnitudes were determined from comparisons with 
Miss Wright’s nearby sequence for SV Centauri.* 

TABLE | 
Comparison Stars 
1900 Spectrum 
R.A. Dec. rs 


h m 8 ° ’ 

60°3278 11 38 20 60 11 Bs: B3 = 

59 37 11 39 17 59 54 FS FS8 9°04 

59 3729 11 3801 60 07 Ar B8 9°45 

59 3799 11 42 18 60 03 Ao Ao 9°68 

60 3310 II 40 33 -60 08 Ao Ao 9°75 

60 3289 11 38 52 60 17 Bg Bg 10°32 
g 60 3281 II 38 30 60 14 A2: Ao 10°47 


Pg. Mag. 


Observations.—Estimates were made on 1090 plates taken at Riverview 
Observatory from 1931 to 1950. Eighty-one of these plates were taken with 
the G camera (Ross—Grubb-—Parsons), the remainder with the R and P cameras 
(Zeiss Astrotriplets). No systematic differences were found between plates 
taken with the different cameras. ‘The brightness of the star remained constant 
from 1931 until about J.D. 2,430,000, when a minimum commenced, which is 
still continuing. ‘The observations, grouped into mean points, are given in 
Table II. During constant light the means are usually of twenty estimates. 
During minimum about ten estimates are combined into a mean. At maximum 
the longest interval covered by a mean point is about 260 days, most of the 
intervals being much shorter than this. During minimum the longest interval 
included in a mean is about 100 days. 

Light curve.—The mean points are plotted as circles in Fig. 1, the area of 
the circles being proportional to the weight, which was taken as equal to the 
number of estimates in the mean point. Unfortunately, fewer plates were taken 
during the war years. Nevertheless, the observations are sufficient to determine 
the shape of the curve reasonably well. The form of the light curve appears to 

* Frances W. Wright, Harvard Ann., 89, 10, 1937. 
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TABLE II 


Mean Points 





J.D. Mag. No. | J.D. Mag. No. .D. Mag. No. 
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Fic. 1.—Light curve of C.P.D. —€0° 3278, 1931-1950. Circles denote mean points, the 
area of the circles being proportional to the weight. 
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leave no reasonable doubt that the minimum is due to the eclipse of one star by 
another. Variation of the R Coronae Borealis type is excluded because of (1) the 
spectrum of the star, B3, (2) the slow descent to minimum and (3) the small 
range. ‘Ihe curve drawn in the figure shows a partial eclipse, with mid-eclipse 
at J.D. 2,433,200. It is possible that a total phase begins about J.D. 2,433,000, 
but it is considered that the most recent observations show a slight brightening 
of the star and indicate that the eclipse is partial. Observations over the next 
year or two will settle the point. 

If the eclipse is partial, its duration is about 17 or 18 years. If it is total, the 
duration is longer still. ‘The period must be very long indeed. In order to carry 
the search for another minimum as far back as can be done here, estimates were 
made on Sydney Astrographic plates and on Franklin-Adams charts. The 
estimates and their dates are as follows: 

J.D. Pg. Mag. 
mi 
2412541 9:60 Sydney plate. 
12865 9°62 os 
18799 9°45 F.—Adams chart. 
22008 9:60 


” ” 


24207 9°78: Sydney plate 
of plate. 


Doubtful estimate : star near edge 


These observations, combined with the Riverview ones, seem to exclude the 
possibility of another eclipse, unless perhaps a very shallow secondary one, 
between 1893, the date of the earliest Sydney plate, and the present eclipse. 
An examination of the Harvard plates would remove all doubt about the matter. 
The latest possible date for the middle of an eclipse of the same depth as the 
present one is 1884, so that the period is certainly longer than 65 years. It is 
probably very much longer. ‘There are seven other eclipsing binaries known 
with periods longer than 1000 days. ‘They are given in the following list, with 
their periods (P) and ratio of period to duration of minimum (D). 

P P/D 

S Doradus 14670 d. 13 

« Aurigae 9883 12°5 

VV Cephei 7430 15 

V 38x Scorpii 6475 

V 532 Ophiuchi >6000 

V 383 Scorpii 4900 

BM Scorpii 1760 
From analogy with these long-period binaries it is probable that the period of 
the new binary is more than ten times the duration of minimum, so that the 
period is probably of the order of 200 years. 

Even though the period of the binary is unknown, and will probably remain 
unknown for a very long time to come, it is possible to deduce some information 
about the component stars from the shape of the light curve. 

The magnitude at maximum, determined from 682 estimates, is 
g™-549 + 0™-002 (p.e.). The average of the measured deviations of the individual 
estimates from the mean is 0™-051. ‘Thus the probable error of a single 
estimate is +0™-043, which is about the probable error usually found by the 

8 
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writer for estimates of eclipsing binaries. The probable error of the mean of 
twenty estimates is then +0™-o10. It is clear that the magnitude at maximum 
is sensibly constant, so that the stars may be taken as spherical. The magnitude 
at minimum is 10™-25. The loss of intensity at mid-eclipse is 0-475. 

The following discussion is based on the assumption that the eclipse is partial. 

Phases during minimum were calculated for periods of 100, 200 and 300 years, 
and the Russell-Shapley y function determined from the observed light curve. 
x (R, %&, $) is 1°82, 1°85, 1°85 for periods of 100, 200 and 300 years respectively. 
It appears that the length of the period has comparatively little influence on this 
function. We may take x (k, %, $) as 1°84. ‘The uncertainty in this value is 
considered to be not so great as 0:05. 

On the assumption of uniformly illuminated disks the range of possible 
values of k, the ratio of the radii, and a», the fraction of the disk obscured at 
mid-eclipse, is from k=1, %=0°52 +0-08 to %=1, R=0-68+0-04. The limits 
indicated for k and % correspond to a range of + 0-05 in y (R, a, 4). 

For complete darkening to the limb the limiting values of k and a» obtained 
from the light curve are k=1, %=0°53, and a=1I, k=0-54. 

It follows that in any case the two stars are fairly comparable in size. In this 
discussion the orbit is taken as circular. This is probably incorrect, but it will 
probably be a very long time before the orbital eccentricity is known. 

The calculated duration of eclipse is about 6300 days for the uniform, and 
6700 days for the darkened, solution, a little over 17 and 18 years respectively. 
The longest eclipse previously known, 1130 days, and the longest period, 
40 years, are those of S Doradus. ‘Thus the new binary is of particular interest. 
It is hoped that spectrographic observations will be made ir the near future, 
while the variable is still at, or near, minimum. It would be very helpful to have 
the spectrum of the companion and the relative luminosity of the stars and, of 
course, radial velocity measures. 

With such a long period it is likely that the orbital eccentricity is high. 
Aitken * finds that binaries with periods of 100 to 150 years have an average 
eccentricity of 0-530, and those with periods over 150 years one of 0-615. ‘Thus 
the eccentricity in this case is probably greater than 0-5. It is even possible that 
the components may be separable telescopically at some part of the period. The 
images of the variable on Sydney Astrographic plates were examined in order 
to detect any sign of duplicity or elongation, but the images appeared normal 
on all the plates. 

Results 
m 
Maximum 9°55. 
Minimum 10°25 about J.D. 2,433,200. 
Duration of minimum Atleast 17 years. 
Period >65 years, probably of the order of 200 years. 


Acknowledgment.—Thanks are due to Mr H. W. Wood, Government 
Astronomer of N.S.W., for permission to use the Sydney Astrographic plates. 


Riverview College Observatory, 
Riverview, N.S.W. : 
1950 September 1. 


*R, G, Aitken, The Binary Stars, 207, 1935. 





ENERGY LEVELS AND TRANSITION PROBABILITIES IN 
p? AND p* CONFIGURATIONS 


R. H. Garstang * 
(Received 1951 January 11) 


Summary 


The energies of the levels of the configurations p* and p‘ are calculated, 
including the Hamiltonian terms arising from the interaction of the spin of 
one electron with the orbit of another and the interaction of the spins of pairs 
of electrons. The spin interaction parameters and unperturbed energies of 
the levels are determined empirically from observational spectroscopic data. 
Comparisons are made with parameters computed from self-consistent field 
wave functions where available. Revised transition probabilities are com- 
puted for the forbidden lines within the configuration p* for the ions C I, 
NII, OI, FTV and Ne V and within the configuration p‘ for the ions O I, 
F Il and Ne Ill. 





1. Introduction.—Previous calculations by Pasternack (1) and by Shortley, 
Aller, Baker and Menzel (2) + of the transition probabilities of forbidden lines 
in the configurations p" were based on the theory of intermediate coupling, 
taking as the perturbation the spin interaction of an electron with its own orbit. 
The results they obtained led to an intensity ratio of the 43727 pair of lines 
4S—*D of OII in disagreement with observations of gaseous nebulae. ‘This 
discrepancy was explained by Aller, Ufford and Van Vleck (3), who showed that 
in the configuration p* the effect of the spin-other-orbit and spin-spin interactions 
between pairs of electrons is large and comparable with the second-order effect 
of the conventional spin-orbit interaction, the latter alone having been included 
by Pasternack and SABM. ‘The former terms are particularly important for 
the p* configuration because the first-order effect of spin-orbit interaction is 
zero for half-closed shells. Similar calculations were performed for NI by 
Ufford and Gilmour (4). In the p? and p* configurations the spin-orbit term is 
greater than the spin-other-orbit and spin-spin terms. ‘The latter are, however, 
partly responsible for the departures from the Landé interval rule in the 
*P term, and it is of interest to use these departures to 1etermine empirically the 
parameters involved and to recompute the transition probabilities of some ions 
of astrophysical importance. 

2. Energy levels in the p* configuration.—We use the perturbation theory 
of Quantum Mechanics to determine the energy levels. ‘The matrix of spin-orbit 
interaction has been given by Condon and Shortley (5) and the matrices of 
spin-other-orbit and spin-spin interactions by Marvin (6). (His M® is our 7.) 
Combining these matrices we form and solve the secular determinant. 

Let ¢ be the atomic parameter involved in spin-orbit interaction and 7 that 
for spin-other-orbit and spin-spin interactions. Let the energies of the 
unperturbed |S, 'D and *P terms be E(1S), E('D) and E(?P) respectively. Write 

* Isaac Newton Student, University of Cambridge. 
+ Henceforth abbreviated to ‘‘ SABM ”’. 





116 R. H. Garstang, Energy levels and Vol. 111 


(SP)=E('S)—E(@P) and (DP)=E(*D)—E(®?P). Then the solutions of the 

secular determinant are 

a ae a 
2(DP) 4(DP)?~ 5(DP)’ 
e " c3 ‘ 14ly 

2.DP) 4(DP)  5(DP)’ 

E(?P,) = E@P) — 46 + 8n, 


E(°P,) = E@P)+ 3o— Fy 





E(?D,) = E(D) + 





E(?P,) = ECP)—£+8n— ces, + cope — Ce 


en. tt. 4 A 
(SP) (SP) (SP) 
Terms in ¢4, ¢*y and 7? are neglected. 

The triplet intervals are 


E(1S,) = E('S) + 








Ce vf “3 ss 14ly 
(DP) 4(DP?  5(DP)’ 
en ee ee 
(SP) (SP) (SP)° 





E(®P,)— E@P,) =¢-2n- 2 


E@P,) — E(°P,) = 36+ 


Departures from the Landé ratio r=2, where 


,_ ECP,)-ECP)) , 
~ EGP,)= BCP s)' " 


are caused partly by terms containing the second and higher powers of ¢ and 
partly by terms containing ». In particular the term “y is often important, 
especially if ¢ is small. 

We shall treat E(}S), E(}D) and E(*®P) as independent parameters to be 
fitted to the observations rather than assuming the relation 6(SP)=15(DP) 
between the inter-multiplet separations predicted by the Slater theory. It is 
well known that the predicted separations do not agree with the observations 
due to the effect of configuration interaction. It is hoped that by using the 
observed rather than the predicted separations at least a portion of the effect of 
configuration interaction will have been taken inte account. ‘here are thus five 
independent equations (1) in the five unknowns {, y, E('S), E('D) and E(?P) 
which can be solved by successive approximation when we insert in the equations 
the observed energy levels. All the observations used in this paper were taken 
from the compilation of Charlotte E. Moore (7). 

Table I contains the results of these calculations for all the atoms with 
p* outer configurations for which the necessary data are available in Miss Moore’s 
tables. 

For hydrogen-like wave functions (" varies linearly with the atomic 
number Z for a given isoelectronic sequence and for non-hydrogenic wave 
functions this is still a very good approximation if a suitable screening constant 
is used. Examination of the values of {4 for the 2p* sequence in Table I shows 
that the values of { for Si 1X and P X are somewhat lower than would be expected 
on the assumption of a linear relation between ('* and Z. Whether this is 
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due to a failure of the theory or to a falling off of the accuracy of the observations 
cannot be decided on the basis of present data. 

It can be shown (Appendix to this paper) that for hydrogen-like wave 
functions 7" varies linearly with Z, for the atoms of an isoelectronic sequence. 
Using the values of » for the 2p? sequence in Table I we get a very nearly linear 
relation except that for PX the value of 7 is much too small. (‘This is clear from 
an inspection of Table I.) No explanation of this can be offered. 


TABLE | 
Empirical Parameters for p* Configurations 
[units : cm.~"] 
n EP) E(:D) ECS) 


97° ‘Ol 89°5 15316 326086 
222° 209 20270 43182 
436° 3°06 419 25237 53536 
789 5°6 759 30283 63880 

1304 8:3 1284 35333 74227 
2055 2052 41405 35544 
3080 18 3140 460570 95950 
4440 2 4610 52740 107380 
5700 13T 6140 59000 118890 


32° 0°38 29°7 10194 21648 
I 
I 


148°5 0°03 151° 6296° 15391° 
314°9 0°37 321 3377 21567 
549°7 ohio) 572 11306 27140 
886 o3 932 13734 32502 
* Observations of 1S and 1D uncertain. 
+ No apparent reason for this low value. See text. 

The values of ¢'4 are linearly related to Z for the 3p? sequence. The 
values of » obtained are very small and uncertain and no reliance should be 
placed on them. 

3. Energy levels in the p* configuration—The matrix of spin-orbit 
interaction is the same as for p* except that —¢ is substituted for ¢. We require 
diagonal matrix elements for spin-other-orbit and spin-spin interactions. 
These have not been calculated previously. ‘The state M=2 of the level *P, 
has as its unperturbed wave function the single Slater determinant (in Condon 
and Shortley’s notation) 

Jrtr-ot 1+]. 
The corresponding diagonal matrix element splits up into the sum of 
24 integrals involving spin-other-orbit interactions and 24 integrals involving 
spin-spin interactions. A typical integral is of the form 

| (1*1-| H|1*r-). 
Each of these integrals is calculated from the formulae of Marvin (6) 
(equations (21) and (22)). If H' is the term in the Hamiltonian due to 
spin-other-orbit interactions and H" that due to spin-spin interactions, then 
we find on adding the 24 integrals in each case 

(P,|H"|*P,) = — 30, 


(4 
(2P,|H™|%P3) = —27. 
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Kramers (8) gave the formulae 
(LSJ|H'| LSJ) =const. [J(J +1)— L(L+1)-—S(S+1)], 
(LSJ |H™| LSJ) =const. {{L(L + 1) + S(S+1)-J(J +1)? 
~([L(L+1)+ S(S+1)-J + 1)]-$ L(L 4+ 1)S(S + 1)} 
for the diagonal matrix elements, where the “ constants”’ 
but may depend on L and S. 
From these we find 
(‘D, | H*|*Dz) =(*S|H*|"So) phe. 
(‘D,|H™|1D,) =(*So|H™|1S,) =o, 


are independent of J 


and using our matrix elements (4) 
(°P,|H"|*P,) = 3) 
(°P,|H"|®P,) =6n, 
(P,|H"|®P,) =27, 
(°P,|H™|8P,) = — 40. 
A study of the detailed computations for p? shows that the effect of the 
non-diagonal elements of H! and H"! is very small. We shall neglect such 
elements in our work on the p* configuration. 

We can now form and solve the secular determinant. Retaining powers of ¢ 
up to the fourth and neglecting ¢y and 7’, we find for the energy levels, in the 
same notation as before, 
hg ha 


E@P,) = E@P)— 36-5 1- spp) + [DPE * BPP’ 





E(?P,) = ECP) + 30+ 5, 

au a | a 
(SP) (SPY * (SP)? 
c rs ha as cf 
2(DP) 4(DP?  &(DP)’ 
2¢° 20° 204 


E('Sp) = E('S) + (SP) + (SP) bi (SP) 


E@P,) =ECAPP)+¢+2n- 








E('D,) = E(!D) + 








The triplet intervals are 





2(DP) 4(DP) 8(DP)*’ 


(6) 


v2 ‘ v4 
E(@P,) —E@P,) =C+ “ n+ ~ by G | 
| 


ECP,)- ECP,)=}{-3)- 25, - 725, + 25 
v7 oat | =$§6-—-3)-— 5 ‘ona Y he a Y =r 
. ye ST (SP) (SPP (SPF 
Fitting these to the observational data for all atoms in the p‘ configuration, for 
which data are available in Miss Moore’s tables (7), we obtain the results in 
Table II. 

We find that ("4 is nearly linear in both the 2p' and 3p' sequences; for the 
higher members of the sequences the values of ¢'* are slightly less than would 
be expected on the basis of a linear law. 73 is nearly linear in Z in the 2p' 
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sequence up to Si VII, the values for P VIII and SIX are too low. ‘The values 
for the lower members of the 3p* sequence are probably of the correct order of 
magnitude but no reliance should be placed on the actual numerical values. 
‘They depend rather critically on the data used. For 7iVII and V VIII the 
calculated » was negative. It is not possible to decide whether observational 


TABLE II 


Empirical Parameters for p* Configuration 
{units : cm.~?] 


n E@P) ECD) EC'S) 


x 


15867 33791 
20870 44914 
25834 55733 
31099 66748 
36310 77638 
41520 88520 
46850 99480 
53180 110410 
57530 121400 


146°6 

320°4 

607 
1041 
1673 
2565 
3760 
535° 
737° 


33 


26 


muumt ONU ES NN 
© Oo Oot 
wMnwin 
° 
w ui 


ttt 


384° " 9231 22168 
670 35 11633 279005 
1065 57 13970 33195 
1588 16399 38408 
2254 2 15694 43597 
3093 4°1 787 21172 48806 
4180 t 23750 54040 
5452 t : 26560 59570 

* Observations of 4S and '!D uncertain. 

+ Calculations give a negative value. 

t Low values. 

§ Estimate by Edlén. 


uncertainties or approximations in the theory are responsible for these various 
discrepancies. It should be remembered that we have used a method of fitting 
theory and observation which leaves no degree of freedom, the solution of five 
equations in five unknowns. 

4. Comparison with theory and discussion.—The parameters ¢ and 7 are 
given, in the usual notation, by the expressions 


(7) 


_ e*h* tel ss y—-3 p2 D2 
: Torimia | |, r, >P*(np, r.)P?(np, r,) dr, drs. (8) 
The values of { and 7» have been obtained from self-consistent field wave functions, 
including exchange, for CI (9), OI and OIII (10) and NelIII (x1) and are 
listed in Table III together with the corresponding empirical values. Aller, 
Ufford and Van Vleck (12) pointed out that an extra term 
om - 
4m P mc On (nxt . Ss; 


t 
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has to be added to the Hamiltonian because spin-other-orbit interaction has 
already been included in the term 
- gia (* 2) 24.0 
Sa*mec? \r dr} 5 
to the extent that the field from the other electron can be regarded as central. 
o,,() is the fraction of charge of an np electron lying interior to radiusr. ‘The 
extra term increases the effective value of { by 47. The empirical values of ¢ 
(Tables I and II) obtained from the observational data should be compared with 
¢'=€+ 4 calculated from equations (7) and (8). This comparison is given in 
Table III. 
Tase III 
Empirical Theoretical 

lon 4 ” 4 ” Ss 

cr 32°6 0°38 42°8 0°40 44°4 

O lll 222°0 1°90 268-7 1°94 276°5 

Oo! 146°6 1°33 184°6 I'ls 189°2 

Ne Ill 607 4°3 690 37 705 


The agreement of theoretical and empirical values of 7 is very good. Although 
the amount of comparison data is not large we seem to have succeeded in 
isolating the portion of the departures from the Landé ratio due to the 
spin-other-orbit and spin-spin interactions. ‘The agreement is not so good 
for ¢. ‘This must be due to the approximations used in the theory. It is 
thought that the empirical values of ¢ are likely to be fairly close to the true 
ones and they will be used in subsequent work. We may compare our values 
with those used by Pasternack (1) and by Robinson and Shortley (13). For 
the 2p” configurations our parameters are greater, for the 2p* smaller than those 
used by these authors. For the 3p and 3p' the differences are very small. 

It is of interest to examine the contributions to the multiplet splitting in °P 
due to the various terms in equations (2) and (6). ‘Table IV gives the ratio (3) 
for a selection of atoms in the p* configurations due to : 

(1) terms in { only (Landé ratio), 

(2) terms in { and higher powers of ¢, 

(3) terms in ¢ and », and 

(4) the complete formulae (the observational ratio). 


The importance of the contribution of » for low values of Z in the 2p4 
configuration is clearly shown. Similar results hold for the p? configurations. 
For very highly ionized atoms in the 2p? and 2p‘ configurations and for atoms 
whose principal quantum number is greater than 2 the effect of the 
spin-other-orbit and spin-spin interactions is quite small. 


TABLE IV 


(1) (2) (3) (4) 
O! 2°00 2°05 2°28 2°33 
Ne Ill 00 2°12 2°21 2°35 
ALVI ‘00 2°34 2°12 2°50 
SI ‘00 2°20 2°04 2°24 


Araki (14) obtained a general formula for multiplet separations including 
spin-orbit, spin-other-orbit and spin-spin interactions and used it to analyse 
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a large number of terms of multiplicity not less than four. He neglected second- 
order effects of ¢. For terms involving d electrons the contribution of ¢ is 
quite large and almost certainly swamps contributions of 7. ‘The writer has 
made an attempt to analyse the d? configuration in 77111 and found it impossible 
to isolate the contribution of », no reasonable values of ¢ and » would give a 
satisfactory fit of the observed levels. It would seem that the neglect of @° at 
least partially vitiates Araki’s work. It seems impossible to say more than this, 
because his method introduces different parameters from ours and as he does not 
determine any of his numerical coefficients explicitly for any atoms, and 
wave functions are not available in any case for atoms with d shells, a comparison 
with theory is not possible. The internal agreement he obtains for most cases 
is, however, remarkable. 

5. Calculation of transition probabilities.-The method has been described 
by SABM, who gave formulae for the line strengths in terms of the coefficients 
in the expressions of the perturbed wave functions in terms of the Russell- Saunders 
wave functions. Using the usual perturbation theory these coefficients have 
been obtained from the spin interaction matrix. 

Let ©, ¥ denote respectively unperturbed and perturbed wave functions. 
We write 

YD.) = aD('D,) + bO(F P,), 
Y(3P,) = —b@('D,) + aM(? P,), 
P(? ») = “D(? Py) + dD("S5), 
YS») = —d®(? Py) + cO('S,) 
and, since there is only one level with J =1, 
Y(*P,) =0(°P,). 
Introducing the abbreviations 
C 7 c 
i= db = ~ . c= 
(DP)’ (DP) (SP)’ 
we find that for the p* configuration 
2 3 
a=1-}0°—-}6..., 
b= y/2[40+ 107-104 :44 
c=I-—a*+207+.,., 
d ‘2[o —o*—20% +7+...], 
where terms in #*, o', 64, or, d* and 7* are neglected. ‘hese expressions 
differ from those of SABM by the inclusion of terms in » (¢ and 7) and by the 
use of the observed values of (DP) and (SP) instead of assuming (DP)=6F, 
and (SP)=15F, of the Slater theory. We prefer to use the observational 
values because the actual ratio (DP):(.SP) does not agree with theory. 
Similarly we find for the p* configuration 
a=I—}07+16+.,.., 
b= /2[— 36+ 407+ 409+ ...], 
¢=I—o"—20°+..., 


d= /2[—a—07 + 207+ . 
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where #4 and o4 are neglected. ‘Terms in ¢ and 7 do not appear in these formulae 
because we have neglected the non-diagonal elements of the spin-other-orbit 
and spin-spin matrices. 

Having obtained a, 6, c and d we use SABM formulae (7) and (8) to obtain 
the line strengths and thence the transition probabilities for the various 
transitions in the configurations. 

Detailed numerical results have been computed for the p* configuration in 
Cl, NU, OUT, FIV and NeV. In the calculation of the quadrupole proba- 


bilities estimates of the integral s, - | 1°P*(2p)dr are required. ‘These have 
5 Jo 

been obtained for CI and OIII directly from the self-consistent field wave 

functions (g, 10). Using these values a relation of the form s>'*=a%(Z—, 

was determined (Z is the atomic number and § a screening constant) and was 

used to estimate s, for the other ions. 

Numerical results have been computed for the p* configuration in Ol, FII 
and NelIll. ‘The values of s, for OI and Nell were calculated from self- 
consistent field wave functions (10, 11) and that for FII estimated in 
the same way as for the p* configuration. 

The transition probabilities are listed in ‘Tables V and VI. ‘The tables also 
include the adopted values of s,,. 


TABLE V 


Transition Probabilities for the 2p? Configuration 


jsec, 3 

Transition ‘Type Ci NU O lll FIV 

1D.'S, e ‘08 6 
P'S, e 0°16 
3P,—S, m 034 
‘PD, m "0030 
e "0°94 

‘PDs, m ‘00103 

e ‘0°14 

*P."D, e 0°42 
‘PP, m 0°74 

e ‘0229 

iP. P; e o13 
‘PP, m ‘0°21 
Sq (atomic units) . 0°735 


50 
‘o'19 
"00260 


~ 


"0023 
I 
‘095 
0713 
"034 
‘0121 
0°64 
0°78 
‘0°12 
‘0°50 
‘o%2¢ 
"323 


- ON 


©00000000000 
eooocgco0o0oe0dodgj0dtd Oo 
02020000000% 


° 
0000000000 


oo00000000 0+ 0 


electric quadrupole. 
magnetic dipole, 
o'o"r denotes r « 10> "~1, where r is the first significant figure. 

Appreciable differences are found between our results and those of 
Pasternack (1). ‘These are due partly to the adoption of new estimates for s,, 
rather smaller than those used by Pasternack, partly to the re-evaluation of ¢ 
and the isolation of », and partly to the use of the observational values of (DP) 
and (SP) instead of the theoretical ones. ‘The estimates of s, should be more 
nearly correct than those used by Pasternack, who based his on screening constants ; 
wave functions are normally too diffuse and the true values are probably still 
smaller than the ones we have used. Examination of the numerical magnitudes 
of the quantities involved shows that the inclusion of spin-other-orbit and spin- 
spin interaction terms in the expression for the perturbed wave functions makes 
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very little ditference to the resulting transition probabilities; the main effect of 
the terms in 7 arises from the new determination of ¢ from the energy levels. 
It has not yet been possible to extend these calculations to include the etfect 
of superposition of configurations. If this could be done we would probably 
obtain an explanation of the difference between the empirical and theoretical 
values of (DP) and (SP) and perhaps improve the agreement of the parameter ¢ 
obtained from the theory with that from the observations. Part of the latter 
difference is, however, probably due to the wave functions used. It seems 
impossible to predict the effect of superposition of configurations on the tran- 
sition probabilities but it is hoped that part of the effect has been taken into account 
by the use of empirical values of (DP) and (SP). 

TaBLe VI P 

Transition Probabilities for the 2p* Configuration 
[sec.~'] 
Transition 'l'ype O! Fil Ne Ill 

28 2°38 
‘0016 0°005! 
49 2°2 
‘044 0°20 
‘0'96 
0138 
0113 


"a 
>°0°3 
‘078 
0°0009 
‘of24 
D°0022 
"0°32 


‘O° 18 
‘0°18 
‘0&9 


“oiis 
O*17 
"0° 103 
-o'8$Q 
0!"13 
Sq (atomic units) *79O 


‘o*2 3 
"433 


2 
° 
° 
° 
° 
° 
° 
0'o° 11 00541 
° 
° 
° 
° 
° 


e=electric quadrupole. 
m-==miagnetic dipole. 
o'o"r denotes r x 10~"~—1, where r is the first significant figure. 


Acknowledgments.—'Vhis work was stimulated by the paper of Aller, Utford 
and Van Vleck (3). ‘he writer would like to thank Professor D. R. Hartree 
and Dr B. Jetfreys for advice and encouragement during the investigation. He 
would also like to thank Dr S. F. Boys for referring him to the work of Araki 
and Dr D. Layzer for a suggestion in connection with the paper of Marvin, 
Thanks are due to the University of Cambridge and to Caius College for the award 
of Studentships, and to the Department of Scientific and Industrial Research 
for a grant. 

APPENDIX 

When the hydrogen-like radial wave function for nuclear charge Z is substituted 

in the formula (8) for » we get, after some simplification, 


ik en el ge ee gE 
aimee eT) wat |, SO La de | atts (LTP de, 





(nl) = 


n+l 

where L**! is the (2/+1)th derivative of the Laguerre polynomial of order 
(n+), n the principal and / the azimuthal quantum numbers of the electron and 
a the radius of the first Bohr orbit of hydrogen. We see at once that for a given 
(nl) electron » is proportional to Z*. ‘The integral above has not been evaluated 
in the general case. In particular cases it is easy to substitute the appropriate 
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Laguerre polynomial and carry out the integrations. The results of this 
calculation, expressed in cm. |, are, for three simple cases, 


n (2p) =0°0133 2%, 
(3p) =0-00113 23, 
9 (3d) =0-00124 Z?. 


If we fit expressions of the form 7 =A(Z—c)* to the observations in ‘lables I 
and II we find that for the 2p” configuration 7 =0-0197 (Z — 3°37)* and for the 
2p' configuration » =0-0117 (Z — 3°17)’. 


Gonville and Caius College, 
Cambridge : 
1951 January 10. 
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